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

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.

Differential loss of embryonic globin genes during the radiation of placental mammals

The differential gain and loss of genes from homologous gene families represents an important source of functional variation among the genomes of different species. Differences in gene content between species are primarily attributable to lineage-specific gene gains via duplication and lineage-specific losses via deletion or inactivation. Here, we use a comparative genomic approach to investigate this process of gene turnover in the beta-globin gene family of placental mammals. By analyzing genomic sequence data from representatives of each of the main superordinal clades of placental mammals, we were able to reconstruct pathways of gene family evolution during the basal radiation of this physiologically and morphologically diverse vertebrate group. Our analysis revealed that an initial expansion of the nonadult portion of the beta-globin gene cluster in the ancestor of placental mammals was followed by the differential loss and retention of ancestral gene lineages, thereby generating variation in the complement of embryonic globin genes among contemporary species. The sorting of epsilon-, gamma-, and eta-globin gene lineages among the basal clades of placental mammals has produced species differences in the functional types of hemoglobin isoforms that can be synthesized during the course of embryonic development.

Platypus globin genes and flanking loci suggest a new insertional model for beta-globin evolution in birds and mammals

BACKGROUND

Vertebrate alpha (alpha)- and beta (beta)-globin gene families exemplify the way in which genomes evolve to produce functional complexity. From tandem duplication of a single globin locus, the alpha- and beta-globin clusters expanded, and then were separated onto different chromosomes. The previous finding of a fossil beta-globin gene (omega) in the marsupial alpha-cluster, however, suggested that duplication of the alpha-beta cluster onto two chromosomes, followed by lineage-specific gene loss and duplication, produced paralogous alpha- and beta-globin clusters in birds and mammals. Here we analyse genomic data from an egg-laying monotreme mammal, the platypus (Ornithorhynchus anatinus), to explore haemoglobin evolution at the stem of the mammalian radiation.

RESULTS

The platypus alpha-globin cluster (chromosome 21) contains embryonic and adult alpha- globin genes, a beta-like omega-globin gene, and the GBY globin gene with homology to cytoglobin, arranged as 5'-zeta-zeta'-alphaD-alpha3-alpha2-alpha1-omega-GBY-3'. The platypus beta-globin cluster (chromosome 2) contains single embryonic and adult globin genes arranged as 5'-epsilon-beta-3'. Surprisingly, all of these globin genes were expressed in some adult tissues. Comparison of flanking sequences revealed that all jawed vertebrate alpha-globin clusters are flanked by MPG-C16orf35 and LUC7L, whereas all bird and mammal beta-globin clusters are embedded in olfactory genes. Thus, the mammalian alpha- and beta-globin clusters are orthologous to the bird alpha- and beta-globin clusters respectively.

CONCLUSION

We propose that alpha- and beta-globin clusters evolved from an ancient MPG-C16orf35-alpha-beta-GBY-LUC7L arrangement 410 million years ago. A copy of the original beta (represented by omega in marsupials and monotremes) was inserted into an array of olfactory genes before the amniote radiation (>315 million years ago), then duplicated and diverged to form orthologous clusters of beta-globin genes with different expression profiles in different lineages.

Genomic evidence for independent origins of beta-like globin genes in monotremes and therian mammals

Phylogenetic reconstructions of the beta-globin gene family in vertebrates have revealed that developmentally regulated systems of hemoglobin synthesis have been reinvented multiple times in independent lineages. For example, the functional differentiation of embryonic and adult beta-like globin genes occurred independently in birds and mammals. In both taxa, the embryonic beta-globin gene is exclusively expressed in primitive erythroid cells derived from the yolk sac. However, the "epsilon-globin" gene in birds is not orthologous to the epsilon-globin gene in mammals, because they are independently derived from lineage-specific duplications of a proto beta-globin gene. Here, we report evidence that the early and late expressed beta-like globin genes in monotremes and therian mammals (marsupials and placental mammals) are the products of independent duplications of a proto beta-globin gene in each of these two lineages. Results of our analysis of genomic sequence data from a large number of vertebrate taxa, including sequence from the recently completed platypus genome, reveal that the epsilon- and beta-globin genes of therian mammals arose via duplication of a proto beta-globin gene after the therian/monotreme split. Our analysis of genomic sequence from the platypus also revealed the presence of a duplicate pair of beta-like globin genes that originated via duplication of a proto beta-globin gene in the monotreme lineage. This discovery provides evidence that, in different lineages of mammals, descendent copies of the same proto beta-globin gene may have been independently neofunctionalized to perform physiological tasks associated with oxygen uptake and storage during embryonic development.

Phylogenetic comparisons suggest that distance from the locus control region guides developmental expression of primate beta-type globin genes

Phylogenetic inferences drawn from comparative data on mammalian beta-globin gene clusters indicate that the ancestral primate cluster contained a locus control region (LCR) and five paralogously related beta-type globin loci (5'-LCR-epsilon-gamma-psieta-delta-beta-3'), with epsilon and gamma expressed solely during embryonic life. A gamma locus tandem duplication (5'-gamma(1)-gamma(2)-3') triggered gamma's evolution toward fetal expression but by a different trajectory in platyrrhines (New World monkeys) than in catarrhines (Old World monkeys and apes, including humans). In platyrrhine (e.g., Cebus) fetuses, gamma(1) at the ancestral distance from epsilon is down-regulated, whereas gamma(2) at increased distance is up-regulated. Catarrhine gamma(1) and gamma(2) acquired longer distances from epsilon (14 and 19 kb, respectively), and both are up-regulated throughout fetal life with gamma(1)'s expression predominating over gamma(2)'s. On enlarging the platyrrhine expression data, we find Aotus gamma is embryonic, Alouatta gamma is inactive at term, and in Callithrix, gamma(1) is down-regulated fetally, whereas gamma(2) is up-regulated. Of eight mammalian taxa now represented per taxon by embryonic, fetal, and postnatal beta-type globin gene expression data, four taxa are primates, and data for three of these primates are from this laboratory. Our results support a model in which a short distance (<10 kb) between epsilon and the adjacent gamma is a plesiomorphic character that allows the LCR to drive embryonic expression of both genes, whereas a longer distance (>10 kb) impedes embryonic activation of the downstream gene.

Gene conversion and functional divergence in the beta-globin gene family

Different models of gene family evolution have been proposed to explain the mechanism whereby gene copies created by gene duplications are maintained and diverge in function. Ohta proposed a model which predicts a burst of nonsynonymous substitutions following gene duplication and the preservation of duplicates through positive selection. An alternative model, the duplication-degeneration-complementation (DDC) model, does not explicitly require the action of positive Darwinian selection for the maintenance of duplicated gene copies, although purifying selection is assumed to continue to act on both copies. A potential outcome of the DDC model is heterogeneity in purifying selection among the gene copies, due to partitioning of subfunctions which complement each other. By using the d(N)/ d(S) (omega) rate ratio to measure selection pressure, we can distinguish between these two very different evolutionary scenarios. In this study we investigated these scenarios in the beta-globin family of genes, a textbook example of evolution by gene duplication. We assembled a comprehensive dataset of 72 vertebrate beta-globin sequences. The estimated phylogeny suggested multiple gene duplication and gene conversion events. By using different programs to detect recombination, we confirmed several cases of gene conversion and detected two new cases. We tested evolutionary scenarios derived from Ohta's model and the DDC model by examining selective pressures along lineages in a phylogeny of beta-globin genes in eutherian mammals. We did not find significant evidence for an increase in the omega ratio following major duplication events in this family. However, one exception to this pattern was the duplication of gamma-globin in simian primates, after which a few sites were identified to be under positive selection. Overall, our results suggest that following gene duplications, paralogous copies of beta-globin genes evolved under a nonepisodic process of functional divergence.

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.

Sequencing and mapping hemoglobin gene clusters in the Australian model dasyurid marsupial Sminthopsis macroura

Comparing globin genes and their flanking sequences across many species has allowed globin gene evolution to be reconstructed in great detail. Marsupial globin sequences have proved to be of exceptional significance. A previous finding of a beta(beta)-like omega(omega) gene in the alpha(alpha) cluster in the tammar wallaby suggested that the alpha and beta cluster evolved via genome duplication and loss rather than tandem duplication. To confirm and extend this important finding we isolated and sequenced BACs containing the alpha and beta loci from the distantly related Australian marsupial Sminthopsis macroura. We report that the alpha gene lies in the same BAC as the beta-like omega gene, implying that the alpha-omega juxtaposition is likely to be conserved in all marsupials. The LUC7L gene was found 3' of the S. macroura alpha locus, a gene order shared with humans but not mouse, chicken or fugu. Sequencing a BAC contig that contained the S. macroura beta globin and epsilon globin loci showed that the globin cluster is flanked by olfactory genes, demonstrating a gene arrangement conserved for over 180 MY. Analysis of the region 5' to the S. macroura epsilon (epsilon) globin gene revealed a region similar to the eutherian LCR, containing sequences and potential transcription factor binding sites with homology to eutherian hypersensitive sites 1 to 5. FISH mapping of BACs containing S. macroura alpha and beta globin genes located the beta globin cluster on chromosome 3q and the alpha locus close to the centromere on 1q, resolving contradictory map locations obtained by previous radioactive in situ hybridization.

Linkage of the beta-like omega-globin gene to alpha-like globin genes in an Australian marsupial supports the chromosome duplication model for separation of globin gene clusters

The structure, function, and evolutionary history of globin genes have been the subject of extensive investigation over a period of more than 40 years, yet new globin genes with highly specialized functions are still being discovered and much remains uncertain about their evolutionary history. Here we investigate the molecular evolution of the beta-globin gene family in a marsupial species, the tammar wallaby, Macropus eugenii. We report the complete DNA sequences of two beta-like globin genes and show by phylogenetic analyses that one of these genes is orthologous to embryonically expressed epsilon-globin genes of marsupials and eutherians and the other is orthologous to adult expressed beta-globin genes of marsupials and eutherians. We show that the tammar wallaby contains a third functional beta-like globin gene, omega-globin, which forms part of the alpha-globin gene cluster. The position of omega-globin on the 3' side of the alpha-globin cluster and its ancient phylogenetic history fit the criteria, originally proposed by Jeffreys et al. (1980), of a "fossil" beta-globin gene and suggest that an ancient chromosome or genome duplication preceded the evolution of unlinked clusters of alpha- and beta-globin genes in mammals and avians. In eutherian mammals, such as humans and mice, omega-globin has been silenced or translocated away from the alpha-globin locus, while in marsupials omega-globin is coordinately expressed with the adult alpha-globin gene just prior to birth to produce a functional hemoglobin (alpha2 omega2).

Hematopoietic gene expression profile in zebrafish kidney marrow

The zebrafish kidney marrow is considered to be the organ of definitive hematopoiesis, analogous to the mammalian bone marrow. We have sequenced 26,143 ESTs and isolated 304 cDNAs with putative full-length ORF from a zebrafish kidney marrow cDNA library. The ESTs formed 7,742 assemblies, representing both previously identified zebrafish ESTs (56%) and recently discovered zebrafish ESTs (44%). About 30% of these EST assemblies have orthologues in humans, including 1,282 disease-associated genes in the Online Mendelian Inheritance in Man (OMIM) database. Comparison of the effective and regulatory molecules related to erythroid functions across species suggests a good conservation from zebrafish to human. Interestingly, both embryonic and adult zebrafish globin genes showed higher homology to the human embryonic globin genes than to the human fetal/adult ones, consistent with evo-devo correlation hypothesis. In addition, conservation of a whole set of transcription factors involved in globin gene switch suggests the regulatory network for such remodeling mechanism existed before the divergence of the teleost and the ancestor of mammals. We also carried out whole-mount mRNA in situ hybridization assays for 493 cDNAs and identified 80 genes (16%) with tissue-specific expression during the first five days of zebrafish development. Twenty-six of these genes were specifically expressed in hematopoietic or vascular tissues, including three previously unidentified zebrafish genes: coro1a, nephrosin, and dab2. Our results indicate that conserved genetic programs regulate vertebrate hematopoiesis and vasculogenesis, and support the role of the zebrafish as an important animal model for studying both normal development and the molecular pathogenesis of human blood diseases.

An orphaned mammalian beta-globin gene of ancient evolutionary origin

Mammals possess multiple, closely linked beta-globin genes that differ in the timing of their expression during development. These genes have been thought to be derived from a single ancestral gene, by duplication events that occurred after the separation of the mammals and birds. We report the isolation and characterization of an atypical beta-like globin gene (omega-globin) in marsupials that appears to be more closely related to avian beta-globin genes than to other mammalian beta-globin genes, including those previously identified in marsupials. Phylogenetic analyses indicate that omega-globin evolved from an ancient gene duplication event that occurred before the divergence of mammals and birds. Furthermore, we show that omega-globin is unlinked to the previously characterized beta-globin gene cluster of marsupials, making this the first report of an orphaned beta-like globin gene expressed in a vertebrate.

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.

Embryonic globins of the marsupial the Tammar Wallaby (Macropus eugenii): bird like and mammal like

1. Tammar Wallaby embryonic blood has been shown to have three alpha-like and two beta-like globin chains in its four haemoglobin components and partial sequences of several chains have been determined. 2. The major embryonic beta-like chain (epsilon) is similar to other mammalian embryonic beta-like chains on the basis of sequencing its first 60 amino acids. 3. There is another embryonic beta-like chain present in one haemoglobin component. It was designated omega and, in its first 54 amino acids, it has features that are more like avian globins than mammalian globins. 4. The one alpha-like embryonic globin sequenced has mammalian rather than avian characteristics. 5. A provisional phylogenetic tree of beta-like globins has been determined. The Tammar epsilon-globin forms a monophyletic group with marsupial and other mammalian embryonic globins; the omega-globin forms a monophyletic group with bird adult and embryonic globins.

Sequences of the primate epsilon-globin gene: implications for systematics of the marmosets and other New World primates

Sequences of the epsilon-globin gene were determined for five species of marmosets, along with approximately 2 kb of 5' flanking sequence. An analysis of these data, compared with those of other primates strongly supports the classification of Callithrix jacchus and C. geoffroyi into the jacchus group, and C. argentata and C. mauesi into the argentata group. The pygmy marmoset, formerly identified as Cebuella pygmaea joined strongly to the argentata group, indicating that without the pygmy marmoset the genus Callithrix would be paraphyletic. Our data support recent studies which indicate that C. pygmaea should be included in the genus Callithrix. Relationships among other primates were as indicated by previous studies of epsilon-globin sequences. Divergence times were estimated according to a local molecular clock. These calculations indicated the divergence of C. mauesi and C. argentata to be approximately 1.6-1.9 Myr (million years ago), and the most recent common ancestor of the marmosets to be between 4.5 and 4.7 Myr. The latter estimate corresponds well to the date of 4.6 Myr calculated from an independent data set.

Characterization of the embryonic globin chains of the marsupial Tammar wallaby, Macropus eugenii

The embryonic hemoglobins of the marsupial Tammar wallaby (Macropus eugenii) are known to aggregate, which was shown by the finding that the Hill coefficient, h, was greater than 4.0 in the upper part of the oxygen equilibrium curve. Here, we have undertaken a detailed primary structure analysis of the Tammar wallaby pouch young hemoglobin complement, which we hoped might provide clues into the residues that cause aggregation and a high embryonic h. The Tammar wallaby embryonic hemoglobin complement is principally four major hemoglobins each with a different isoelectric point. Two early expressed hemoglobins contain the same embryonic beta-like chain, epsilon (epsilon), but two separate alpha-like chains, termed zeta and zeta prime (zeta and zeta') both of which are N-terminally blocked. The later two expressed hemoglobins contain the same adult alpha-chain, but different beta-like chains. The latest expressed hemoglobin contains the same beta-like chain, epsilon, as the two early expressed forms, but the third expressed hemoglobin contains a unique beta-like chain which we have termed omega (omega). A protein database similarity search using the first 54 N-terminal amino acids of the omega-chain showed a range of sequence identities of 57-72% to all known mammalian beta-like chains, including the other marsupial epsilon-chains. The closest identity, reflected by both the highest percentage identity and Smith-Waterman score, was with the embryonic beta-chains of the aves. While the primary structures of the hemoglobins reported here do not explain the low hemoglobin-oxygen affinity in embryonic marsupial blood, the finding of the similarity with the bird globin-like sequence with one of the marsupial chains has implications on mammalian globin evolution. How many other marsupials and placental mammals are harboring a bird-like globin in their embryos?