Nucleotide sequence analysis of the delta beta-globin gene region in humans

The continuous DNA sequence of a 16.5-kilobase pair region encompassing the linked delta beta-globin gene cluster in humans is presented with a detailed restriction endonuclease map. There are 38 differences (0.5%) in comparison with published sequence data, corrected for errors in sequencing, resulting in polymorphic rates of 0.2% in exons and 0.76% in 5'-gene flanking regions. Fifteen changes result in the generation or elimination of restriction sites which may be useful in linkage disequilibrium studies. Two pairs of inverted Alu repeats, a pyrimidine-rich region 5' to delta, and (TG)n, (Pu/Py)n, and (ATTTT)n tracts 5' to beta are described. Dinucleotide frequencies and deviation from expected values approximated those found in total human genomic DNA. Regions of less than 50% A + T content were found associated with Alu sequences, a 150-base pair region immediately 5' to the beta gene, exon regions from both genes, and an area 3' to the beta gene. These regions also contained significantly lower than expected CpG levels compared to other regions, suggesting a possible relationship between DNA organizational patterns and functionally important regions. In addition, strand asymmetries in base composition in this region differ from those associated with the fetal globin genes.     

The globin gene family of the cephalochordate amphioxus: implications for chordate globin evolution

Background  

The lancelet amphioxus (Cephalochordata) is a close relative of vertebrates and thus may enhance our understanding of vertebrate gene and genome evolution. In this context, the globins are one of the best studied models for gene family evolution. Previous biochemical studies have demonstrated the presence of an intracellular globin in notochord tissue and myotome of amphioxus, but the corresponding gene has not yet been identified. Genomic resources of Branchiostoma floridae now facilitate the identification, experimental confirmation and molecular evolutionary analysis of its globin gene repertoire.     

The sequence of the gorilla fetal globin genes: evidence for multiple gene conversions in human evolution

Two fetal globin genes (G gamma and A gamma) from one chromosome of a lowland gorilla (Gorilla gorilla gorilla) have been sequenced and compared to three human loci (a G gamma-gene and two A gamma-alleles). A comparison of regions of local homology among these five sequences indicates that long after the duplication that produced the two nonallelic gamma-globin loci of catarrhine primates, about 35 million years (Myr) ago, at least one gene conversion event occurred between these loci. This conversion occurred not long before the ancestral divergence (about 6 Myr ago) of Homo and Gorilla. After this ancestral divergence, a minimum of three more gene conversion events occurred in the human lineage. Each human A gamma-allele shares specific sequence features with the gorilla A gamma-gene; one such distinctive allelic feature involves the simple repeated sequence in IVS 2. This suggests that early in the human lineage the A gamma-genes may have undergone a crossing-over event mediated by this simple repeated sequence. The DNA sequences from coding regions of both G gamma- and A gamma-loci, a comparison of 292 codons in the corresponding gorilla and human genes, show an unusually low evolutionary rate, with only two nonsilent differences and, surprisingly, not even one silent substitution. The two nonsynonymous substitutions observed predict a glycine at codon 73 and an arginine at codon 104 in the gorilla A gamma-sequence rather than aspartic acid and lysine, respectively, in human A gamma. Because only arginine has been found at position 104 in gamma-chains of Old World monkeys, it may represent the ancestral residue lost in gorilla and human G gamma-chains and in the human A gamma-chain. Possibly the arginine codon (AGG) was replaced by the lysine codon (AAG) in the G gamma-gene of a common ancestor of Homo and Gorilla and then was transferred to the A gamma-gene by subsequent conversions in the human lineage. DNA sequence conversions, similar to that attributed to the fetal gamma-globin genes, appear to be relatively frequent phenomena and, if widespread throughout the genome, may have profound evolutionary consequences.      

Nucleotide sequence analysis of human beta-globin gene by the quantification method: mutations in 3'-splice junction sequence and beta-thalassemia

The nucleotide sequence at the intron-exon junction in the human beta-globin gene was analyzed by the quantification method (categorical discriminant analysis) proposed previously. Using the sample score of a 16-nucleotide sequence at a 3'-splice junction, we studied to what extent such a sequence contains the 3'-splice signal. To examine the applicability of our method, we further studied several mutants of beta-thalassemia, where nucleotide changes exist at 3'-splice junction sequences of the first and second introns. Other mutants involve point mutations which generate new 3'-splice signals within the first intron. Experimental results on the abnormal splicing in those mutants could be explained in terms of the sample scores of 16-nucleotide sequences and their locations relative to the branch point.     

The evolution of a plant globin gene family

Summary We have analyzed the sequences of soybean leghemoglobin genes as an initial step toward understanding their mode of evolution. Alignment of the sequences of plant globin genes with those of animals reveals that (i) based on the proportion of nucleotide substitutions that have occurred at the first, second, and third codon positions, the time of divergence of plant and animal globin gene families appears to be extremely remote (between 900 million and 1.4 billion years ago, if one assumes constancy of evolutionary rate in both the plant and animal lineages) and (ii) in addition to the normal regulatory sequences on the 5 end, an approximately 30-base-pair sequence, specific to globin genes, that surrounds the cap site is conserved between the plant and animal globin genes. Comparison of the leghemoglobin sequences with one another shows that (i) the relative amount of sequence divergence in various coding and noncoding regions is roughly similar to that found for animal globin genes and (ii) as in animal globin genes, the positions of insertions and deletions in the intervening sequences often coincide with the locations of direct repeats. Thus, the mode of evolution of the plant globin genes appears to resemble, in many ways, that of their animal counterparts. We contrast the overall intergenic organization of the plant globin genes with that of animal genes, and discuss the possibility of the concerted evolution of the leghemoglobin genes.     

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.     

Nucleotide polymorphism in colicin E2 gene clusters: evidence for nonneutral evolution

To explore the molecular mechanisms behind the diversification of colicin gene clusters, we examined DNA sequence polymorphism for the colicin gene clusters of 14 colicin E2 (ColE2) plasmids obtained from natural isolates of Escherichia coli. Two types of ColE2 plasmids are revealed, with type II gene clusters generated by recombination between type I ColE2 and ColE7 gene clusters. The levels and patterns of DNA polymorphism are different between the two types. Type I polymorphism is distributed evenly along the gene cluster, while type II accumulates polymorphism at an elevated rate in the 5' end of the colicin gene. These differences may be explained by recombinational origins of type II gene clusters. The pattern of divergence between the ColE2 gene cluster and its close relative ColE9 is not correlated with the pattern of polymorphism within ColE2, suggesting that this gene cluster is not evolving in a neutral fashion. A statistical test confirms significant departures from the predictions of neutrality. These data lend further support to the hypothesis that colicin gene clusters may evolve under the influence of nonneutral forces.      

Nucleotide sequence of the BALB/c mouse beta-globin complex

The nucleotide sequence of 55,856 base-pairs containing all seven beta-globin homologous structures from chromosome 7 of the BALB/c mouse is reported. This sequence links together previously published sequences of the beta-globin genes, pseudogenes and repetitive elements. Using low stringency computer searches, we found no additional beta-globin homologous sequences, but did find many more long interspersed repetitive sequences (L1) than predicted by hybridization. L1 is a major component of the mouse beta-globin complex with at least 15 elements comprising about 22% of the reported sequence. Most open reading frames greater than 300 base-pairs in the cluster overlap with L1 repeats or globin genes. Polypurine, polypyrimidine and alternating purine/pyrimidine tracts are not evenly dispersed throughout the complex, but they do not appear to be excluded from or restricted to particular regions. Several regions of intergenic homology were detected in dot-plot comparisons of the mouse sequence with itself and with the human beta-globin sequence. The significance of these homologies is unclear, but these regions are candidates for further study in functional assays in erythroid cell lines or transgenic animals.     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

Nucleotide sequence of the soybean mitochondrial 18S rRNA gene: evidence for a slow rate of divergence in the plant mitochondrial genome

Summary The nucleotide sequence of the 18S rRNA gene from soybean mitochondria was determined and is presented here in comparison to the 18S rRNA genes from wheat and maize mitochondria. All three genes exhibit remarkable sequence similarity supporting the proposal that there is a slower rate of nucleotide divergence in plant mitochondrial DNA (mtDNA) as compared to the mtDNA of animals. A lower degree of sequence similarity is observed between the dicotyledenous plant soybean and either wheat (84%) or maize (85%) than between the two monocots (96%). A possible secondary structure for the soybean 18S rRNA is presented that is analogous to the proposed structure for the E. coli 16S rRNA.     

Conformational changes in the globin family during evolution

Summary Conformational restrictions imposed on the fixation of insertions and deletions by the three-dimensional structure of globins during evolution are analyzed. An evolutionary tree for the primary structure of 38 taxonomically distant globins was constructed by a computer method. Based on this tree, the calculated fixation frequency of point mutations was 50-fold higher than that established jointly for deletions and insertions, and the fixation frequency of deletions was more than three times that established for insertions. It was also found that deletions and insertions are predominantly fixed in the interhelical sections and at the ends of the -helices of the globin molecules. Conformational analysis of the packing of the -helices in the spatial structure of globins demonstrated that the fixation of deletions and insertions in the center of an -helix produces a marked distortion of their normal packing. The possible role of deletions and insertions in the evolution of protein families is discussed.     

Molecular cloning and sequence analysis of two rat major globin cDNAs

Two cDNA clones for globins of the adult Wistar rat were isolated from a reticulocyte cDNA library and the nucleotide sequences of the inserts were determined. One clone contained a cDNA insert consisting of 556 bp and the other contained one of 577 bp, both covering the entire coding sequences for rat globins. Comparisons of their predicted amino acid sequences with known sequences of rat globins revealed that these cDNAs coded for a rat major alpha- and a major beta-globin, I alpha and II beta, respectively. The cause of diversity of rat globins was discussed in terms of the nucleotide sequences of cDNAs and known amino acid sequences of globins.     

Displacement synthesis of globin complementary DNA: evidence for sequence amplification

We have examined a cDNA displacement synthesis procedure in which the extent of precursor incorporation and the unusual kinetics of displacement synthesis suggest a unique replicative form of DNA and the occurrence of multiple rounds of displacement synthesis, leading to amplification of mRNA sequences. Globin double-stranded DNA containing a hairpin loop was extended by the addition of a homopolymer to the 3' end. This was followed by displacement synthesis with the Klenow fragment of DNA polymerase I that was primed by an oligonucleotide hybridized to the homopolymer. Thus, the hairpin cDNA was copied to form an open duplex with an inverted repetition of globin sequences. These molecules can then serve as templates for additional synthesis which would be primed from oligomers bound the homopolymer. Globin cDNA sequences appear to be amplified 10-fold or more by this procedure. Globin cDNA obtained by displacement synthesis was similar in size to the original template. However, displaced molecules associate to the extent that they are not readily resolved by electrophoresis or sedimentation under nondenaturing conditions. Restriction endonuclease digests of 32P-labeled displaced strands gave fragment patterns similar to rabbit globin cDNA hairpin molecules. S1 nuclease studies demonstrated that displaced complexes and replication intermediates are partially single stranded, which might account for their aggregation properties.