Molecular history of gene conversions in the primate fetal gamma-globin genes. Nucleotide sequences from the common gibbon, Hylobates lar

Comparative and phylogenetic analyses of homologous sequences from closely related species reveal genetic events which have happened in the past and thus provide considerable insight into molecular genetic processes. One such process which has been especially important in the evolution of multigene families is gene conversion. The fetal gamma 1 and gamma 2-globin genes of catarrhine primates (humans, apes, and Old World monkeys) underwent numerous gene conversion events after they arose from a gene duplication event 25-35 million years ago. By including the gamma 1- and gamma 2-globin gene sequences from the common gibbon, Hylobates lar, the present work expands the gamma-globin data set to represent all major groups of hominoid primates. A computer-assisted algorithm is introduced which reveals converted DNA segments and provides results very similar to those obtained by site-by-site evolutionary reconstruction. Both methods provide strong evidence for at least 14 different converted stretches in catarrhine primates as well as five conversions in ancestral lineages. Features of gene conversions generalized from this molecular history are 1) conversions are restricted to regions maintaining high degrees of sequence similarity, 2) one gene may dominate in converting another gene, 3) sequences involved in conversions may accumulate changes more rapidly than expected, and 4) certain elements, such as polypurine/polypyrimidine [Y)n) and (TG)n elements, appear to be hotspots for initiating or terminating conversion events.     

The phylogeny of human globin genes investigated by the maximum parsimony method

Gene phylogenetic trees were constructed by the maximum parsimony method for various sets of ninety six globin chain amino acid sequences spanning plant and animal kingdoms. The method, executed by several computer programs, constructed ancestor and descendant globin messengers on tree topologies which required the least number of nucleotide replacements to account for the evolution of the globins. The human myoglobin-hemoglobin divergence was traced to a gene duplication which occurred either in the first vertebrates or earlier yet in the common ancestor of chordates and annelids, the alpha-beta divergence to a gene duplication in the common ancestor of teleosts and tetrapods, the gamma divergence from typical beta chains to a gene duplication in basal therian mammals, and the delta separation from beta to a duplication in the basal catarrhine primates. Evidence was provided by the globin phylogenies for the hominoid affinities of the gibbon and the close phyletic relationship of the African apes to man. Over the period of teleos-tetrapod divergence the globin messengers evolved at an average rate of 18.5 nucleotide replacements per 100 codons per 10⁸ years, a faster rate than most previous estimates. Very fast and very slow rates were encountered in different globin lineages and at different stages of descent, reducing the effectiveness of globins as molecular clocks. Rates increased with gene duplication and decreased after selection discovered useful specializations in the products of genes which had previously been freer to accept mutations. The early eutherian radiation was characterized by rapid rates of globin evolution, but the later hominoid radiation by extremely slow rates. This pattern was related to more complicated grades of internal organization evolving in human ancestors. The types of nucleotide replacements in the globin messengers over the long course of globin evolution did not seem indicative of any special mutational mechanisms.     

Tarsius delta- and beta-globin genes: conversions, evolution, and systematic implications

Comparisons between duplicated genes have shown that gene conversions play an important role in the evolution of multigene families. Previous comparisons have documented in the recently duplicated gamma-fetal globin genes of catarrhine primates, over 15 separate conversions affecting extensive stretches of coding and noncoding sequences. In the present study, delta- and beta- globin genes from a lower primate Tarsius syrichta, and the delta-globin gene of the Asian great ape, Pongo pygmaeus, have been isolated and sequenced. Comparisons of these sequences with other primate delta and beta sequences confirmed a previously reported conversion in an anthropoid ancestor and revealed additional conversions in basal primate, stem haplorhine, tarsier, and early lemur lineages. Conversions found between primate delta- and beta-globin genes contrast with those found in the gamma-genes in that delta-beta conversions appear much less frequently and are more restricted to regions conserved by selection (i.e. coding and 5'-regulatory sequences). These differences indicate that soon after a duplication occurs, conversions can be quite frequent and encompass extensive portions of the duplicated region. With time, sequence differences accumulate, particularly in noncoding regions, and limit both the frequency and size of the conversions. Sequences conserved by selection accumulate differences more slowly and are therefore subject to gene conversions for a longer period of time. Both unconverted and converted sequences were consistent in supporting the placement of tarsier with anthropoids.     

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

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

Concerted evolution of the primate immunoglobulin alpha-gene through gene conversion.

We determined four nucleotide sequences of the hominoid immunoglobulin alpha (C alpha) genes (chimpanzee C alpha 2, gorilla C alpha 2, and gibbon C alpha 1 and C alpha 2 genes), which made possible the examination of gene conversions in all hominoid C alpha genes. The following three methods were used to detect gene conversions: 1) phenetic tree construction; 2) detection of a DNA segment with extremely low variability between duplicated C alpha genes; and 3) a site by site search of shared nucleotide changes between duplicated C alpha genes. Results obtained from method 1 indicated a concerted evolution of the duplicated C alpha genes in the human, chimpanzee, gorilla, and gibbon lineages, while results obtained from method 2 suggested gene conversions in the human, gorilla, and gibbon C alpha genes. With method 3 we identified clusters of shared nucleotide changes between duplicated C alpha genes in human, chimpanzee, gorilla, and gibbon lineages, and in their hypothetical ancestors. In the present study converted regions were identified over the entire C alpha gene region excluding a few sites in the coding region which have escaped from gene conversion. This indicates that gene conversion is a general phenomenon in evolution, that can be clearly observed in non-functional regions.     

Intergenic DNA sequences flanking the pseudo alpha globin genes of human and chimpanzee.

We have determined the sequence of 2400 base pairs upstream from the human pseudo alpha globin (psi alpha) gene, and for comparison, 1100 base pairs of DNA within and upstream from the chimpanzee psi alpha gene. The region upstream from the promoter of the psi alpha gene shows no significant homology to the intergenic regions of the adult alpha 2 and alpha 1 globin genes. The chimpanzee gene has a coding defect in common with the human psi alpha gene, showing that the product of this gene, if any, was inactivated before the divergence of human and chimpanzee. However the chimpanzee gene contains a normal ATG initiation codon in contrast to the human gene which has GTG as the initiation codon. The psi alpha genes of both human and chimpanzee are flanked by the same Alu family member. The structure and position of this repeat have not been altered since the divergence of human and chimpanzee, and it is at least as well conserved as its immediate flanking sequence. Comparing human and chimpanzee, the 300 bp Alu repeat has accumulated only two base substitutions and one length mutation; the adjacent 300 bp flanking region has accumulated five base substitutions and twelve length mutations.     

Nucleotide sequence analysis of the lemur beta-globin gene family: evidence for major rate fluctuations in globin polypeptide evolution

Lemur beta-related globin genes have been isolated and sequenced. Orthology of prosimian and human epsilon-, gamma-, and beta-related globin genes was established by dot-matrix analysis. All of these lemur globin genes potentially encode functional beta-related globin polypeptides, though precisely when the gamma-globin gene is expressed remains unknown. The organization of the 18-kb brown lemur beta-globin gene cluster (5' epsilon-gamma-[psi eta-delta]-beta 3') is consistent with its evolution by contraction via unequal crossing-over from the putative ancestral mammalian beta-globin gene cluster (5' epsilon-gamma- eta-delta-beta 3'). The dwarf lemur nonadult globin genes are arranged as in the brown lemur. Similar levels of synonymous (silent) nucleotide substitutions and noncoding DNA sequence differences have accumulated between species in all of these genes, suggesting a uniform rate of noncoding DNA divergence throughout primate beta-globin gene clusters. These differences are comparable with those observed in the nonfunctional psi eta pseudogene and have therefore accumulated at the presumably maximal neutral rate. In contrast, nonsynonymous (replacement) nucleotide substitutions show a significant heterogeneity in distribution for both the same gene in different lineages and different genes in the same lineage. These major fluctuations in replacement but not silent substitution rates cannot be attributed to changes in mutation rate, suggesting that changes in the rate of globin polypeptide evolution in primates is not governed solely by variable mutation rates.      

Nucleotide sequence of the goat embryonic alpha globin gene (zeta) and linkage and evolutionary analysis of the complete alpha globin cluster

In previous studies we identified and sequenced clones containing two adult alpha globin genes of the goat. Additional studies have revealed the presence of an embryonic alpha globin gene termed zeta. Sequence analysis of the gene shows that it is the largest mammalian or avian globin gene cloned to date. Its unusual size is mainly due to a 14 base-pair tandem repeat sequence in its first intron. A similar sequence is also found in the first intron of the human zeta gene. The goat zeta coding sequence differs greatly from that of the adult alpha, particularly at amino acid position 38, where it codes for the amino acid replacement of Gln for Thr. This change may confer a higher intrinsic O2 affinity on the zeta globin protein, ensuring a sufficient O2 supply for the developing goat embryo. The cloning and sequencing of this gene completes the alpha globin locus of the goat, composed of three genes in the following order 5'-zeta-I alpha-II alpha-3'. Evolutionary comparisons of the goat alpha locus with other amphibian, avian and mammalian loci reveal several interesting features. Statistical analysis confirms the hypothesis that the embryonic alpha gene is much older (400 million years) than the embryonic beta gene (200 million years), and that it is descended from a primordial gene, whose present-day counterpart is the Xenopus larval alpha globin gene. Our results also suggest that after the divergence of the avian line, the alpha A gene converted the alpha D gene during the evolution of the pre-mammalian line. The alpha D globin gene remains unconverted in the avian line, potentially because of insertion/deletion sequences that may prevent any gene conversion event. The divergence rates of specific globin genes have been analyzed and found to form an essentially straight line, in agreement with the neutralist view of evolution.     

Sequences from the 5′ flanking region of the ϵ‐globin gene support the relationship of Callicebus with the pitheciins

The purpose of this study was to determine nucleotide sequences from the 5′ flanking region of the ϵ‐globin gene of selected platyrrhine primates and to analyze the data for phylogenetic information and estimated times of divergence. We report new sequence data for two species of New World monkeys, Callicebus torquatus and Pithecia irrorata. We analyzed these data in conjunction with homologous sequences from other primate species. The data support the hypothesis that the titi monkeys (Callicebus) and seed predators (Tribe Pitheciini) form a clade (Subfamily Pitheciinae), and also provide limited support for that subfamily being allied with the atelines. We also present estimated dates of divergence for the Callicebus and pitheciin lineages. Am. J. Primatol. 48:69–75, 1999. © 1999 Wiley‐Liss, Inc.     

Species status and phylogeography of two closely related Coptolabrus species (Coleoptera: Carabidae) in South Korea inferred from mitochondrial and nuclear gene sequences

We investigated the species status and intraspecific phylogeography in South Korea of two ground beetle species, Coptolabrus jankowskii and Coptolabrus smaragdinus (Coleoptera: Carabidae), using statistical parsimony networks and nested clade analyses based on sequences from the mitochondrial cytochrome oxidase subunit I (COI) and nuclear phosphoenolpyruvate carboxykinase (PepCK) and wingless (Wg) genes. Although traditional parsimony tree construction generally failed to resolve interspecific relationships and construct biologically meaningful genealogies, analysis using statistical parsimony networks yielded statistically significant inter- and intraspecific genealogical structures. We found that although these two species represent a notable case of trans-species polymorphisms in both mitochondrial and nuclear gene sequences, their status as separate species was evidenced by the nonrandom association between species and nested clades at various nesting levels. The exceptional occurrence of shared identical or very similar COI sequences was considered to be the result of introgressive hybridization. In addition, range expansion and fragmentation events across the Korean Peninsula and adjacent islands were inferred from nested clade phylogeographical analyses. The COI gene revealed the geographical divergence of major eastern and western clades and historical biogeographical events within each major clade, whereas the nuclear PepCK gene, which did not reveal corresponding east-west clades, indicated past fragmentation and range expansion across wide areas that may have been the result of older biogeographical events. Thus, phylogeographical inferences drawn from analyses of mitochondrial and nuclear genes can reveal different and potentially complementary information about phylogeographical processes.     

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

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

Genomic and phylogenic comparisons of the alpha-globin and beta-globin intergenic sequences between zebra fish (Danio rerio) and six closely related Cyprinindae species

The alpha-globin and beta-globin genes of the zebrafish are tightly linked on the same chromosome in a 3'-5' and 5'-3' configuration, respectively. Although the location of the controlling sequences has been mapped to the intergenic region, analysis determining the uniqueness of this unusual arrangement to zebrafish has not been undertaken. To explore this, we isolated, sequenced, and phylogenetically analyzed the intergenic region between globin gene families of seven Cyprinindae species including zebrafish. These species were grouped into an in group (immediate relatives, not so distant relatives), and an out group (distant relative). Cellulose acetate electrophoresis of hemoglobin (Hb) detected multiple variants in each species, but a band with electrophoretic mobility (EM) of 6.7 x 10(-5) cm(2).volt(-1).sec(-1) was shared between species. Polymerase chain reaction (PCR) amplification of the intergenic globin gene region also detected a 1.0-kb fragment that was repeated in the in group and a 1.2-kb fragment in the out group. Sequence comparison confirmed that the genetic orientation and controlling sequences location were conserved throughout this region in all seven species. This phylogenic footprinting indicated that the configuration was not exclusive to zebrafish. To confirm sequence alignment, maximum parsimony phylogenic analysis, was performed. Only one member of that group the giant danio, was not closely clustered, being located almost equidistance between the immediate relative and the species of the other clusters. This may represent an ancestral configuration prior to transposition of the alpha globin and beta-globin genes families to nonsynteny.     

Phylogenetic Analysis of the Hoplolaiminae Inferred from Combined D2 and D3 Expansion Segments of 28S rDNA

DNA sequences of the D2-D3 expansion segments of the 28S gene of ribosomal DNA from 23 taxa of the subfamily Hoplolaiminae were obtained and aligned to infer phylogenetic relationships. The D2 and D3 expansion regions are G-C rich (59.2%), with up to 20.7% genetic divergence between Scutellonema brachyurum and Hoplolaimus concaudajuvencus. Molecular phylogenetic analysis using maximum likelihood and maximum parsimony was conducted using the D2-D3 sequence data. Of 558 characters, 254 characters (45.5%) were variable and 198 characters (35.4%) were parsimony informative. All phylogenetic methods produced a similar topology with two distinct clades: One clade consists of all Hoplolaimus species while the other clade consists of the rest of the studied Hoplolaiminae genera. This result suggests that Hoplolaimus is monophyletic. Another clade consisted of Aorolaimus, Helicotylenchus, Rotylenchus, and Scutellonema species. Phylogenetic analysis using the outgroup species Globodera rostocheinsis suggests that Hoplolaiminae is paraphyletic. In this study, the D2-D3 region had levels of DNA sequence divergence sufficient for phylogenetic analysis and delimiting species of Hoplolaiminae.     

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

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

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25–50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events.