Rhesus fetal globin genes. Concerted gene evolution in the descent of higher primates

A combined top-down and bottom-up MS approach for the characterization of hemoglobin variants in Rhesus monkeys

Sickle cell disease is caused by one of the 1200 known hemoglobin variations. A single-point mutation β6(A3)Glu→Val leads to sickling of red blood cells, which in turn causes a lack of oxygen supply to tissue and organs. Although sickle cell disease is well understood, treatment options are currently underdeveloped. The only Food and Drug Administration-approved drug is hydroxyurea, an inducer of fetal γ-hemoglobin, which is known to have a higher oxygen affinity than adult hemoglobins and thus alleviates symptoms. In the search for better cures, Rhesus monkeys (Macaca mulatta) serve as models for monitoring success of induction of fetal γ-hemoglobins and with recent advances in proteomics, MS has become the leading technique to determine globin expression. Similar to humans, Rhesus monkeys possess hemoglobin variants that have not been sufficiently characterized to initiate such a study. Therefore, we developed a combined bottom-up and top-down approach to identify and characterize novel hemoglobin variants of the umbilical cord blood of Rhesus monkeys. A total of four different variants were studied: α, β, γ1 and γ2. A new α- and β-hemoglobin variant was identified, and the two previously hypothesized γ-hemoglobins were identified. In addition, glutathionylation of both γ-hemoglobin variants at their cysteines has been characterized. The combined approach outperformed either bottom-up or top-down alone and can be used for characterization of unknown hemoglobin variants and their PTMs.

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.

Humans and old world monkeys have similar patterns of fetal globin expression

The expression of epsilon- and gamma-globin mRNA and protein has been determined in three Old World monkey species (Macaca mulatta, Macaca nemestrina, and Cercopithecus aethiops). Using RT-PCR with primers for epsilon- and gamma-globin, both mRNAs were detected in early fetal stages, whereas at 128 days (85% of full term), only gamma was expressed. High-performance liquid chromatography was used for separation and quantitation, and matrix-assisted laser desorption/ionization mass spectrometry was used for identification of globin polypeptides. An alpha-globin polymorphism was observed in all of the species examined. During fetal life, gamma-globin was the predominant expressed beta-type globin. The red blood cells of infants still contained substantial amounts of gamma-globin, which declined to negligible levels in 14 weeks as beta-globin expression reached adult values. The ratio of gamma1- to gamma2-globins (equivalent to Ggamma/Agamma in humans) was approximately 2.5, similar to the Ggamma/Agamma ratio observed in humans. Thus, gamma-globin gene expression in these Old World monkeys species has three features in common with human expression: expression of both duplicated gamma genes, the relative preponderance of gamma1 over gamma2 expression, and the delay of the switch from gamma- to beta-globin until the perinatal period. Thus, the catarrhines seem to share a common pattern of developmental switching in the beta-globin gene cluster, which is distinct from the timing of expression in either prosimians or the New World monkeys. Our results indicate that an Old World monkey, such as Rhesus, could serve as a model organism (resembling humans) for experimentally investigating globin gene expression patterns during the embryonic, fetal, and postnatal stages.

Phylogenetic relationships among Brazilian howler monkeys, genus Alouatta (Platyrrhini, Atelidae), based on g1-globin pseudogene sequences

The genus Alouatta (howler monkeys) is the most widely distributed of New World primates, and has been arranged in three species groups: the Central American Alouatta palliata group and the South American Alouatta seniculus and Alouatta caraya groups. While the latter is monotypic, the A. seniculus group encompasses at least three species (A. seniculus, A. belzebul and A. fusca). In the present study, approximately 600 base pairs of the g1-globin pseudogene were sequenced in the four Brazilian species (A. seniculus, A. belzebul, A. fusca and A. caraya). Maximum parsimony and maximum likelihood methods yielded phylogenetic trees with the same arrangement: {A. caraya [A. seniculus (A. fusca, A. belzebul)]}. The most parsimonious tree had bootstrap values greater than 82% for all groupings, and strength of grouping values of at least 2, supporting the sister clade of A. fusca and A. belzebul. The study also confirmed the presence of a 150-base pair Alu insertion element and a 1.8-kb deletion in the g1-globin pseudogene in A. fusca, features found previously in the remaining three species. The cladistic classification based on molecular data agrees with those of morphological studies, with the monospecific A. caraya group being clearly differentiated from the A. seniculus group.

Molecular phylogeny of ateline new world monkeys (Platyrrhini, atelinae) based on gamma-globin gene sequences: evidence that brachyteles is the sister group of lagothrix

Nucleotide sequences, each spanning approximately 7 kb of the contiguous gamma1 and gamma2 globin genomic loci, were determined for seven species representing all extant genera (Ateles, Lagothrix, Brachyteles, and Alouatta) of the New World monkey subfamily Atelinae. After aligning these seven ateline sequences with outgroup sequences from several other primate (non-ateline) genera, they were analyzed by maximum parsimony, maximum likelihood, and neighbor-joining algorithms. All three analyzes estimated the same phylogenetic relationships: [Alouatta [Ateles (Brachyteles, Lagothrix)]]. Brachyteles and Lagothrix are sister-groups supported by 100% of bootstrap replications in the parsimony analyses. Ateles joins this clade, followed by the basal genus Alouatta; these joinings were strongly supported, again with 100% bootstrap values. This cladistic pattern for the four ateline genera is congruent with that obtained in previous studies utilizing epsilon-globin, IRBP, and G6PD nuclear genomic sequences as well as mitochondrial COII sequences. Because the number of aligned nucleotide positions is much larger in the present datasetoff than in any of these other datasets, much stronger support was obtained for the cladistic classification that divides subfamily Atelinae into tribes Alouattini (Alouatta) and Atelini, while the latter divides into subtribes Atelina (Ateles) and Brachytelina (Brachyteles and Lagothrix).

Dynamics of regulatory evolution in primate beta-globin gene clusters: cis-mediated acquisition of simian gamma fetal expression patterns

Phylogenetic reconstructions by parsimony were carried out on an enlarged body of primate gamma1 and gamma2-globin sequences. The results confirm that gamma1 and gamma2 arose from a tandem duplication in an ancient simian lineage ancestral to both platyrrhines (New World monkeys) and catarrhines (Old World monkeys and hominoids). Gene conversions between the two gamma homologs were frequent over the gamma gene proper but less frequent over the 5' flanking and very infrequent over the 3' flanking regions. The ancient platyrrhine conversion in the most distal 5' flanking region had the polarity of gamma2-->gamma1. Recent platyrrhine conversions between 5' regulatory sequences were very infrequent, in striking contrast to catarrhines which have large, uninterrupted stretches of converted 5' regulatory sequences. Comparisons of reconstructed ancestral primate and simian gamma promoter sequences revealed an accumulation of 21 nucleotide substitutions concentrated in or near cis-elements that may have mediated the change from embryonic to fetal gamma expression. Almost all 21 substitutions were retained in the lineages leading to functional gamma genes of extant catarrhines (both gamma1 and gamma2) and platyrrhines (most often gamma2). Fewer of these simian specific substitutions were retained in the platyrrhine gamma1 genes and new mutations occurred more often in the platyrrhine gamma1 than gamma2 promoters.

Evolution and selection of primate T cell antigen receptor BV8 gene subfamily

The set of potential T cell receptor specificities is highly diverse. The relative contributions of T cell receptor (TCR) V beta gene segment polymorphisms, duplications, deletions, and gene conversions to this final T cell receptor protein diversity are unknown. To study these mechanisms, we sequenced and compared closely related primate TCR gene segments from BV8S1, S2, and S5. Interspecies comparisons show that these gene segments have sustained multiple duplication, gene conversion, and deletion events during the last 35 million years of anthropoid primate evolution. BV8 coding sequences are generally conserved with respect to their flanking noncoding sequences, but we find no evidence for positive or negative selection in sequences coding for the first two putative complementarity-determining (ligand-binding) regions. Sequences of TCRBV8 gene segments from unrelated humans demonstrate no nonsynonymous substitutions in nonleader regions of either the BV8S1 or S2 gene segments. We conclude that gene duplication, deletion, and conversion mechanism contribute in a substantial way to the overall diversity of the TCRBV8 gene segment repertoire in primate evolution and that germline substitutions and consequent polymorphisms in CDRs 1 and 2 of these gene segments probably do not play an active role in generating TCR beta chain protein variation.

Fetal globin expression in New World monkeys

Reverse phase chromatography of the globin chains of adult, newborn, and fetal erythrocytes from three species of New World monkeys (Cebus apella, Aotus azarae, and Callithrix jacchus) representing three of the seven platyrrhine clades showed that gamma-globin expression was fetal in these animals. The globins were identified by a combination of chemical sequencing and mass spectrometric analysis. Since gamma-globin expression is fetal in the other major simian branch, the catarrhines, but embryonic in prosimian primates and nonprimate placental mammals, the evolution of fetal recruitment can now be assigned to the period between the simian-prosimian divergence (55 million years ago) and the platyrrhine-catarrhine divergence (35 million years ago). The gamma-globin gene underwent tandem duplication during the same evolutionary epoch, in accord with a model that suggests that the downstream duplicated gamma-gene (gamma2) was free to acquire the mutations necessary for fetal recruitment. Mass spectrometric analysis of tryptic digests of the gamma-globins verified the amino acid sequences deduced from genomic sequencing. Detailed analysis of high performance liquid chromatography and matrix-assisted laser desorption/ionization mass spectrometry data showed that gamma2-globin in Cebus was expressed to a far greater extent than gamma1-globin, supporting inferences drawn from a study of the promoter sequences. A "pre-gamma"-globin was observed in C. apella and shown to be primarily the glutathionyl adduct. The other species, A. azarae and C. jacchus, also express only one gamma-globin polypeptide. This work provides biochemical evidence of an evolutionary trend in the platyrrhines to alter the duplicated gamma-globin gene locus so that only one gamma-globin polypeptide is expressed.

Variation in short tandem repeat sequences--a survey of twelve microsatellite loci for use as forensic identification markers

Alleles at 12 Short Tandem Repeat loci have been sequenced to investigate candidate loci for a multiplex Short Tandem Repeat system for forensic identification, and for single-locus amplification of Short Tandem Repeat loci. Variation from the consensus sequence was found at 6 loci, while one locus, D21S11, was found to be complex in sequence. The presence of non-consensus alleles does not rule out loci for inclusion as forensic identification markers, but size differences between alleles of 1 base pair require very precise sizing. We suggest criteria for the suitability of Short Tandem Repeat loci as forensic identification markers, and propose a universal allele nomenclature for simple and compound Short Tandem Repeats. The effect of the repeat unit sequence of the evolution of Short Tandem Repeats is discussed.

Primary sequence, evolution, and repetitive elements of the Gallus gallus (chicken) beta-globin cluster

The DNA sequence of the Gallus gallus (chicken) beta-globin cluster was completed and analyzed. This G + C-rich region is 23.7 kb in length and includes the rho-, beta H-, beta A-, and epsilon-globin genes, the enhancer found between the beta A and epsilon genes, and three upstream DNase I hypersensitive sites. The CpG dinucleotides are nonrandomly distributed, being present at an increased relative frequency near the promoters and upstream hypersensitive sites. The cluster has an unusually low TA dinucleotide frequency. The upstream hypersensitive sites (5'HS1, 5'HS2, and 5'HS3) contain DNA sequence motifs recognized by erythroid transcription factors. However, no significant sequence similarity was found among the upstream hypersensitive sites and the beta A/epsilon enhancer. The G. gallus upstream site sequences were not similar to the upstream sites of the mammalian globin clusters, probably due to the small size of the functional regions and large evolutionary distance between the classes. The avian cluster evolved by gene duplication from an ancestor beta-globin gene, first producing the epsilon and the rho/beta H/beta A ancestor genes, then the rho and the beta H/beta A ancestor genes, and finally the beta H- and beta A-globins. Four probable gene conversions can be documented: beta A to beta H, epsilon to beta H, and rho/epsilon (twice). The cluster shows a massive overrepresentation of a non-LTR retrotransposon, CR1, which accounts for 16% of the DNA. We suggest that the locus is a preferred site for CR1 insertion.

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

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

Fetal recruitment of anthropoid gamma-globin genes. Findings from phylogenetic analyses involving the 5'-flanking sequences of the psi gamma 1 globin gene of spider monkey Ateles geoffroyi

We determined the nucleotide sequence of a 2.5 kb DNA fragment (1 kb is 10(3) base-pairs) that includes exon 1, intron 1 and about 1.4 kb of 5'-flanking DNA of the spider monkey gamma 1-globin pseudogene locus and compared this sequence to its homologous from other primates and rabbit. This region of the gamma 1 locus of spider monkey still retains conserved regulatory elements, suggesting that it became a pseudogene late in New World monkey phylogeny. In the 250 base-pair region immediately 5' from the transcription start site where many known regulatory elements are located, a higher rate of nucleotide substitutions occurred in the ancestral anthropoid (human, ape and monkey) lineage than in the prosimian (galago) lineage, as was also the case for non-synonymous substitutions in the coding region. The opposite pattern was observed for most other non-coding regions and for synonymous substitutions. These substitution patterns correlate with the embryonic-to-fetal transformation of the gamma-globin genes of the ancestral anthropoids. Analysis of the 5'-flanking sequences suggests that 11 gene conversion events have occurred in the anthropoid gamma-gene lineages. In the parts of the 5'-flanking region where no gene conversions have been detected, gamma 2-gene sequences have accumulated more nucleotide changes than gamma 1, which suggests that the gamma 2 gene was the more redundant duplicate that may have accumulated first the nucleotide changes responsible for the anthropoid fetal pattern of gamma-globin gene expression.

Primate evolution at the DNA level and a classification of hominoids

The genetic distances among primate lineages estimated from orthologous noncoding nucleotide sequences of beta-type globin loci and their flanking and intergenic DNA agree closely with the distances (delta T50H values) estimated by cross hybridization of total genomic single-copy DNAs. These DNA distances and the maximum parsimony tree constructed for the nucleotide sequence orthologues depict a branching pattern of primate lineages that is essentially congruent with the picture from phylogenetic analyses of morphological characters. The molecular evidence, however, resolves ambiguities in the morphological picture and provides an objective view of the cladistic position of humans among the primates. The molecular data group humans with chimpanzees in subtribe Hominina, with gorillas in tribe Hominini, orangutans in subfamily Homininae, gibbons in family Hominidae, Old World monkeys in infraorder Catarrhini, New World monkeys in semisuborder Anthropoidea, tarsiers in suborder Haplorhini, and strepsirhines (lemuriforms and lorisiforms) in order Primates. A seeming incongruency between organismal and molecular levels of evolution, namely that morphological evolution appears to have speeded up in higher primates, especially in the lineage to humans, while molecular evolution has slowed down, may have the trivial explanation that relatively small genetic changes may sometimes result in marked phenotypic changes.