Structure and chromosomal arrangement of leghemoglobin genes in kidney bean suggest divergence in soybean leghemoglobin gene loci following tetraploidization

We have determined the structure of one of the leghemoglobin (Lb) genes of Phaseolus vulgaris (kidney bean) and deduced the chromosomal arrangement of leghemoglobin genes by genomic hybridizations with Lb and two other sequences, each specific to the 5' or 3' region of the soybean leghemoglobin loci. By comparing this organization with two other species of legumes, Glycine max (soybean) and G. soja (wild soybean), a phylogeny of leghemoglobin gene loci was traced. The intragenic structure of the kidney bean leghemoglobin gene shows the same intron/exon arrangement as that of soybean leghemoglobin genes and extensive sequence homologies in both coding as well as 5' and 3' non-coding regions. The presence in the kidney bean genome of four leghemoglobin genes suggests that tandem duplications of a single primordial plant globin gene had occurred to generate four leghemoglobin genes in an `Lb-locus' before Glycine and Phaseolus species diverged. Chromosome duplication by tetraploidization in Glycine generated two loci containing four genes each. A large deletion in one of the two four-gene loci in soybean resulted in the generation of the Lbc₂ locus containing two leghemoglobin genes, one functional and another pseudo (LbΨ₂). This pseudogene, unlike that present on the main locus, is represented by only two and a half exons and appears to be truncated. The two other truncated genes (LbT₁ and LbT₂) were probably generated similarly in the genome of Glycine spp. following tetraploidization before the divergence of G. max and G. soja.     

Evolutionary rate variation among vertebrate beta globin genes: implications for dating gene family duplication events

A comprehensive dataset of 62 beta globin gene sequences from various vertebrates was compiled to test the molecular clock and to estimate dates of gene duplications. We found that evolution of the beta globin family of genes is not clock-like, a result that is at odds with the common use of this family as an example of a constant rate of evolution over time. Divergence dates were estimated either with or without assuming the molecular clock, and both analyses produced similar date estimates, which are also in general agreement with estimates reported previously. In addition we report date estimates for seven previously unexamined duplication events within the beta globin family. Despite multiple sources of rate variation, the average rate across the beta globin phylogeny yielded reasonable estimates of divergence dates in most cases. Exceptions were cases of gene conversion, where it appears to have led to underestimates of divergence dates. Our results suggest (i) the major duplications giving rise to the paralogous beta globin genes are associated with significant evolutionary rate variation among gene lineages; and (ii) genes arising from more recent gene duplications (e.g., tandem duplications within lineages) do not appear to differ greatly in rate. We believe this pattern reflects a complex interplay of evolutionary forces where natural selection for diversifying paralogous functions and lineage-specific effects contribute to rate variation on a long-term basis, while gene conversion tends to increase sequence similarity. Gene conversion effects appear to be stronger on recent gene duplicates, as their sequences are highly similar. Lastly, phylogenetic analyses do not support a previous report that avian globins are members of a relic lineage of omega globins.     

The primary structures of two leghemoglobin genes from soybean

We present the complete nucleotide sequences of two leghemoglobin genes isolated from soybean DNA. Both genes contain three intervening sequences which interrupt the two coding sequences in identical positions. The 5' and 3' flanking sequences in both genes contain conserved sequences similar to those found in corresponding positions in other eukaryotic genes. Thus, the general DNA sequence organization of these plant genes is similar to that of other eukaryotic genes.     

Primary structure and evolutionary relationship between the adult α-globin genes and their 5′-flanking regions ofXenopus laevis andXenopus tropicalis

Summary To investigate the evolution of globin genes in the genusXenopus, we have determined the primary structure of the related adult ₁- and _II genes ofX. laevis and of the adult -globin gene ofX. tropicalis, including their 5-flanking regions. All three genes are comprised of three exons and two introns at homologous positions. The exons are highly conserved and code for 141 amino acids. By contrast, the corresponding introns vary in length and show considerable divergence. Comparison of 900 bp of the 5-flanking region revealed that theX. tropicalis gene contains a conserved proximal 310-bp promoter sequence, comprised of the canonical TATA and CCAAT motifs at homologous positions, and five conserved elements in the same order and at similar positions as previously shown for the corresponding genes ofX. laevis. We therefore conclude that these conserved upstream elements may represent regulatory sequences for cell-specific regulation of the adultXenopus globin genes.     

Leghemoglobin-like sequences in the DNA of four actinorhizal plants

Summary A cloned cDNA partial copy of a soybean leghemoglobin mRNA was used to probe genomic DNA of four species of actinorhizal plants. Southern blot hybridization revealed the presence of sequences with homology to the leghemoglobin probe in DNA from Alnus glutinosa, Casuarina glauca, Ceanothus americanus and Elaeagnus pungens. The hybridization patterns of the restriction fragments revealed some fragment size conservation between the DNA of soybean and the DNA of four actinorhizal plants which are taxonomically unrelated to soybean or to each other. The results presented here indicate that globin gene sequences are much more widely distributed in the plant kingdom than has previously been thought. Furthermore, if sequence conservation is actually as high as the restriction fragment patterns suggest, the evolution of the DNA surrounding the globin sequences has been highly constrained.     

A nematode hemoglobin gene contains an intron previously thought to be unique to plants

Summary Hemoglobin genes from plants and animals both have a characteristic chromosomal organization. Plant hemoglobin genes contain a unique intron inserted into the heme-binding domain of exon 2. This intron has not been previously reported in animal globin genes, and its loss was hypothesized to have occurred early in the evolution of hemoglobins. We report here a unique six-intron, seven-exon internally duplicated nematode hemoglobin gene that contains an intron equivalent to the plant central intron in its first repeat. This nematode hemoglobin gene has lost both the central and the normal third intron in its second repeat. The nematode globin also contains a unique intron between its secretory peptide leader sequence and its coding sequence, which is absent in other extracellular invertebrate globin genes. Possible models to explain the head-to-tail duplication of this gene are discussed. Offprint requests to: B. Pohajdak     

The Molecular Evolution of Actin

We have investigated the molecular evolution of plant and nonplant actin genes comparing nucleotide and amino acid sequences of 20 actin genes. Nucleotide changes resulting in amino acid substitutions (replacement substitutions) ranged from 3-7% for all pairwise comparisons of animal actin genes with the following exceptions. Comparisons between higher animal muscle actin gene sequences and comparisons between higher animal cytoplasmic actin gene sequences indicated less than 3% divergence. Comparisons between plant and nonplant actin genes revealed, with two exceptions, 11-15% replacement substitution. In the analysis of plant actins, replacement substitution between soybean actin genes SAc1, SAc3, SAc4 and maize actin gene MAc1 ranged from 8-10%, whereas these members within the soybean actin gene family ranged from 6-9% replacement substitution. The rate of sequence divergence of plant actin sequences appears to be similar to that observed for animal actins. Furthermore, these and other data suggest that the plant actin gene family is ancient and that the families of soybean and maize actin genes have diverged from a single common ancestral plant actin gene that originated long before the divergence of monocots and dicots. The soybean actin multigene family encodes at least three classes of actin. These classes each contain a pair of actin genes that have been designated kappa (SAc1, SAc6), lambda (SAc2, SAc4) and mu (SAc3, SAc7). The three classes of soybean actin are more divergent in nucleotide sequence from one another than higher animal cytoplasmic actin is divergent from muscle actin. The location and distribution of amino acid changes were compared between actin proteins from all sources. A comparison of the hydropathy of all actin sequences, except from Oxytricha, indicated a strong similarity in hydropathic character between all plant and nonplant actins despite the greater number of replacement substitutions in plant actins. These protein sequence comparisons are discussed with respect to the demonstrated and implicated roles of actin in plants and animals, as well as the tissue-specific expression of actin.     

The beta globin gene cluster of the prosimian primate Galago crassicaudatus: nucleotide sequence determination of the 41-kb cluster and comparative sequence analyses.

The nucleotide sequence of the beta globin gene cluster of the prosimian Galago crassicaudatus has been determined. A total sequence spanning 41,101 bp contains and links together previously published sequences of the five galago beta-like globin genes (5'-epsilon-gamma-psi eta-delta-beta-3'). A computer-aided search for middle interspersed repetitive sequences identified 10 LINE (L1) elements, including a 5' truncated repeat that is orthologous to the full-length L1 element found in the human epsilon-gamma intergenic region. SINE elements that were identified included one Alu type I repeat, four Alu type II repeats, and two methionine tRNA-derived Monomer (type III) elements. Alu type II and Monomer sequences are unique to the galago genome. Structural analyses of the cluster sequence reveals that it is relatively A+T rich (about 62%) and regions with high G+C content are associated primarily with globin coding regions. Comparative analyses with the beta globin cluster sequences of human, rabbit, and mouse reveal extensive sequence homologies in their genic regions, but only human, galago, and rabbit sequences share extensive intergenic sequence homologies. Divergence analyses of aligned intergenic and flanking sequences from orthologous human, galago, and rabbit sequences show a gradation in the rate of nucleotide sequence evolution along the cluster where sequences 5' of the epsilon globin gene region show the least sequence divergence and sequences just 5' of the beta globin gene region show the greatest sequence divergence.     

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.      

Aspects of the Origin and Evolution of Non-Vertebrate Hemoglobins

Hemoglobins, found in members of almost all invertebrate phyla,display an extraordinary diversity of form and function. Althoughsome are intracellular with chains and assemblies similar insize to those of vertebrates, others are giant extracellularproteins with masses as large as 8,000 kilodaltons. Two verydifferent strategies have evolved for the stabilization of theselarge molecules. The first is the formation of both intra- andinterchain disulfide bonds that effectively immobilize segmentsof the protein, and the second is the evolution by gene duplicationof multi-domain chains with from two to eighteen myoglobin-like,heme-containing domains in a single polypeptide that may beas large as {small tilde}260 kilodaltons. The genes for vertebrateglobins have a characteristic two-intron, threeexon structure.The gene encoding chain c of the hemoglobin of the earthwormLumbricus terrestris has precisely the same organization andsplice junction positions. This shows that these positions havebeen conserved for at least 600 million years, the estimatedtime of divergence of annelids and the ancestor to chordates.The gene encoding a hemoglobin (leghemoglobin) of higher plantsshows exactly the same splice junctions as in the globin genesof vertebrates except that the middle exon is split by an additionalintron which is believed to have been lost early in animal evolution.The occurrence of hemoglobins in diverse higher plants suggeststhat they might be present in all plants albeit at very lowconcentrations and perhaps serving an enzymatic function. Hemoglobin,broadly defined as a heme-containing protein capable of reversiblecombination with oxygen, also occurs in bacteria and fungi.The discovery of a bacterial hemoglobin that is 26% identicalwith lupin leghemoglobin indicates a procaryotic origin. Thepossibility that hemoglobin may have evolved from a cytochromeis suggested by the presence of hemoglobins in the yeast, Candida,and the bacterium, Alcaligenes, that contain both heme and flavin.They must therefore have evolved by the fusion of the genesfor two different proteins. However, the possible homology ofthe heme domains with plant or animal hemoglobins remains tobe determined. The lactic dehydrogenase of yeast, cytochromeb2, is also a soluble flavoheme protein with a heme domain thatis homologous with mammalian cytochrome b5. Globin may haveevolved in part from a member of the cytochrome b5 family, butif so, the event must have occurred so early that only a borderlineperhaps random correspondence of amino acid sequences remains     

Slow mitochondrial DNA sequence evolution in the Anthozoa (Cnidaria)

Mitochondrial genes have been used extensively in population genetic and phylogeographical analyses, in part due to a high rate of nucleotide substitution in animal mitochondrial DNA (mtDNA). Nucleotide sequences of anthozoan mitochondrial genes, however, are virtually invariant among conspecifics, even at third codon positions of protein-coding sequences. Hence, mtDNA markers are of limited use for population-level studies in these organisms. Mitochondrial gene sequence divergence among anthozoan species is also low relative to that exhibited in other animals, although higher level relationships can be resolved with these markers. Substitution rates in anthozoan nuclear genes are much higher than in mitochondrial genes, whereas nuclear genes in other metazoans usually evolve more slowly than, or similar to, mitochondrial genes. Although several mechanisms accounting for a slow rate of sequence evolution have been proposed, there is not yet a definitive explanation for this observation. Slow evolution and unique characteristics may be common in primitive metazoans, suggesting that patterns of mtDNA evolution in these organisms differ from that in other animal systems.     

Phylogenetic analysis of polymorphic DNA sequences at the Adh locus inDrosophila melanogaster and its sibling species

Summary Recent sequencing of over 2300 nucleotides containing the alcohol dehydrogenase (Adh) locus in each of 11Drosophila melanogaster lines makes it possible to estimate the approximate age of the electrophoretic fast-slow polymorphism. Our estimates, based on various possible patterns of evolution, range from 610,000 to 3,500,000 years, with 1,000,000 years as a reasonable point estimate. Furthermore, comparison of these sequences with those of the homologous region ofD. simulans andD. mauritiana allows us to infer the pattern of evolutionary change of theD. melanogaster sequences. The integrity of the Adh-f electrophoretic alleles as a single lineage is supported by both unweighted pair-group method (UPGMA) and parsimony analyeses. However, considerable divergence among the Adh-s lines seems to have preceded the origin of the Adh-f allele. Comparisons of the sequences ofD. melanogaster genes with those ofD. simulans andD. mauritiana genes suggest that the split between the latter two species occurred more recently than the divergence of some of the present-day Adh-s genes inD. melanogaster. The phylogenetic analyses of theD. melanogaster sequences show that the fastslow distinction is not perfect, and suggest that intragenic recombination or gene conversion occurred in the evolution of this locus. We extended conventional phylogenetic analyses by using a statistical technique for detecting and characterizing recombination events. We show that the pattern of differentiation of DNA sequences inD. melanogaster is roughly compatible with the neutral theory of molecular evolution.     

Independent gene evolution in the potato actin gene family demonstrated by phylogenetic procedures for resolving gene conversions and the phylogeny of angiosperm actin genes

Summary Nine different actin DNA sequences were isolated from the common potato,Solanum tuberosum, and the nucleotide sequence of five actin loci and of two allelic variants are presented. Unlike the wide variation in intron position among animal actin genes, the potato actin genes have three introns situated in the same positions as reported for all other angiosperm actin genes. Using a novel combination of analytical procedures (G-test and compatibility analysis), we could not find evidence of frequent large or small nonreciprocal exchanges of genetic material between the sequenced loci, although there were a few candidates. Resolution of such gene conversion events and the quantification of independence of gene evolution in multigene families is critical to the inference of phylogenetic relationships. Comparison with actin genes in other angiosperm species suggests that the actin multigene family can be divided into a number of subfamilies, evolved by descent rather than gene conversion, which are of possible functional origin, with one major subfamily diversification occurring before the divergence of monocots and dicots. The silent rate of nucleotide substitution was estimated to be similar to that suggested for a number of other plant nuclear genes, whereas the replacement rate was extremely slow, suggestive of selective constraints.     

Chromosomal arrangement of leghemoglobin genes in soybean.

A cluster of four different leghemoglobin (Lb) genes was isolated from AluI-HaeIII and EcoRI genomic libraries of soybean in a set of overlapping clones which together include 45 kilobases (kb) of contiguous DNA. These four genes, including a pseudogene, are present in the same orientation and are arranged in the order: 5'-Lba-Lbc1-Lb psi-Lbc3-3'. The intergenic regions average 2.5 kb. In addition to this main Lb locus, there are other Lb genes which do not appear to be contiguous to this locus. A sequence probably common to the 3' region of Lb loci was found flanking the Lbc3 gene. The 3' flanking region of the main Lb locus also contains a sequence that appears to be expressed more abundantly in root tissue. Another sequence which is primarily expressed in root and leaf is found 5' to two Lb loci. Overall, the main leghemoglobin locus is similar in structure to the mammalian globin gene loci.     

Molecular evolutionary analysis of theywvz/7B globin gene cluster of the insectChironomus thummi

We report the sequence of 8.1 kb of DNA containing the 3 end of one and seven other complete intronless globin genes from theywvz/7B locus of the dipteranChironomus thummi thummi. One of these (cttv) appears to be a pseudogene by virtue of a premature termination codon, whereas the others encode apparently functional globin polypeptides. Taken together with previously published data, theC. th. thummi ywvz/7B locus codes for at least 11 globins, five of which differ from one another by no more than two amino acids. In contrast, only nine globin genes are found in a comparable genomic clone isolated fromC. th. piger. As indicated by sequence alignment, this difference in copy number can be attributed to a loss of one gene (fusion of globin genes 7B8 and 7B10) in thepiger lines, coupled with a gain (globin gene 7139) in thethummi lineage. Comparisons between thethummi andpiger sequences showed thatywvz/713 intergenic regions have maintained a level of 91 % similarity since thethummi/piger divergence: most differences are simply due to single base substitutions or insertion/deletion events in either thethummi or thepiger DNA, but three instances of partially overlapping deletions were also detected. A phylogenetic analysis ofywvz/713 gene products was conducted, from which a plausible reconstruction of the evolutionary history of the locus was obtained. In addition, alignment of globin 7B amino acid sequences suggested that globin genes 7B2 and 7B3 (reported at the protein and cDNA level, respectively, but not contained on theC. th. thummi orC. th. piger genomic clones) are possibly chimeric genes. Given the trend toward expansion of theC. thummi globin gene family in general and of the globin 7B subfamily in particular, we propose that increased copy number of these genes has been positively selected as a mechanism to achieve a high Hb concentration in the larval hemolymph. Correspondence to: G. Bergtrom     

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.