Distinctive Patterns of Evolution of the δ-Globin Gene (HBD) in Primates

Ana Moleirinho; Alexandra M. Lopes; Susana Seixas; Ramiro Morales-Hojas; Maria J. Prata; António Amorim

doi:10.1371/journal.pone.0123365

Abstract

In most vertebrates, hemoglobin (Hb) is a heterotetramer composed of two dissimilar globin chains, which change during development according to the patterns of expression of α- and β-globin family members. In placental mammals, the β-globin cluster includes three early-expressed genes, ε(HBE)-γ(HBG)-ψβ(HBBP1), and the late expressed genes, δ (HBD) and β (HBB). While HBB encodes the major adult β-globin chain, HBD is weakly expressed or totally silent. Paradoxically, in human populations HBD shows high levels of conservation typical of genes under strong evolutionary constraints, possibly due to a regulatory role in the fetal-to-adult switch unique of Anthropoid primates. In this study, we have performed a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a set of diverse mammalian taxa, focusing on the evolution and functional divergence of HBD in primates. Our analysis revealed that anthropoids are an exception to a general pattern of concerted evolution in placental mammals, showing a high level of sequence conservation at HBD, less frequent and shorter gene conversion events. Moreover, this lineage is unique in the retention of a functional GATA-1 motif, known to be involved in the control of the developmental expression of the β-like globin genes. We further show that not only the mode but also the rate of evolution of the δ-globin gene in higher primates are strictly associated with the fetal/adult β-cluster developmental switch. To gain further insight into the possible functional constraints that have been shaping the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, as shown by different ω values between apes and Old World Monkeys + New World Monkeys (0.06 versus 0.43, respectively), these estimates corroborated a constrained evolution for HBD in Anthropoid lineages, which is unlikely to be related to protein function. Collectively, these findings suggest that sequence change at the δ-globin gene has been under strong selective constraints over 65 Myr of primate evolution, likely due to a regulatory role in ontogenic switches of gene expression.

Citation: Moleirinho A, Lopes AM, Seixas S, Morales-Hojas R, Prata MJ, Amorim A (2015) Distinctive Patterns of Evolution of the δ-Globin Gene (HBD) in Primates. PLoS ONE 10(4): e0123365. https://doi.org/10.1371/journal.pone.0123365

Academic Editor: Francesc Calafell, Universitat Pompeu Fabra, SPAIN

Received: June 6, 2014; Accepted: March 2, 2015; Published: April 8, 2015

Copyright: © 2015 Moleirinho et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work was supported by the Portuguese Foundation for Science and Technology (FCT) (PTDC/SAU-MET/110323/2009). AM and AML are supported by fellowships from FCT (SFRH / BD / 73508 / 2010 and SFRH / BPD / 73366 / 2010, respectively). SS is supported by POPH-QREN – Promotion of scientific employment, by the European Social Fund, and by national funds of the Ministry of Education and Science. IPATIMUP is an Associate Laboratory of the Portuguese Ministry of Education and Science and is partially supported by FCT. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Hemoglobin (Hb), found in the circulating red blood cells of all vertebrates, is the major oxygen-transporting molecule, playing a key role in the cellular aerobic metabolism [27]. In mammals, Hb is a heterotetramer composed of two α-like and two β-like globin chains that are differentially expressed during development, such that functionally distinct Hb isoforms are synthesized in embryonic and adult erythroid cells [27–29].

These globin chains are encoded by members of the α- and β-globin gene families, which arose via tandem duplication of an ancestral, single-copy globin gene approximately 450–500 Mya, in the common ancestor of jawed vertebrates [14,23,25,35,78]. The two paralogous gene families exhibit a number of significant differences in gene content among jawed vertebrate taxa. These differences are especially pronounced in the case of the β-globin cluster, in which distinct repertoires of mammalian β-like globin genes originated by independent lineage-specific duplications followed by functional divergence [33–35,52–54,57,78,79]. In both monotremes and marsupials, the β-globin gene cluster contains a single pair of genes, the early expressed ε-globin and the late expressed β-globin [53]. In contrast, within the eutherian stem, further tandem duplications gave rise to a cluster of five β-like globin genes, containing early-expressed genes, located at the 5’ end of the cluster ε-(HBE)-γ(HBG)-ψβ(HBBP1), and late expressed genes, δ (HBD) and β (HBB), at the 3’ end, consistent with the orientation in contemporary species [26,31,53]. The fine tuning of the level and timing of expression of each of these genes relies on interactions with the locus control region (LCR), located from approximately 6 to 18 kb upstream of HBE [4,11,84].

Over the course of eutherian evolution the structure of the β-globin gene cluster has been dynamic and the late-expressed HBD and HBB paralogs have experienced different evolutionary fates. In the majority of mammals, the adult form of Hb (α₂β₂) contains similar β-chain subunits which are encoded by one or more copies of the HBB gene [52]. Contrastingly, the δ-globin gene, although present in almost all eutherian species examined to date, is frequently pseudogenized [19,24,26,30,54]. In a few species, a transcriptionally active but weakly expressed copy of the δ-globin gene was maintained, encoding the δ-globin chain of the minor fraction of the adult Hb (α₂δ₂), known as HbA_2, which is thus assumed to be physiologically irrelevant [45,46,76]. In fact, the δ-globin chain is absent in Old World Monkeys (OWM) [45,46] and ranges from 1% concentration in hominoids [10] to 40% in the galago [80], reaching 6% in New World Monkeys (NWM) [74] and 18% in tarsiers [40]. Surprisingly, in some eutherians HBD shows a level of sequence conservation typical of genes under strong evolutionary constraints [91]. In humans for example, HBD was found to have lower nucleotide diversity than HBB, suggesting that purifying selection has shaped the evolutionary history of HBD [49,87] through an unrecognized role not associated with oxygen transport [46,74,86].

The involvement of HBD in the fetal/adult Hb switch was proposed decades ago [5,55] and since then some studies have provided evidence supporting this hypothesis [6,49,67]. In fact, the fetal to adult Hb switch of anthropoid primates is unique. Furthermore, while both Anthropoids and Prosimians possess a γ-globin gene, its switch after birth only takes place in the major anthropoid branch, the catarrhines, occurring earlier in NWM, whereas in Prosimians it is only expressed at the embryonic stage [38]. Therefore, phylogenetic and comparative genomic analysis across placental mammals with distinct repertoires of β-like genes and corresponding expression programs should provide clues to the evolution and putative functional divergence of the δ-globin gene.

However, the evolutionary history of the eutherian HBD is quite complex due to unusually frequent sequence exchanges through extensive gene conversion and unequal recombination with its neighbor, β-globin [19,31,40,54,80], resulting in extensive sequence homogenization and hampering the assignment of orthologous relationships among HBD and HBB genes. Indeed, HBD was initially thought to be the result of a recent duplication in primate evolution, approximately 40 MY [18], but recently an older origin has been proposed [30,53]. Under such scenario of controversy, a revisit of the evolutionary history of the adult δ-globin gene can help elucidate the origin of this gene family, which in spite of many efforts is not yet fully understood [31,37,40,45,46,53,60,74,80].

Here, we perform a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a wide range of mammalian taxa, with a special focus on primates. Our results further document reticulation in the topology of the evolutionary history of δ-globin gene, demonstrating that it has behaved as an evolutionary palimpsest, with repeated and partially overlapping β / δ sequence transfers obscuring orthology. Additionally, we show that the δ-globin gene is highly conserved in Anthropoids, with a particularly strong signal of purifying selection in Great Apes. Sequence conservation at this locus is unlikely related to protein function and may reflect mutational constraints on regulatory regions involved in the fetal-to-adult developmental expression switch of the β-globin cluster.

Materials and Methods

DNA Sequence data and gene identification

To obtain DNA sequences spanning the entire HBD and HBB genes, we used Blat queries to interrogate the genome assemblies of several mammalian species available in the UCSC Genome Browser website (http://genome-euro.ucsc.edu/index.html). Whenever the HBB sequence was available, we used it to identify its paralogous sequences (the first hit); alternatively we used the human HBB and HBD sequences to identify their corresponding orthologs. Due to a history of concerted evolution, in some cases high sequence identity between HBD and HBB produced ambiguous results. In these cases the genome coordinates were used to distinguish between the two genes, since the order of the β-like globin genes has been maintained throughout mammalian evolution. Additional genomic data was obtained either from the High Throughput Genomic Sequences database (HTSG), Trace Archives or by direct sequencing to complete sequence gaps and to include further mammalian species. Detailed information on sampling is listed in S1 Table. Following all these steps of manual curation, our sample included 29 sequences from the three major subclasses of mammals: 1 Prototheria (Ornithorhynchus anatinus), 1 Methateria (Monodelphis domestica) and 27 Eutheria, including representatives of the following superordinal groups (1 Xenarthra, 7 Laurasiatheria and 19 Euarchontoglires). We also included one avian species (Gallus gallus) as outgroup.

Evolutionary analysis

The reconstruction of the evolutionary history of the HBD and HBB genes was carried out across the mammalian phylogenetic tree with the chicken sequence as outgroup. The identified coding sequences were translated into the corresponding protein and aligned using the Expresso mode of T-Coffee in order to take into account any structural information of the protein available in the databases [2]. Each of the non-coding sequences, including the 5’ and 3’ untranslated regions (UTRs) plus 2 introns, were aligned separately as nucleotide sequences using the accurate mode of T-Coffee [51]. Phylogenetic reconstruction was performed using Maximum Likelihood (ML) and Bayesian Inference (BI) as the optimality criteria. Trees were estimated using the complete gene sequence and including only the coding sequence (CDS). The partition strategies and models of evolution implemented in these analyses were identified using PartitionFinder v 1.1.1 [41]. PartitionFinder was run specifying the raxml and mrbayes models of evolution for computational capacity reasons.

ML analyses were performed with RAxML [75] run in the CIPRES Science Gateway [48]. Phylogenetic analysis of the complete gene was run with the following partitions: 1) 5’ UTR + Exon 1, 2 and 3 + Intron 1; and 2) Intron 2 + 3’ UTR. The GTR+G+I model was implemented for each of these. ML analysis of the CDS was performed with three partitions corresponding to the three codon positions. The GTR+G model was implemented for each partition. The resulting trees were evaluated with 1000 bootstrap replicates.

BI analyses were run in MrBayes 3.2 [64] in the CIPRES Science Gateway. The complete gene analysis was performed with three partitions: 1) 5’ UTR + Exon 1, 2 and 3 + Intron 1; 2) Intron 2; and 3) 3’ UTR. The implemented models of evolution were the K80+G, GTR+G+I and HKY+G+I for the first, second and third partitions, respectively. BI with the CDS sequence was performed with the data matrix partitioned according to codon position; the models of evolution implemented were the K80+G for the first codon position and the GTR+G for the second and third positions. Two independent runs of 10 million (complete gene analysis) and 5 million (CDS analysis) generations with 8 chains each (7 heated and one cold) were set up. Trees were sampled every 200 (complete gene) or 100 (CDS) generation and the first 12500 trees (25% of the sample) were discarded as burn-in. Convergence and burn-in were assessed using Tracer 1.6 [61], MCMC Trace Analysis Package http://tree.bio.ed.ac.uk/software/tracer/).

Additionally, a third phylogenetic analysis was performed using models of codon substitution with the coding sequence only as input. This method allows us to take into consideration any potential divergence in codon usage across the mammalian lineage and differences in selective pressure among gene copies, which included pseudogenised copies. Analyses were run using CodonPhyML [21] with the GY CF3x4 [22] model of codon substitution and specifying the M3 model of selective pressure [90]. An initial tree was obtained using BioNJ + GYECMK07 and the topology was searched using the subtree pruning and regrafting (SPR) heuristic search. Branch support was obtained using the approximate Likelihood Ration Test (aLRT) as implemented in the software. To identify potential recombination events in our primate data set we used the recombination detection package (RDP3) [43], applying a set of seven statistical methods, which includes RDP [42], GENECONV [56], Bootscan [44], Maxchi [72], Chimaera [59], SiSscan [20] and 3Seq [8]. Briefly, two of these are phylogenetic methods, which infer recombination when different parts of the genome result in discordant topologies (RDP and Bootscan), while the other five are nucleotide substitution methods, which examine the sequences either for a significant clustering of substitutions or for a fit to an expected statistical distribution: MaxChi, Chimaera, Geneconv, 3Seq and SiScan. The latter primarily uses genetic similarity estimates but also takes some phylogenetic information into account. To uncover both species specific and ancient events of gene conversion, we conducted the analysis with the two paralogous genes for all species simultaneously. Settings for all methods executed in RDP3 were as follows: Sequences were considered to be linear, the p-value cutoff was set to 0.05, and the standard Bonferroni correction was used. In addition, phylogenetic relationships were recovered for each fragment showing signs of gene conversion, in Anthropoids and Catarrhines, following the same procedure described above. ML analyses were performed using the GTR+G model without sequence partitioning. The resulting trees were evaluated with 1000 bootstrap replicates.

Search for open reading frames and promoter analysis

The genomic sequence of the HBB and HBD genes was used to predict the locations and exon-intron structures using the program Genscan [12], followed by sequence alignment of known exon sequences and manual inspection of homology. Protein sequence was obtained using the translation tool implemented in the ExPASy web server (http://www.expasy.org) [3]. Sequence alignment for eutherian HBB proteins and HBD open reading frame were performed using ClustalW [83] implemented in Geneious version 5.5 created by Biomatters (available from http://www.geneious.com/). To identify conserved motifs previously shown to be essential for HBB and HBD expression, promoter sequences located 5’ to the two β-like globin genes (~200bp) were aligned using the same approach as above. To confirm transcription factor binding site matches identified within the alignments, we used MatInspector implemented in the Genomatix Software Suite (http://www.genomatix.de/index.html).

Evolutionary rate estimates and selection tests

Pairwise sequence divergence was deduced from Jukes and Cantor distance calculated with DnaSP v.5.10 [65]. In our estimates we scored each insertion or deletion, regardless of length, as one difference as in [13]. Divergence times between species were obtained with TimeTree [32]. Maximum-likelihood estimates of dN/dS (ω; dS—synonymous substitution rate and dN—non-synonymous substitution rate) were carried out using the codeml program from the software package Phylogenetic Analysis by Maximum Likelihood—PAML version 4.8 [89]. To investigate the selective pressures that have shaped the evolution of HBB and HBD genes we first calculated dN/dS ratios (M0 model) for each gene separately in the entire mammalian phylogeny. Next, to test the hypothesis of variable selective pressures among HBD in primates, we performed nested branch models using either the one-ratio model calculated for the whole anthropoid phylogeny, the two-ratio estimated for Great Apes and other primates and three-ratio inferred for Great Apes, OWM and NWM [7,88]. Although ω values below 1 (ω < 1) are generally considered as an evidence of purifying selection, to reject the hypothesis of neutral evolution all models were compared with a null model where ω was fixed to 1 (ω = 1). The significance was obtained with likelihood ratio tests (LRT) which were calculated as twice the variation of the likelihoods (-2Δl) with a χ² distribution. For the calculation of HBD ω values, pseudogene sequences were only included after the removal of positions affected by premature stop codons and frameshift mutations. In the specific case of lemur species, the HBD sequences were excluded from the analysis given their hybrid ψβ/δ nature [37].

Results

Evolutionary history of the two adult β-like genes in mammals

We conducted a phylogenetic analysis of the adult β-like globin genes, HBD and HBB, in a diverse dataset, including monotremes, marsupials and placental mammals, and one avian species, which was used as outgroup. The phylogenies obtained with ML, BI and the codon model approach were similar when either the complete gene sequence or the CDS were used (Fig 1A–1D and S4 Fig).

Download:

Fig 1. Phylograms depicting relationships among adult β-like genes in mammals.

The phylogeny reconstructions were performed using two methods: A), C) Maximum Likelihood and B), D) Bayesian Inference; trees were estimated using the A), B) coding sequence and C), D) complete gene sequence. Branch support values are given on the internodes. Red branches represent the Great Apes, green the OWM and blue the NWM; the pink and orange branches represent the common branch of Catarrhines and Anthropoids, respectively.

https://doi.org/10.1371/journal.pone.0123365.g001

However, trees obtained using the complete gene were different from those estimated with the CDS only (Fig 1 and S4 Fig). In the phylogenies based on the CDS the two β-like paralogous genes from the same species cluster together, consistent with a process of interparalogue gene conversion, referred to as “concerted evolution” and previously described in several taxa [19,31,40,52,80]. The exceptions to this general pattern occur in Anthropoid primates (monkeys and apes) where HBD and HBB are grouped into two reciprocally monophyletic groups, and in the lemur species. This could lead to the erroneous interpretation that HBD arose recently through a duplication in primate evolution, as initially thought [18], but more likely reflects a gene conversion event in the common ancestor of those primate lineages. In the lemur species, the exception is easily explained by the fact that in the ancestry of lemurs a hybrid ψβ/δ pseudogene was created by unequal crossing-over between misaligned HBD and HBBP1 sequences [37]. Finally, an HBD ortholog is absent from rat, and in mouse we found no evidence for interparalog gene conversion between HBD and HBB orthologs (HBD-T1 and HBB-T1), in agreement with previous results [33]. Phylogenetic analyses based on the entire gene sequence lead to better-supported trees, which more reliably replicate the species tree, recovering the true evolutionary history of gene duplication. This phylogeny is consistent with previous results supporting a duplication event of HBB and HBD genes after the marsupial/eutherian split [53]. Both ML and BI trees contain two major sister clades, one comprising the HBD gene copies of most species and a second clade mainly with the HBB genes. Within the HBB clade, the HBD genes from galago, cow, pig and dolphins are clustered with each of their HBB paralogs. In galago, it is well documented that the ancestral HBD gene was almost completely converted by the HBB gene, explaining the monophyletic pattern observed [80]. According to a recent analysis the phylogenetic incongruence seen in the latter three species is likely due to independent gene conversion events that followed an extensive unequal cross-over spanning HBD intron 2 in the stem lineage of cetartiodactyls [19]. Our results also confirm that the interparalog conversion is largely restricted to the coding regions of globin genes [19,26,30,31,40,45,53], because the incorporation of the phylogenetic information of non-coding sequence, including intronic and 5’ and 3’ flanking sequence, proved to be especially useful in assigning orthologous relationships.

Analysis of gene structure

In all species examined the HBB gene retains an intact ORF with conserved donor/acceptor splice sites encoding a polypeptide of 144–146 amino acids (S1 Fig). By contrast, the HBD gene has accumulated several inactivating mutations throughout the mammalian phylogeny. These include the introduction of premature stop codons, small insertions and deletions and mutations in the consensus donor (GT) or acceptor (AG) splice sites (S2 Fig). Moreover, HBB and HBD also display different across species conservation patterns (S3 Fig). The HBB promoter is highly conserved and contains conserved consensus TATA, CAAT and EKLF (Erythroid Krüppel-like factor) binding motifs in most species examined. Only cat and cow show substitutions in the EKLF consensus sequence (S3A Fig) but present an upstream EKLF binding element, identified by MatInspector, which is likely to replace the possible disrupted EKLF motif. A lower conservation in HBD promoter region is readily apparent from the difficulty in obtaining a good multiple alignment for all species. We detected three species, galago, cat and microbat, in which HBD has a β-like promoter, acquired through independent gene conversion events (S3A Fig), as previously shown [19,80]. In anthropoids, high sequence homology was observed at the HBD promoter region (S3B Fig), with a conserved consensus TATA binding motif, a functional GATA-1 motif [47] close to the mutated CAAC box, and lack of the EKLF binding element in all these species. These features of the HBD promoter region have been shown to be responsible for the low expression of the adult δ-globin gene [63,81,82]. The remaining species (S3C Fig) all share the major defect in the proximal δ-promoter, the absence of a consensus EKLF-binding motif, but do not have the GATA-1 motif common to all other δ-like promoters. The lack of various conserved motifs in the HBD promoters that are crucial for β-like globin expression suggests that they are transcriptionally inefficient in tarsier, mouse, rabbit, dolphin, cow, pig, horse, megabat and armadillo.

Recombination events in primates

In the phylogenetic analyses we did not detect evidence of recombination events between HBD and HBB in anthropoid primates. This result is in agreement with other studies that proposed a δ-β gene conversion in the anthropoid stem [40]. However, it has been suggested that further gene conversions occurred independently in catarrhine and platyrrhine lineages [37,45,60]. Although we did not find phylogenetic evidence for gene conversions within these primate lineages, we cannot exclude the possibility of short-tract gene conversion events not detectable by phylogenetic analysis. Therefore, we sought to re-examine the possibility of lineage-specific gene conversion taking advantage of a more comprehensive sample of primate species and the use of multiple methods implemented in the software package RDP3 [43] for detecting recombination signals and putative recombinant sequences. We found robust signals for four independent interparalog gene conversion events, in which portions of HBB were copied onto HBD (summarized in Table 1).

Download:

Table 1. Summary of gene conversion analysis for primate HBD and HBB paralogues.

https://doi.org/10.1371/journal.pone.0123365.t001

Two of these gene conversion events, 1 and 2, were detected by all seven methods, and corroborate independent gene conversions previously identified in tarsier and galago [40,80]. A third event was detected, corresponding to the conversion of the first HBD exon and intron by HBB sequences, in catarrhines (OWM and Great Apes). Since this event was detected in orthologous HBD copies of all catarrhines represented, it most likely took place before the divergence of OWM and Great Apes, along the pink branch in Fig 1A. The event number 4 corresponds to a more extensive conversion tract present in all anthropoids (platyrrhines and catarrhines), suggesting that those sequences have all descended from an anthropoid ancestor sequence in which the recombination event occurred (orange branch in Fig 1A). This δ-β conversion appears to extend from the 5’ promoter region till the end of the second exon, however in catarrhines the signal for this older gene conversion is restricted to the second exon, because a subsequent event (identified as number 3) overprinted part of the older one.

In order to assess whether the conversion events 3 and 4 have been correctly identified, we constructed and compared two phylogenetic trees, one with the region containing evidence for an older gene conversion in both platyrrhines and catarrhines (nucleotide 367–632), and a second one with the portion of the alignment between the inferred breakpoints in event 3 (nucleotide 97–385) where a second, more recent conversion event occurred in the catarrhine stem (S5 Fig). The topology of the first tree (S5A Fig) is consistent with a gene conversion in the common ancestor of Anthropoids, in agreement with results obtained with other methods [31,40,45,73]. In contrast, in the second tree (S5B Fig) the HBD and HBB genes group together within the catarrhine and platyrrhine lineages, as expected under the hypothesis of a conversion event restricted to catarrhines. This tree topology could also indicate that parallel gene conversion events have occurred in the stem of catarrhines and platyrrhines; however, no gene conversion event in platyrrhines has been detected with any of the methods used in our analysis. Nevertheless, this alternative hypothesis remains disputable given that an older conversion event occurring between still very closely related sequences would be very difficult to detect. Therefore, although our results do not fully support previous evidence that HBD has been involved in a conversion in platyrrhines [60], it cannot be ruled out. It is noteworthy that albeit two independent gene conversion events occurred in different Anthropoid lineages, the GATA-1 motif in the HBD promoter region remained intact.

Evolutionary rates and functional constraints in primates

From our previous analysis, it is apparent that while HBD has diverged at markedly different rates in different primate lineages, as shown by the variable branch lengths in the phylogenetic trees (Fig 1), anthropoid HBDs share a high sequence identity not only in their coding region but also at the promoter. As a first measure of the rate of HBD evolution, we compared the genetic distance between humans and 13 other primate species, for both adult β-like globin genes. HBBP1 was also included in the analysis due to the unusually slow substitution rates previously reported for this pseudogene estimated by comparison of human, gorilla and chimpanzee sequences [13]. Genetic distances were then plotted against the corresponding divergence times for each pairwise comparison, and the linear regression trend line was estimated for each group, as shown in Fig 2.

Download:

Fig 2. Genetic distance vs divergence times between human and different primate species (Anthropoids and Prosimians) for β-like genes.

Circles and dotted lines correspond to introns while diamonds and solid lines correspond to exons. Divergence times between humans and other species were obtained with TimeTree [32] and are as follows: Ptr: 6,3 Myr; Ggo 8,8 Myr; Ppy: 15.7 Myr; Nle 20,4 Myr; OWM (Mcc, Panu, Cgue and Caa): 29 Myr; NWM (Sbol, Cjac and Anan): 42,6 Myr; Tsyr: 65,5 Myr; Ogar: 74 Myr. Linear regression trend lines were set to intercept the origin.

https://doi.org/10.1371/journal.pone.0123365.g002

From the slope of the trend lines and the r² values we are able to compare the rate of intron and exon evolution and its constancy over time, respectively. The results presented in Fig 2 and S2 Table show that, overall, exons evolved at a lower rate than introns, except for HBBP1, in which nucleotide differences are more homogeneously distributed between exons and introns, as expected for a pseudogene. In Prosimians, there is a trend towards increasing evolutionary rates, more pronounced in non-coding regions (introns and HBPP1 exons). Remarkably, the rate of evolution of HBD exons has remained relatively constant across primate evolution (r² = 0.95) and is comparable to that of HBB exons, even though in higher primates HBD is either silent or contributes to only a very small fraction of adult Hb. To gain further insight into the possible functional constraints that have shaped the evolutionary history of HBD in primates, we calculated dN/dS (ω) ratios under alternative models of gene evolution (table 2).

Download:

Table 2. Parameter Estimates and Likelihood Scores under Different Branch Models.

https://doi.org/10.1371/journal.pone.0123365.t002

First, we estimated single ω values for the entire mammalian phylogeny (M0 model), for either HBB or HBD genes. The observed ω values were significantly lower than 1 (ω_HBB = 0.26 and ω_HBD = 0.40) pointing to an overall conservation of HBB and HBD. The ω obtained for HBB is the expected for a functional gene and is in agreement with previous estimates [1]. In the case of HBD orthologs, the detected signal of purifying selection may be an outcome of gene conversion by the HBB gene in multiple lineages, mostly in the coding region. As illegitimate recombination between HBD and HBB has been noticeably reduced since the last common ancestor of Anthropoids 65.5 Mya, we examined the extent of the selective pressures exerted in this specific clade. Overall, the significance of the low value obtained under the one-ratio model (ω_{HBD_Anthropoids} = 0.31) rejects the hypothesis of a neutral evolution of HBD in anthropoids. Then, to examine whether HBD has been subject to variable selective constraints among different anthropoid clades (Fig 1A), we applied the two-ratio model separating the phylogenetic group of Great Apes from all remaining primates (OWM and NWM). In addition, we also applied the three-ratio model, in which ω was allowed to vary across Great Apes, OWM and NWM clades. The comparison of the one-ratio and the two-ratio models showed that the two-ratio presents an improved fit to anthropoid HBD evolution and that it did not differ significantly from the three-ratio model. Nevertheless, in both two-ratio and three-ratio models the Great Apes clade showed an extremely low ω value (ω_{HBD_Apes} = 0.06), while the OWM+NWM branches present a higher ω (ω_{HBD_OWM+NWM} = 0.43), but still suggesting a constrained evolution. Although our results indicate that HBD might have experienced different selective pressures throughout primate evolution, these estimates corroborated a high conservation of HBD in Anthropoid lineages that is unlikely related to protein function, since in most primate species this gene is either weakly expressed or not transcribed at all.

Discussion

In humans and in chimpanzees, unusually high levels of HBD sequence conservation, when compared to functional paralogs, have been described [49]. Such pattern of conservation has been difficult to reconcile with the negligible expression of HbA₂. Moreover, the evolutionary history of HBD is complex and orthologous relationships among HBD and its paralog gene (HBB) have been obscured by a history of recurrent gene conversion and unequal crossing overs, throughout eutherian evolution [19,52–54]. Here we gained insight into the evolutionary history of HBD and its likely regulatory role in the fetal-to-adult switch unique of Anthropoids, by performing a comprehensive phylogenetic and comparative analysis of the two adult β-like globin genes in a wide range of mammalian taxa. The results from our phylogenetic reconstruction are in agreement with previous findings which demonstrated that HBD duplication occurred before the radiation of Eutheria [53]. The obtained tree topology is also consistent with a history of concerted evolution between HBD and HBB that has created chimeric β/δ fusion genes in multiple, independent lineages [19,31,40,54,80]. However, our results show that primates represent an exception to this common trend, given that phylogenetic relationships are maintained throughout this lineage, suggesting that illegitimate recombination between HBD and HBB has been noticeably reduced since the last common ancestor of Anthropoids 65,5 Mya. Indeed, the recombination analyses here presented demonstrate that even though a certain level of gene conversion has occurred in some Anthropoid lineages, these events have taken place in the stem of the major branches (platyrrhines and catarrhines), as previously inferred from smaller datasets [31,40,45,60]. Moreover, gene conversion in Anthropoids was restricted to shorter regions (exon and intron 1) than those most frequently identified in other species (exons 1, 2, and 3 and intron 1) [19,30,31,40,80]. The distinct phylogenetic patterns between Anthropoid primates and nearly all other species, obtained when using the coding sequence, reflect the relative extent of these events. Although most mammals retain both adult β-like gene copies, only the HBB gene is functional and essential. HBD has apparently become dispensable and in several lineages has pseudogenized. In fact, we confirmed that HBD inactivation occurred in a wide range of eutherian species by the accumulation of loss-of-function mutations that disrupted either the ORF or the promoter region. It is noteworthy that when considering extant mammalian species it is only possible to establish bona-fide homology relationships between HBD orthologs among Anthropoid primates. Moreover, evolutionary tests suggest that HBD is evolving under selective constraints across different Anthropoid species, as hypothesized several decades ago [13,24,37,74]. Such high level of conservation is at odds with the variable expressivity of the δ-globin chain, which is absent in OWM [45,46] and ranges from 1% concentration in hominoids [10] to 6% in NWM [74]. Selective constraints on protein evolution do not seem a plausible explanation for the signal of purifying selection detected, since to date, HBA₂ has no recognized physiological function [62,76,77].

Interestingly, a role for HBD and HBBP1 in the regulatory mechanisms coordinating the fetal-to adult switch has been proposed in early independent studies [5,13,24,26,55]. Taking into account that the mechanism of Hb switch is common to all simian primates [38], we might expect to find similar patterns of conservation and diversity in ortholog HBD sequences for a 65.5 Myr time frame. Accordingly, the patterns of conservation we have now uncovered perfectly overlap with HBG duplication and the acquisition of a fetally expressed hemoglobin in anthropoid primates. Noteworthy, we detected in all anthropoid primates a conserved functional GATA-1 motif in the promoter of HBD, which has remained intact despite recurrent gene conversion events overlapping the promoter region among these lineages. Considering that HBD has very low expression levels in anthropoids, the conservation of a functional GATA-1 binding motif suggests other functional constraints rather than positive regulation of δ-globin gene expression. Indeed Gaudry, et al. [19] demonstrated that only late expressed β-like globin genes retaining an HBB-like promoter are efficiently transcribed. Developmental regulation of gene expression at the β-globin cluster involves the formation of chromatin loops mediated by several transcription factors and cofactors [66,68]. It has been shown that GATA-1, along with other cofactors, is required for efficient long-range chromatin interactions between LCR and β-like globin genes, namely at the time of γ- to β- globin switch [9,39,85]. Importantly, is has also been demonstrated that the HBD upstream region harbors a binding site for BCL11A, which is a biochemically validated and fundamental switching factor necessary for fetal hemoglobin silencing [67,69]. Remarkably, strong interactions between the LCR and the region encompassing both HBD and HBBP1 were uncovered by chromosome conformation (3C and 5C) analyses at the β-globin locus [6,17,70]. Collectively, these findings suggest that HBD and HBBP1 might be involved in chromatin looping in the human-globin cluster, a crucial mechanism for temporal coordination of gene expression [16,36]. Interestingly, the observed differences in the rate of evolution between the branches leading to Great Apes and the common branch of OWM and NWM suggest different selective pressures, which may reflect alternative mechanisms of controlling expression in the β-globin cluster among Anthropoids. Selective constraints on the protein function cannot be completely ruled out, although evidence of functional relevance of HbA₂ is lacking. HbA₂ has features that are nearly identical to those of HbA [15] but, even though in the absence of β-chain production in β-thalassemia major it becomes the predominant oxygen carrier, it never effectively replaces HbA function. Concordantly, in humans mutations in HBD are per se clinically silent [71,76]. The evolutionary history of HBD in mammals has been shaped by concerted evolution, however our results show that in the extant species of Anthropoids gene conversion events have not been frequent, and when they do occur the exchanged sequence tract is short. The low frequency of gene conversion as well as the conservation of a motif involved in chromatin remodeling could be an outcome of stronger selective constraints acting on HBD in anthropoids, as previously reported for other functionally important regions [50,58,92]. These results are also consistent with previous evidence of purifying selection reducing HBD genetic diversity in human populations and in chimpanzees [49]. We have now characterized the evolution of HBD in mammals and found that the unusual high levels of conservation in this genomic regions are shared across several primate species. In the light of recent advances in the understanding of the β-globin cluster regulation, we propose that the similar evolutionary trajectory of HBD in Anthropoids is due to functional constraints related with the intricate process of chromatin and protein interactions coordinating the developmental expression of β-like globin genes.

Supporting Information

S1 Fig. Sequence alignment for eutherian HBB proteins.

https://doi.org/10.1371/journal.pone.0123365.s001

(PDF)

S2 Fig. Sequence alignment of eutherian HBD open reading frame.

Blue and purple filled boxes mark the exons and donnor/acceptor splice sites, respectively. Dots represent nucleotide identities to the human sequence that was set as reference. Coloured nucleotides indicate changes to the human sequence and aminoacid alterations are marked by filled coloured boxes. The lemur species were excluded from the analysis given that their hybrid ψβ/δ pseudogene [37] generates multiple misalignments.

https://doi.org/10.1371/journal.pone.0123365.s002

(PDF)

S3 Fig. Sequence alignment of eutherian HBB and HBD promoters.

A) HBB and HBB-like HBD promoters; B) Anthropoid HBD promoters and C) HBD promoters lacking the TF binding motifs which are conserved in HBB-like and HBD-like promoters. Conserved binding motifs are indicated in grey boxes. Again, the lemur species were excluded from the analysis.

https://doi.org/10.1371/journal.pone.0123365.s003

(PDF)

S4 Fig. Phylograms depicting relationships among adult β-like genes in mammals.

The phylogenic tree, based on the coding sequence, was constructed using the Goldman–Yang codon model. Branch support values, obtained using the approximate Likelihood Ration Test (aLRT), are given on the internodes.

https://doi.org/10.1371/journal.pone.0123365.s004

(PDF)

S5 Fig. Maximum Likelihood phylograms depicting relationships among β-like genes of Anthropoids.

The phylogeny reconstructions were based on A) the portion of the alignment that contain evidence of the anthropoid gene conversion (nucleotide 367–632) and B) the portion of the alignment between the inferred breakpoint in event 3 (nucleotide 97–385). Bootstrap branch support (1000 replicates) are given on the internodes.

https://doi.org/10.1371/journal.pone.0123365.s005

(PDF)

S1 Table. Listing of HBB and HBD sequences used for phylogeny reconstructions.

https://doi.org/10.1371/journal.pone.0123365.s006

(PDF)

S2 Table. Evolutionary rates based on Jukes Cantor distance.

https://doi.org/10.1371/journal.pone.0123365.s007

(PDF)

Author Contributions

Conceived and designed the experiments: AM AA. Performed the experiments: AM. Analyzed the data: AM AML SS MJP AA RMH. Contributed reagents/materials/analysis tools: AM AML SS RMH. Wrote the paper: AM AML SS MJP AA RMH.

References

1. Aguileta G, Bielawski JP, Yang Z. Gene conversion and functional divergence in the beta-globin gene family. J Mol Evol. 2004;59: 177–189. pmid:15486692
- View Article
- PubMed/NCBI
- Google Scholar
2. Armougom F, Moretti S, Poirot O, Audic S, Dumas P, Schaeli B, et al. Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee. Nucleic Acids Res. 2006;34: W604–608. pmid:16845081
- View Article
- PubMed/NCBI
- Google Scholar
3. Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Research. 2012;40: W597–W603. pmid:22661580
- View Article
- PubMed/NCBI
- Google Scholar
4. Bank A. Regulation of human fetal hemoglobin: new players, new complexities. Blood. 2006;107: 435–443. pmid:16109777
- View Article
- PubMed/NCBI
- Google Scholar
5. Bank A, Mears JG, Ramirez F. Disorders of human hemoglobin. Science. 1980;207: 486–493. pmid:7352255
- View Article
- PubMed/NCBI
- Google Scholar
6. Beauchemin H, Trudel M. Evidence for a bigenic chromatin subdomain in regulation of the fetal-to-adult hemoglobin switch. Mol Cell Biol. 2009;29: 1635–1648. pmid:19114559
- View Article
- PubMed/NCBI
- Google Scholar
7. Bielawski JP, Yang Z. Maximum likelihood methods for detecting adaptive evolution after gene duplication. J Struct Funct Genomics. 2003;3: 201–212. pmid:12836699
- View Article
- PubMed/NCBI
- Google Scholar
8. Boni MF, Posada D, Feldman MW. An Exact Nonparametric Method for Inferring Mosaic Structure in Sequence Triplets. Genetics. 2007;176: 1035–1047. pmid:17409078
- View Article
- PubMed/NCBI
- Google Scholar
9. Bottardi S, Ross J, Bourgoin V, Fotouhi-Ardakani N, Affar el B, Trudel M, et al. Ikaros and GATA-1 combinatorial effect is required for silencing of human gamma-globin genes. Mol Cell Biol. 2009;29: 1526–1537. pmid:19114560
- View Article
- PubMed/NCBI
- Google Scholar
10. Boyer S, Crosby E, Noyes A, Fuller G, Leslie S, Donaldson L, et al. Primate hemoglobins: Some sequences and some proposals concerning the character of evolution and mutation. Biochemical Genetics. 1971;5: 405–448. pmid:4999925
- View Article
- PubMed/NCBI
- Google Scholar
11. Bulger M, Groudine M. Looping versus linking: toward a model for long-distance gene activation. Genes Dev. 1999;13: 2465–2477. pmid:10521391
- View Article
- PubMed/NCBI
- Google Scholar
12. Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997;268: 78–94. pmid:9149143
- View Article
- PubMed/NCBI
- Google Scholar
13. Chang LY, Slightom JL. Isolation and nucleotide sequence analysis of the beta-type globin pseudogene from human, gorilla and chimpanzee. J Mol Biol. 1984;180: 767–784. pmid:6098690
- View Article
- PubMed/NCBI
- Google Scholar
14. Czelusniak J, Goodman M, Hewett-Emmett D, Weiss ML, Venta PJ, Tashian RE. Phylogenetic origins and adaptive evolution of avian and mammalian haemoglobin genes. Nature. 1982;298: 297–300. pmid:6178039
- View Article
- PubMed/NCBI
- Google Scholar
15. de Bruin SH, Janssen LHM. Comparison of the oxygen and proton binding behavior of human hemoglobin A and A2. Biochimica et Biophysica Acta (BBA)—Protein Structure. 1973;295: 490–494.
- View Article
- Google Scholar
16. Dean A. On a chromosome far, far away: LCRs and gene expression. Trends Genet. 2006;22: 38–45. pmid:16309780
- View Article
- PubMed/NCBI
- Google Scholar
17. Dostie J, Richmond TA, Arnaout RA, Selzer RR, Lee WL, Honan TA, et al. Chromosome Conformation Capture Carbon Copy (5C): a massively parallel solution for mapping interactions between genomic elements. Genome Res. 2006;16: 1299–1309. pmid:16954542
- View Article
- PubMed/NCBI
- Google Scholar
18. Efstratiadis A, Posakony JW, Maniatis T, Lawn RM, O'Connell C, Spritz RA, et al. The structure and evolution of the human beta-globin gene family. Cell. 1980;21: 653–668. pmid:6985477
- View Article
- PubMed/NCBI
- Google Scholar
19. Gaudry MJ, Storz JF, Butts GT, Campbell KL, Hoffmann FG. Repeated Evolution of Chimeric Fusion Genes in the β-Globin Gene Family of Laurasiatherian Mammals. Genome Biology and Evolution. 2014.
- View Article
- Google Scholar
20. Gibbs MJ, Armstrong JS, Gibbs AJ. Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000;16: 573–582. pmid:11038328
- View Article
- PubMed/NCBI
- Google Scholar
21. Gil M, Zanetti MS, Zoller S, Anisimova M. CodonPhyML: Fast Maximum Likelihood Phylogeny Estimation under Codon Substitution Models. Molecular Biology and Evolution. 2013.
- View Article
- Google Scholar
22. Goldman N, Yang Z. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994;11: 725–736. pmid:7968486
- View Article
- PubMed/NCBI
- Google Scholar
23. Goodman M, Czelusniak J, Koop BF, Tagle DA, Slightom JL. Globins: a case study in molecular phylogeny. Cold Spring Harb Symp Quant Biol. 1987;52: 875–890. pmid:3454296
- View Article
- PubMed/NCBI
- Google Scholar
24. Goodman M, Koop BF, Czelusniak J, Weiss ML. The eta-globin gene. Its long evolutionary history in the beta-globin gene family of mammals. J Mol Biol. 1984;180: 803–823. pmid:6527390
- View Article
- PubMed/NCBI
- Google Scholar
25. Goodman M, Moore GW, Matsuda G. Darwinian evolution in the genealogy of haemoglobin. Nature. 1975;253: 603–608. pmid:1089897
- View Article
- PubMed/NCBI
- Google Scholar
26. Hardies SC, Edgell MH, Hutchison CA 3rd. Evolution of the mammalian beta-globin gene cluster. J Biol Chem. 1984;259: 3748–3756. pmid:6706976
- View Article
- PubMed/NCBI
- Google Scholar
27. Hardison R. Evolution of hemoglobin and its genes. Cold Spring Harb Perspect Med. 2012;2: a011627. pmid:23209182
- View Article
- PubMed/NCBI
- Google Scholar
28. Hardison R. Hemoglobins from bacteria to man: evolution of different patterns of gene expression. J Exp Biol. 1998;201: 1099–1117. pmid:9510523
- View Article
- PubMed/NCBI
- Google Scholar
29. Hardison R. Organization, evolution and regulation of the globin genes. In: Steinberg MH, Forget BG, Higgs DR, Nagel RL, editors. Disorders of Hemoglobin: Genetics, Pathophysiology, and Clinical Management. Cambridge: Cambridge University Press; 2001. pmid:12779271
30. Hardison RC. 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. Mol Biol Evol. 1984;1: 390–410. pmid:6599973
- View Article
- PubMed/NCBI
- Google Scholar
31. Hardison RC, Margot JB. Rabbit globin pseudogene psi beta 2 is a hybrid of delta- and beta-globin gene sequences. Mol Biol Evol. 1984;1: 302–316. pmid:6599969
- View Article
- PubMed/NCBI
- Google Scholar
32. Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006;22: 2971–2972. pmid:17021158
- View Article
- PubMed/NCBI
- Google Scholar
33. Hoffmann FG, Opazo JC, Storz JF. New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents. Mol Biol Evol. 2008;25: 2589–2600. pmid:18780876
- View Article
- PubMed/NCBI
- Google Scholar
34. Hoffmann FG, Storz JF. The alphaD-globin gene originated via duplication of an embryonic alpha-like globin gene in the ancestor of tetrapod vertebrates. Mol Biol Evol. 2007;24: 1982–1990. pmid:17586601
- View Article
- PubMed/NCBI
- Google Scholar
35. Hoffmann FG, Storz JF, Gorr TA, Opazo JC. Lineage-specific patterns of functional diversification in the alpha- and beta-globin gene families of tetrapod vertebrates. Mol Biol Evol. 2010;27: 1126–1138. pmid:20047955
- View Article
- PubMed/NCBI
- Google Scholar
36. Holwerda S, de Laat W. Chromatin loops, gene positioning, and gene expression. Front Genet. 2012;3: 217. pmid:23087710
- View Article
- PubMed/NCBI
- Google Scholar
37. Jeffreys AJ, Barrie PA, Harris S, Fawcett DH, Nugent ZJ, Boyd AC. Isolation and sequence analysis of a hybrid delta-globin pseudogene from the brown lemur. J Mol Biol. 1982;156: 487–503. pmid:6214636
- View Article
- PubMed/NCBI
- Google Scholar
38. Johnson RM, Gumucio D, Goodman M. Globin gene switching in primates. Comp Biochem Physiol A Mol Integr Physiol. 2002;133: 877–883. pmid:12443943
- View Article
- PubMed/NCBI
- Google Scholar
39. Keys JR, Tallack MR, Zhan Y, Papathanasiou P, Goodnow CC, Gaensler KM, et al. A mechanism for Ikaros regulation of human globin gene switching. Br J Haematol. 2008;141: 398–406. pmid:18318763
- View Article
- PubMed/NCBI
- Google Scholar
40. Koop BF, Siemieniak D, Slightom JL, Goodman M, Dunbar J, Wright PC, et al. Tarsius delta- and beta-globin genes: conversions, evolution, and systematic implications. J Biol Chem. 1989;264: 68–79. pmid:2491855
- View Article
- PubMed/NCBI
- Google Scholar
41. Lanfear R, Calcott B, Ho SY, Guindon S. Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29: 1695–1701. pmid:22319168
- View Article
- PubMed/NCBI
- Google Scholar
42. Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000;16: 562–563. pmid:10980155
- View Article
- PubMed/NCBI
- Google Scholar
43. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26: 2462–2463. pmid:20798170
- View Article
- PubMed/NCBI
- Google Scholar
44. Martin DP, Posada D, Crandall KA, Williamson C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses. 2005;21: 98–102. pmid:15665649
- View Article
- PubMed/NCBI
- Google Scholar
45. Martin SL, Vincent KA, Wilson AC. Rise and fall of the delta globin gene. J Mol Biol. 1983;164: 513–528. pmid:6188843
- View Article
- PubMed/NCBI
- Google Scholar
46. Martin SL, Zimmer EA, Kan YW, Wilson AC. Silent delta-globin gene in Old World monkeys. Proc Natl Acad Sci U S A. 1980;77: 3563–3566. pmid:6251467
- View Article
- PubMed/NCBI
- Google Scholar
47. Matsuda M, Sakamoto N, Fukumaki Y. Delta-thalassemia caused by disruption of the site for an erythroid-specific transcription factor, GATA-1, in the delta-globin gene promoter. Blood. 1992;80: 1347–1351. pmid:1515647
- View Article
- PubMed/NCBI
- Google Scholar
48. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees; 2010. pp. 1–8.
- View Article
- Google Scholar
49. Moleirinho A, Seixas S, Lopes AM, Bento C, Prata MJ, Amorim A. Evolutionary constraints in the beta-globin cluster: the signature of purifying selection at the delta-globin (HBD) locus and its role in developmental gene regulation. Genome Biol Evol. 2013;5: 559–571. pmid:23431002
- View Article
- PubMed/NCBI
- Google Scholar
50. Noonan JP, Grimwood J, Schmutz J, Dickson M, Myers RM. Gene conversion and the evolution of protocadherin gene cluster diversity. Genome Res. 2004;14: 354–366. pmid:14993203
- View Article
- PubMed/NCBI
- Google Scholar
51. Notredame C, Higgins DG, Heringa J. T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302: 205–217. pmid:10964570
- View Article
- PubMed/NCBI
- Google Scholar
52. Opazo JC, Hoffmann FG, Storz JF. Differential loss of embryonic globin genes during the radiation of placental mammals. Proceedings of the National Academy of Sciences. 2008;105: 12950–12955. pmid:18755893
- View Article
- PubMed/NCBI
- Google Scholar
53. Opazo JC, Hoffmann FG, Storz JF. Genomic evidence for independent origins of β-like globin genes in monotremes and therian mammals. Proceedings of the National Academy of Sciences. 2008;105: 1590–1595. pmid:18216242
- View Article
- PubMed/NCBI
- Google Scholar
54. Opazo JC, Sloan AM, Campbell KL, Storz JF. Origin and Ascendancy of a Chimeric Fusion Gene: The β/δ-Globin Gene of Paenungulate Mammals. Molecular Biology and Evolution. 2009;26: 1469–1478. pmid:19332641
- View Article
- PubMed/NCBI
- Google Scholar
55. Ottolenghi S, Giglioni B, Comi P, Gianni AM, Polli E, Acquaye CT, et al. Globin gene deletion in HPFH, delta (o) beta (o) thalassaemia and Hb Lepore disease. Nature. 1979;278: 654–657. pmid:450068
- View Article
- PubMed/NCBI
- Google Scholar
56. Padidam M, Sawyer S, Fauquet CM. Possible emergence of new geminiviruses by frequent recombination. Virology. 1999;265: 218–225. pmid:10600594
- View Article
- PubMed/NCBI
- Google Scholar
57. Patel VS, Cooper SJ, Deakin JE, Fulton B, Graves T, Warren WC, et al. Platypus globin genes and flanking loci suggest a new insertional model for beta-globin evolution in birds and mammals. BMC Biol. 2008;6: 34. pmid:18657265
- View Article
- PubMed/NCBI
- Google Scholar
58. Petronella N, Drouin G. Purifying selection against gene conversions in the folate receptor genes of primates. Genomics. 2014;103: 40–47. pmid:24184359
- View Article
- PubMed/NCBI
- Google Scholar
59. Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci U S A. 2001;98: 13757–13762. pmid:11717435
- View Article
- PubMed/NCBI
- Google Scholar
60. Prychitko T, Johnson RM, Wildman DE, Gumucio D, Goodman M. The phylogenetic history of New World monkey beta globin reveals a platyrrhine beta to delta gene conversion in the atelid ancestry. Mol Phylogenet Evol. 2005;35: 225–234. pmid:15737593
- View Article
- PubMed/NCBI
- Google Scholar
61. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6; 2014. Available: http://beast.bio.ed.ac.uk/Tracer.
62. Ranney HM, Lam R, Rosenberg G. Some properties of hemoglobin A2. American Journal of Hematology. 1993;42: 107–111. pmid:8416283
- View Article
- PubMed/NCBI
- Google Scholar
63. Ristaldi MS, Casula S, Porcu S, Marongiu MF, Pirastu M, Cao A. Activation of the delta-globin gene by the beta-globin gene CACCC motif. Blood Cells Mol Dis. 1999;25: 193–209. pmid:10575545
- View Article
- PubMed/NCBI
- Google Scholar
64. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61: 539–542. pmid:22357727
- View Article
- PubMed/NCBI
- Google Scholar
65. Rozas J. DNA sequence polymorphism analysis using DnaSP. Methods Mol Biol. 2009;537: 337–350. pmid:19378153
- View Article
- PubMed/NCBI
- Google Scholar
66. Sankaran VG, Orkin SH. The switch from fetal to adult hemoglobin. Cold Spring Harb Perspect Med. 2013;3: a011643. pmid:23209159
- View Article
- PubMed/NCBI
- Google Scholar
67. Sankaran VG, Xu J, Byron R, Greisman HA, Fisher C, Weatherall DJ, et al. A Functional Element Necessary for Fetal Hemoglobin Silencing. New England Journal of Medicine. 2011;365: 807–814. pmid:21879898
- View Article
- PubMed/NCBI
- Google Scholar
68. Sankaran VG, Xu J, Orkin SH. Advances in the understanding of haemoglobin switching. Br J Haematol. 2010;149: 181–194. pmid:20201948
- View Article
- PubMed/NCBI
- Google Scholar
69. Sankaran VG, Xu J, Ragoczy T, Ippolito GC, Walkley CR, Maika SD, et al. Developmental and species-divergent globin switching are driven by BCL11A. Nature. 2009;460: 1093–1097. pmid:19657335
- View Article
- PubMed/NCBI
- Google Scholar
70. Sanyal A, Lajoie BR, Jain G, Dekker J. The long-range interaction landscape of gene promoters. Nature. 2012;489: 109–113. pmid:22955621
- View Article
- PubMed/NCBI
- Google Scholar
71. Schechter AN. Hemoglobin research and the origins of molecular medicine. 3927–3938 p; 2008
- View Article
- Google Scholar
72. Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34: 126–129. pmid:1556748
- View Article
- PubMed/NCBI
- Google Scholar
73. Song G, Hsu CH, Riemer C, Zhang Y, Kim HL, Hoffmann F, et al. Conversion events in gene clusters. BMC Evol Biol. 2011;11: 226. pmid:21798034
- View Article
- PubMed/NCBI
- Google Scholar
74. Spritz RA, Giebel LB. The structure and evolution of the spider monkey delta-globin gene. Mol Biol Evol. 1988;5: 21–29. pmid:2833675
- View Article
- PubMed/NCBI
- Google Scholar
75. Stamatakis A. RAxML Version 8: A tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics. 2014.
- View Article
- Google Scholar
76. Steinberg MH, Adams JG 3rd. Hemoglobin A2: origin, evolution, and aftermath. Blood. 1991;78: 2165–2177. pmid:1932737
- View Article
- PubMed/NCBI
- Google Scholar
77. Steinberg MH, et al. Disorders of Hemoglobin. Cambridge University Press; 2009.
78. Storz JF, Opazo JC, Hoffmann FG. Gene duplication, genome duplication, and the functional diversification of vertebrate globins. Mol Phylogenet Evol. 2013;66: 469–478. pmid:22846683
- View Article
- PubMed/NCBI
- Google Scholar
79. Storz JF, Opazo JC, Hoffmann FG. Phylogenetic diversification of the globin gene superfamily in chordates. IUBMB Life. 2011;63: 313–322. pmid:21557448
- View Article
- PubMed/NCBI
- Google Scholar
80. Tagle DA, Slightom JL, Jones RT, Goodman M. Concerted evolution led to high expression of a prosimian primate delta globin gene locus. J Biol Chem. 1991;266: 7469–7480. pmid:2019578
- View Article
- PubMed/NCBI
- Google Scholar
81. Tang DC, Ebb D, Hardison RC, Rodgers GP. Restoration of the CCAAT box or insertion of the CACCC motif activates [corrected] delta-globin gene expression. Blood. 1997;90: 421–427. pmid:9207479
- View Article
- PubMed/NCBI
- Google Scholar
82. Tang DC, Rodgers GP. Activation of the human delta-globin gene promoter in primary adult erythroid cells. Br J Haematol. 1998;103: 835–838. pmid:9858241
- View Article
- PubMed/NCBI
- Google Scholar
83. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22: 4673–4680. pmid:7984417
- View Article
- PubMed/NCBI
- Google Scholar
84. Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W. Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell. 2002;10: 1453–1465. pmid:12504019
- View Article
- PubMed/NCBI
- Google Scholar
85. Vakoc CR, Letting DL, Gheldof N, Sawado T, Bender MA, Groudine M, et al. Proximity among distant regulatory elements at the beta-globin locus requires GATA-1 and FOG-1. Mol Cell. 2005;17: 453–462. pmid:15694345
- View Article
- PubMed/NCBI
- Google Scholar
86. Vincent KA, Wilson AC. Evolution and transcription of old world monkey globin genes. J Mol Biol. 1989;207: 465–480. pmid:2760921
- View Article
- PubMed/NCBI
- Google Scholar
87. Webster MT, Clegg JB, Harding RM. Common 5' beta-globin RFLP haplotypes harbour a surprising level of ancestral sequence mosaicism. Hum Genet. 2003;113: 123–139. pmid:12736816
- View Article
- PubMed/NCBI
- Google Scholar
88. Yang Z. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998;15: 568–573. pmid:9580986
- View Article
- PubMed/NCBI
- Google Scholar
89. Yang Z. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Molecular Biology and Evolution. 2007;24: 1586–1591. pmid:17483113
- View Article
- PubMed/NCBI
- Google Scholar
90. Yang Z, Nielsen R, Goldman N, Pedersen AM. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000;155: 431–449. pmid:10790415
- View Article
- PubMed/NCBI
- Google Scholar
91. Zhang J. Evolution by gene duplication: an update. Trends in Ecology & Evolution. 2003;18: 292–298.
- View Article
- Google Scholar
92. Zhao Z, Hewett-Emmett D, Li WH. Frequent gene conversion between human red and green opsin genes. J Mol Evol. 1998;46: 494–496. pmid:9541545
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Aguileta G, Bielawski JP, Yang Z. Gene conversion and functional divergence in the beta-globin gene family. J Mol Evol. 2004;59: 177–189. pmid:15486692
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Armougom F, Moretti S, Poirot O, Audic S, Dumas P, Schaeli B, et al. Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee. Nucleic Acids Res. 2006;34: W604–608. pmid:16845081
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Research. 2012;40: W597–W603. pmid:22661580
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Bank A. Regulation of human fetal hemoglobin: new players, new complexities. Blood. 2006;107: 435–443. pmid:16109777
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Bank A, Mears JG, Ramirez F. Disorders of human hemoglobin. Science. 1980;207: 486–493. pmid:7352255
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Beauchemin H, Trudel M. Evidence for a bigenic chromatin subdomain in regulation of the fetal-to-adult hemoglobin switch. Mol Cell Biol. 2009;29: 1635–1648. pmid:19114559
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Bielawski JP, Yang Z. Maximum likelihood methods for detecting adaptive evolution after gene duplication. J Struct Funct Genomics. 2003;3: 201–212. pmid:12836699
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Boni MF, Posada D, Feldman MW. An Exact Nonparametric Method for Inferring Mosaic Structure in Sequence Triplets. Genetics. 2007;176: 1035–1047. pmid:17409078
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Bottardi S, Ross J, Bourgoin V, Fotouhi-Ardakani N, Affar el B, Trudel M, et al. Ikaros and GATA-1 combinatorial effect is required for silencing of human gamma-globin genes. Mol Cell Biol. 2009;29: 1526–1537. pmid:19114560
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Boyer S, Crosby E, Noyes A, Fuller G, Leslie S, Donaldson L, et al. Primate hemoglobins: Some sequences and some proposals concerning the character of evolution and mutation. Biochemical Genetics. 1971;5: 405–448. pmid:4999925
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Bulger M, Groudine M. Looping versus linking: toward a model for long-distance gene activation. Genes Dev. 1999;13: 2465–2477. pmid:10521391
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Burge C, Karlin S. Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997;268: 78–94. pmid:9149143
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Chang LY, Slightom JL. Isolation and nucleotide sequence analysis of the beta-type globin pseudogene from human, gorilla and chimpanzee. J Mol Biol. 1984;180: 767–784. pmid:6098690
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Czelusniak J, Goodman M, Hewett-Emmett D, Weiss ML, Venta PJ, Tashian RE. Phylogenetic origins and adaptive evolution of avian and mammalian haemoglobin genes. Nature. 1982;298: 297–300. pmid:6178039
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. de Bruin SH, Janssen LHM. Comparison of the oxygen and proton binding behavior of human hemoglobin A and A2. Biochimica et Biophysica Acta (BBA)—Protein Structure. 1973;295: 490–494.
View Article
Google Scholar

[58] View Article

[59] Google Scholar

[ref16] 16. Dean A. On a chromosome far, far away: LCRs and gene expression. Trends Genet. 2006;22: 38–45. pmid:16309780
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref17] 17. Dostie J, Richmond TA, Arnaout RA, Selzer RR, Lee WL, Honan TA, et al. Chromosome Conformation Capture Carbon Copy (5C): a massively parallel solution for mapping interactions between genomic elements. Genome Res. 2006;16: 1299–1309. pmid:16954542
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref18] 18. Efstratiadis A, Posakony JW, Maniatis T, Lawn RM, O'Connell C, Spritz RA, et al. The structure and evolution of the human beta-globin gene family. Cell. 1980;21: 653–668. pmid:6985477
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref19] 19. Gaudry MJ, Storz JF, Butts GT, Campbell KL, Hoffmann FG. Repeated Evolution of Chimeric Fusion Genes in the β-Globin Gene Family of Laurasiatherian Mammals. Genome Biology and Evolution. 2014.
View Article
Google Scholar

[73] View Article

[74] Google Scholar

[ref20] 20. Gibbs MJ, Armstrong JS, Gibbs AJ. Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000;16: 573–582. pmid:11038328
View Article
PubMed/NCBI
Google Scholar

[76] View Article

[77] PubMed/NCBI

[78] Google Scholar

[ref21] 21. Gil M, Zanetti MS, Zoller S, Anisimova M. CodonPhyML: Fast Maximum Likelihood Phylogeny Estimation under Codon Substitution Models. Molecular Biology and Evolution. 2013.
View Article
Google Scholar

[80] View Article

[81] Google Scholar

[ref22] 22. Goldman N, Yang Z. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994;11: 725–736. pmid:7968486
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref23] 23. Goodman M, Czelusniak J, Koop BF, Tagle DA, Slightom JL. Globins: a case study in molecular phylogeny. Cold Spring Harb Symp Quant Biol. 1987;52: 875–890. pmid:3454296
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref24] 24. Goodman M, Koop BF, Czelusniak J, Weiss ML. The eta-globin gene. Its long evolutionary history in the beta-globin gene family of mammals. J Mol Biol. 1984;180: 803–823. pmid:6527390
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref25] 25. Goodman M, Moore GW, Matsuda G. Darwinian evolution in the genealogy of haemoglobin. Nature. 1975;253: 603–608. pmid:1089897
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref26] 26. Hardies SC, Edgell MH, Hutchison CA 3rd. Evolution of the mammalian beta-globin gene cluster. J Biol Chem. 1984;259: 3748–3756. pmid:6706976
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref27] 27. Hardison R. Evolution of hemoglobin and its genes. Cold Spring Harb Perspect Med. 2012;2: a011627. pmid:23209182
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref28] 28. Hardison R. Hemoglobins from bacteria to man: evolution of different patterns of gene expression. J Exp Biol. 1998;201: 1099–1117. pmid:9510523
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref29] 29. Hardison R. Organization, evolution and regulation of the globin genes. In: Steinberg MH, Forget BG, Higgs DR, Nagel RL, editors. Disorders of Hemoglobin: Genetics, Pathophysiology, and Clinical Management. Cambridge: Cambridge University Press; 2001. pmid:12779271

[ref30] 30. Hardison RC. 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. Mol Biol Evol. 1984;1: 390–410. pmid:6599973
View Article
PubMed/NCBI
Google Scholar

[112] View Article

[113] PubMed/NCBI

[114] Google Scholar

[ref31] 31. Hardison RC, Margot JB. Rabbit globin pseudogene psi beta 2 is a hybrid of delta- and beta-globin gene sequences. Mol Biol Evol. 1984;1: 302–316. pmid:6599969
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref32] 32. Hedges SB, Dudley J, Kumar S. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 2006;22: 2971–2972. pmid:17021158
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref33] 33. Hoffmann FG, Opazo JC, Storz JF. New genes originated via multiple recombinational pathways in the beta-globin gene family of rodents. Mol Biol Evol. 2008;25: 2589–2600. pmid:18780876
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref34] 34. Hoffmann FG, Storz JF. The alphaD-globin gene originated via duplication of an embryonic alpha-like globin gene in the ancestor of tetrapod vertebrates. Mol Biol Evol. 2007;24: 1982–1990. pmid:17586601
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref35] 35. Hoffmann FG, Storz JF, Gorr TA, Opazo JC. Lineage-specific patterns of functional diversification in the alpha- and beta-globin gene families of tetrapod vertebrates. Mol Biol Evol. 2010;27: 1126–1138. pmid:20047955
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref36] 36. Holwerda S, de Laat W. Chromatin loops, gene positioning, and gene expression. Front Genet. 2012;3: 217. pmid:23087710
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref37] 37. Jeffreys AJ, Barrie PA, Harris S, Fawcett DH, Nugent ZJ, Boyd AC. Isolation and sequence analysis of a hybrid delta-globin pseudogene from the brown lemur. J Mol Biol. 1982;156: 487–503. pmid:6214636
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref38] 38. Johnson RM, Gumucio D, Goodman M. Globin gene switching in primates. Comp Biochem Physiol A Mol Integr Physiol. 2002;133: 877–883. pmid:12443943
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref39] 39. Keys JR, Tallack MR, Zhan Y, Papathanasiou P, Goodnow CC, Gaensler KM, et al. A mechanism for Ikaros regulation of human globin gene switching. Br J Haematol. 2008;141: 398–406. pmid:18318763
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref40] 40. Koop BF, Siemieniak D, Slightom JL, Goodman M, Dunbar J, Wright PC, et al. Tarsius delta- and beta-globin genes: conversions, evolution, and systematic implications. J Biol Chem. 1989;264: 68–79. pmid:2491855
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref41] 41. Lanfear R, Calcott B, Ho SY, Guindon S. Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29: 1695–1701. pmid:22319168
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref42] 42. Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000;16: 562–563. pmid:10980155
View Article
PubMed/NCBI
Google Scholar

[160] View Article

[161] PubMed/NCBI

[162] Google Scholar

[ref43] 43. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26: 2462–2463. pmid:20798170
View Article
PubMed/NCBI
Google Scholar

[164] View Article

[165] PubMed/NCBI

[166] Google Scholar

[ref44] 44. Martin DP, Posada D, Crandall KA, Williamson C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses. 2005;21: 98–102. pmid:15665649
View Article
PubMed/NCBI
Google Scholar

[168] View Article

[169] PubMed/NCBI

[170] Google Scholar

[ref45] 45. Martin SL, Vincent KA, Wilson AC. Rise and fall of the delta globin gene. J Mol Biol. 1983;164: 513–528. pmid:6188843
View Article
PubMed/NCBI
Google Scholar

[172] View Article

[173] PubMed/NCBI

[174] Google Scholar

[ref46] 46. Martin SL, Zimmer EA, Kan YW, Wilson AC. Silent delta-globin gene in Old World monkeys. Proc Natl Acad Sci U S A. 1980;77: 3563–3566. pmid:6251467
View Article
PubMed/NCBI
Google Scholar

[176] View Article

[177] PubMed/NCBI

[178] Google Scholar

[ref47] 47. Matsuda M, Sakamoto N, Fukumaki Y. Delta-thalassemia caused by disruption of the site for an erythroid-specific transcription factor, GATA-1, in the delta-globin gene promoter. Blood. 1992;80: 1347–1351. pmid:1515647
View Article
PubMed/NCBI
Google Scholar

[180] View Article

[181] PubMed/NCBI

[182] Google Scholar

[ref48] 48. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees; 2010. pp. 1–8.
View Article
Google Scholar

[184] View Article

[185] Google Scholar

[ref49] 49. Moleirinho A, Seixas S, Lopes AM, Bento C, Prata MJ, Amorim A. Evolutionary constraints in the beta-globin cluster: the signature of purifying selection at the delta-globin (HBD) locus and its role in developmental gene regulation. Genome Biol Evol. 2013;5: 559–571. pmid:23431002
View Article
PubMed/NCBI
Google Scholar

[187] View Article

[188] PubMed/NCBI

[189] Google Scholar

[ref50] 50. Noonan JP, Grimwood J, Schmutz J, Dickson M, Myers RM. Gene conversion and the evolution of protocadherin gene cluster diversity. Genome Res. 2004;14: 354–366. pmid:14993203
View Article
PubMed/NCBI
Google Scholar

[191] View Article

[192] PubMed/NCBI

[193] Google Scholar

[ref51] 51. Notredame C, Higgins DG, Heringa J. T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302: 205–217. pmid:10964570
View Article
PubMed/NCBI
Google Scholar

[195] View Article

[196] PubMed/NCBI

[197] Google Scholar

[ref52] 52. Opazo JC, Hoffmann FG, Storz JF. Differential loss of embryonic globin genes during the radiation of placental mammals. Proceedings of the National Academy of Sciences. 2008;105: 12950–12955. pmid:18755893
View Article
PubMed/NCBI
Google Scholar

[199] View Article

[200] PubMed/NCBI

[201] Google Scholar

[ref53] 53. Opazo JC, Hoffmann FG, Storz JF. Genomic evidence for independent origins of β-like globin genes in monotremes and therian mammals. Proceedings of the National Academy of Sciences. 2008;105: 1590–1595. pmid:18216242
View Article
PubMed/NCBI
Google Scholar

[203] View Article

[204] PubMed/NCBI

[205] Google Scholar

[ref54] 54. Opazo JC, Sloan AM, Campbell KL, Storz JF. Origin and Ascendancy of a Chimeric Fusion Gene: The β/δ-Globin Gene of Paenungulate Mammals. Molecular Biology and Evolution. 2009;26: 1469–1478. pmid:19332641
View Article
PubMed/NCBI
Google Scholar

[207] View Article

[208] PubMed/NCBI

[209] Google Scholar

[ref55] 55. Ottolenghi S, Giglioni B, Comi P, Gianni AM, Polli E, Acquaye CT, et al. Globin gene deletion in HPFH, delta (o) beta (o) thalassaemia and Hb Lepore disease. Nature. 1979;278: 654–657. pmid:450068
View Article
PubMed/NCBI
Google Scholar

[211] View Article

[212] PubMed/NCBI

[213] Google Scholar

[ref56] 56. Padidam M, Sawyer S, Fauquet CM. Possible emergence of new geminiviruses by frequent recombination. Virology. 1999;265: 218–225. pmid:10600594
View Article
PubMed/NCBI
Google Scholar

[215] View Article

[216] PubMed/NCBI

[217] Google Scholar

[ref57] 57. Patel VS, Cooper SJ, Deakin JE, Fulton B, Graves T, Warren WC, et al. Platypus globin genes and flanking loci suggest a new insertional model for beta-globin evolution in birds and mammals. BMC Biol. 2008;6: 34. pmid:18657265
View Article
PubMed/NCBI
Google Scholar

[219] View Article

[220] PubMed/NCBI

[221] Google Scholar

[ref58] 58. Petronella N, Drouin G. Purifying selection against gene conversions in the folate receptor genes of primates. Genomics. 2014;103: 40–47. pmid:24184359
View Article
PubMed/NCBI
Google Scholar

[223] View Article

[224] PubMed/NCBI

[225] Google Scholar

[ref59] 59. Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci U S A. 2001;98: 13757–13762. pmid:11717435
View Article
PubMed/NCBI
Google Scholar

[227] View Article

[228] PubMed/NCBI

[229] Google Scholar

[ref60] 60. Prychitko T, Johnson RM, Wildman DE, Gumucio D, Goodman M. The phylogenetic history of New World monkey beta globin reveals a platyrrhine beta to delta gene conversion in the atelid ancestry. Mol Phylogenet Evol. 2005;35: 225–234. pmid:15737593
View Article
PubMed/NCBI
Google Scholar

[231] View Article

[232] PubMed/NCBI

[233] Google Scholar

[ref61] 61. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6; 2014. Available: http://beast.bio.ed.ac.uk/Tracer.

[ref62] 62. Ranney HM, Lam R, Rosenberg G. Some properties of hemoglobin A2. American Journal of Hematology. 1993;42: 107–111. pmid:8416283
View Article
PubMed/NCBI
Google Scholar

[236] View Article

[237] PubMed/NCBI

[238] Google Scholar

[ref63] 63. Ristaldi MS, Casula S, Porcu S, Marongiu MF, Pirastu M, Cao A. Activation of the delta-globin gene by the beta-globin gene CACCC motif. Blood Cells Mol Dis. 1999;25: 193–209. pmid:10575545
View Article
PubMed/NCBI
Google Scholar

[240] View Article

[241] PubMed/NCBI

[242] Google Scholar

[ref64] 64. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61: 539–542. pmid:22357727
View Article
PubMed/NCBI
Google Scholar

[244] View Article

[245] PubMed/NCBI

[246] Google Scholar

[ref65] 65. Rozas J. DNA sequence polymorphism analysis using DnaSP. Methods Mol Biol. 2009;537: 337–350. pmid:19378153
View Article
PubMed/NCBI
Google Scholar

[248] View Article

[249] PubMed/NCBI

[250] Google Scholar

[ref66] 66. Sankaran VG, Orkin SH. The switch from fetal to adult hemoglobin. Cold Spring Harb Perspect Med. 2013;3: a011643. pmid:23209159
View Article
PubMed/NCBI
Google Scholar

[252] View Article

[253] PubMed/NCBI

[254] Google Scholar

[ref67] 67. Sankaran VG, Xu J, Byron R, Greisman HA, Fisher C, Weatherall DJ, et al. A Functional Element Necessary for Fetal Hemoglobin Silencing. New England Journal of Medicine. 2011;365: 807–814. pmid:21879898
View Article
PubMed/NCBI
Google Scholar

[256] View Article

[257] PubMed/NCBI

[258] Google Scholar

[ref68] 68. Sankaran VG, Xu J, Orkin SH. Advances in the understanding of haemoglobin switching. Br J Haematol. 2010;149: 181–194. pmid:20201948
View Article
PubMed/NCBI
Google Scholar

[260] View Article

[261] PubMed/NCBI

[262] Google Scholar

[ref69] 69. Sankaran VG, Xu J, Ragoczy T, Ippolito GC, Walkley CR, Maika SD, et al. Developmental and species-divergent globin switching are driven by BCL11A. Nature. 2009;460: 1093–1097. pmid:19657335
View Article
PubMed/NCBI
Google Scholar

[264] View Article

[265] PubMed/NCBI

[266] Google Scholar

[ref70] 70. Sanyal A, Lajoie BR, Jain G, Dekker J. The long-range interaction landscape of gene promoters. Nature. 2012;489: 109–113. pmid:22955621
View Article
PubMed/NCBI
Google Scholar

[268] View Article

[269] PubMed/NCBI

[270] Google Scholar

[ref71] 71. Schechter AN. Hemoglobin research and the origins of molecular medicine. 3927–3938 p; 2008
View Article
Google Scholar

[272] View Article

[273] Google Scholar

[ref72] 72. Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34: 126–129. pmid:1556748
View Article
PubMed/NCBI
Google Scholar

[275] View Article

[276] PubMed/NCBI

[277] Google Scholar

[ref73] 73. Song G, Hsu CH, Riemer C, Zhang Y, Kim HL, Hoffmann F, et al. Conversion events in gene clusters. BMC Evol Biol. 2011;11: 226. pmid:21798034
View Article
PubMed/NCBI
Google Scholar

[279] View Article

[280] PubMed/NCBI

[281] Google Scholar

[ref74] 74. Spritz RA, Giebel LB. The structure and evolution of the spider monkey delta-globin gene. Mol Biol Evol. 1988;5: 21–29. pmid:2833675
View Article
PubMed/NCBI
Google Scholar

[283] View Article

[284] PubMed/NCBI

[285] Google Scholar

[ref75] 75. Stamatakis A. RAxML Version 8: A tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics. 2014.
View Article
Google Scholar

[287] View Article

[288] Google Scholar

[ref76] 76. Steinberg MH, Adams JG 3rd. Hemoglobin A2: origin, evolution, and aftermath. Blood. 1991;78: 2165–2177. pmid:1932737
View Article
PubMed/NCBI
Google Scholar

[290] View Article

[291] PubMed/NCBI

[292] Google Scholar

[ref77] 77. Steinberg MH, et al. Disorders of Hemoglobin. Cambridge University Press; 2009.

[ref78] 78. Storz JF, Opazo JC, Hoffmann FG. Gene duplication, genome duplication, and the functional diversification of vertebrate globins. Mol Phylogenet Evol. 2013;66: 469–478. pmid:22846683
View Article
PubMed/NCBI
Google Scholar

[295] View Article

[296] PubMed/NCBI

[297] Google Scholar

[ref79] 79. Storz JF, Opazo JC, Hoffmann FG. Phylogenetic diversification of the globin gene superfamily in chordates. IUBMB Life. 2011;63: 313–322. pmid:21557448
View Article
PubMed/NCBI
Google Scholar

[299] View Article

[300] PubMed/NCBI

[301] Google Scholar

[ref80] 80. Tagle DA, Slightom JL, Jones RT, Goodman M. Concerted evolution led to high expression of a prosimian primate delta globin gene locus. J Biol Chem. 1991;266: 7469–7480. pmid:2019578
View Article
PubMed/NCBI
Google Scholar

[303] View Article

[304] PubMed/NCBI

[305] Google Scholar

[ref81] 81. Tang DC, Ebb D, Hardison RC, Rodgers GP. Restoration of the CCAAT box or insertion of the CACCC motif activates [corrected] delta-globin gene expression. Blood. 1997;90: 421–427. pmid:9207479
View Article
PubMed/NCBI
Google Scholar

[307] View Article

[308] PubMed/NCBI

[309] Google Scholar

[ref82] 82. Tang DC, Rodgers GP. Activation of the human delta-globin gene promoter in primary adult erythroid cells. Br J Haematol. 1998;103: 835–838. pmid:9858241
View Article
PubMed/NCBI
Google Scholar

[311] View Article

[312] PubMed/NCBI

[313] Google Scholar

[ref83] 83. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22: 4673–4680. pmid:7984417
View Article
PubMed/NCBI
Google Scholar

[315] View Article

[316] PubMed/NCBI

[317] Google Scholar

[ref84] 84. Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W. Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell. 2002;10: 1453–1465. pmid:12504019
View Article
PubMed/NCBI
Google Scholar

[319] View Article

[320] PubMed/NCBI

[321] Google Scholar

[ref85] 85. Vakoc CR, Letting DL, Gheldof N, Sawado T, Bender MA, Groudine M, et al. Proximity among distant regulatory elements at the beta-globin locus requires GATA-1 and FOG-1. Mol Cell. 2005;17: 453–462. pmid:15694345
View Article
PubMed/NCBI
Google Scholar

[323] View Article

[324] PubMed/NCBI

[325] Google Scholar

[ref86] 86. Vincent KA, Wilson AC. Evolution and transcription of old world monkey globin genes. J Mol Biol. 1989;207: 465–480. pmid:2760921
View Article
PubMed/NCBI
Google Scholar

[327] View Article

[328] PubMed/NCBI

[329] Google Scholar

[ref87] 87. Webster MT, Clegg JB, Harding RM. Common 5' beta-globin RFLP haplotypes harbour a surprising level of ancestral sequence mosaicism. Hum Genet. 2003;113: 123–139. pmid:12736816
View Article
PubMed/NCBI
Google Scholar

[331] View Article

[332] PubMed/NCBI

[333] Google Scholar

[ref88] 88. Yang Z. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998;15: 568–573. pmid:9580986
View Article
PubMed/NCBI
Google Scholar

[335] View Article

[336] PubMed/NCBI

[337] Google Scholar

[ref89] 89. Yang Z. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Molecular Biology and Evolution. 2007;24: 1586–1591. pmid:17483113
View Article
PubMed/NCBI
Google Scholar

[339] View Article

[340] PubMed/NCBI

[341] Google Scholar

[ref90] 90. Yang Z, Nielsen R, Goldman N, Pedersen AM. Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics. 2000;155: 431–449. pmid:10790415
View Article
PubMed/NCBI
Google Scholar

[343] View Article

[344] PubMed/NCBI

[345] Google Scholar

[ref91] 91. Zhang J. Evolution by gene duplication: an update. Trends in Ecology & Evolution. 2003;18: 292–298.
View Article
Google Scholar

[347] View Article

[348] Google Scholar

[ref92] 92. Zhao Z, Hewett-Emmett D, Li WH. Frequent gene conversion between human red and green opsin genes. J Mol Evol. 1998;46: 494–496. pmid:9541545
View Article
PubMed/NCBI
Google Scholar

[350] View Article

[351] PubMed/NCBI

[352] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

DNA Sequence data and gene identification

Evolutionary analysis

Search for open reading frames and promoter analysis

Evolutionary rate estimates and selection tests

Results

Evolutionary history of the two adult β-like genes in mammals

Analysis of gene structure

Recombination events in primates

Evolutionary rates and functional constraints in primates

Discussion

Supporting Information

S1 Fig. Sequence alignment for eutherian HBB proteins.

S2 Fig. Sequence alignment of eutherian HBD open reading frame.

S3 Fig. Sequence alignment of eutherian HBB and HBD promoters.

S4 Fig. Phylograms depicting relationships among adult β-like genes in mammals.

S5 Fig. Maximum Likelihood phylograms depicting relationships among β-like genes of Anthropoids.

S1 Table. Listing of HBB and HBD sequences used for phylogeny reconstructions.

S2 Table. Evolutionary rates based on Jukes Cantor distance.

Author Contributions

References