The BC200 RNA Gene and Its Neural Expression Are Conserved in Anthropoidea (Primates)

The gene encoding BC200 RNA arose from a monomeric Alu element. Subsequently, the RNA had been recruited or exapted into a function of the nervous system. Here we confirm the presence of the BC200 gene in several primate species among the Anthropoidea. The period following the divergence of New World monkeys and Old World monkeys from their common ancestor is characterized by a significantly higher substitution rate in the examined 5′ flanking region than in the BC200 RNA coding region itself. Furthermore, the conservation of CpG dimers in the RNA coding region (200 bp) is drastically increased compared to the 5′ flanking region (∼400 bp) over all 12 species examined. Finally, the brain-specific expression pattern of BC200 RNA and its presence as a ribonucleoprotein particle (RNP) are conserved in Old World and New World monkeys. Our studies indicate that the gene encoding BC200 RNA was created at least 35–55 million years ago and its presence, mode of expression, and association with protein(s) as an RNP are under selective pressure. Received: 1 December 1997 / Accepted: 3 June 1998     

Mhc Allelic Diversity and Modern Human Origins

Thirty complete coding sequences of human major histocompatibility complex (Mhc) class II DRB alleles, spanning 237 codons, were analyzed for phylogenetic information using distance, parsimony, and likelihood approaches. Allelic genealogies derived from different parts of the coding sequence (exon 2, the 5′ and 3′ ends of exon 2, respectively, and exons 3–6) were compared. Contrary to prior assertions, a rigorous analysis of allelic genealogies in this gene family cannot be used to justify the claim that the lineage leading to modern humans contained on average at least 100,000 individuals. Phylogenetic inferences based upon the exon 2 region of the DRB loci are complicated by selection and recombination, so this part of the gene does not provide a complete and accurate view of allelic relationships. Attempts to reconstruct human history from genetic data must use realistic models which consider the complicating factors of nonequilibrium populations, recombination, and different patterns of selection. Received: 19 February 1997 / Accepted: 12 June 1997     

Duplication and Diversification of the Apolipoprotein CI (APOCI) Genomic Segment in Association with Retroelements

We have previously shown that several multicopy gene families within the major histocompatibility complex (MHC) arose from a process of segmental duplication. It has also been observed that retroelements play a role in generating diversity within these duplicated segments. The objective of this study was to compare the genomic organization of a gene duplication within another multicopy gene family outside the MHC. Using new continuous genomic sequence encompassing the APOE-CII gene cluster, we show that APOCI and its pseudogene, APOCI′, are contained within large duplicated segments which include sequences from the hepatic control region (HCR). Flanking Alu sequences are observed at both ends of the duplicated unit, suggesting a possible role in the integration of these segments. As observed previously within the MHC, the major differences between the segments are the insertion of sequences (approximately 200–1000 bp in length), consisting predominantly of Alu sequences. Ancestral retroelements also contribute to the generation of sequence diversity between the segments, especially within the 3′ poly(A) tract of Alu sequences. The exonic and regulatory sequences of the APOCI and HCR loci show limited sequence diversity, with exon 3 being an exception. Finally, the typing of pre- and postduplication Alus from both segments indicates an estimated time of duplication of approximately 37 million years ago (mya), some time prior to the separation of Old and New World monkeys. Received: 17 July 1999 / Accepted: 6 November 1999     

Phylogenetic Relationships Within the Genus Equus and the Evolution of α and θ Globin Genes

Sequences of the α1, α2 and θ globin genes from six equid species have been determined to investigate relationships within the genus Equus. Analyses using standard phylogenetic methods, or an approach designed to account for the effects of gene conversion between the α genes, gave broadly similar results and show that the horses diverged from the zebra/ass ancestor ∼2.4 million years ago and that the zebra and ass species arose in a rapid radiation ∼0.9 million years ago. These results from the α genes are corroborated by θ gene data and are in contrast to mitochondrial DNA studies of the phylogeny of this genus, which suggest a more gradual set of speciation events. Received: 22 April 1997 / Accepted: 20 July 1998     

Phylogenetic Position of Mammoth and Steller's Sea Cow Within Tethytheria Demonstrated by Mitochondrial DNA Sequences

Here we report DNA sequences from mitochondrial cytochrome b gene segments (1,005 base pairs per species) for the extinct woolly mammoth (Mammuthus primigenius) and Steller's sea cow (Hydrodamalis gigas) and the extant Asian elephant (Elephas maximus), the Western Indian manatee (Trichechus manatus), and the hyrax (Procavia capensis). These molecular data have allowed us to construct the phylogeny for the Tethytheria. Our molecular data resolve the trichotomy between the two species of living elephants and the mammoth and confirm that the mammoth was more closely related to the Asian elephant than to the African elephant. Our data also suggest that the sea cow–dugong divergence was likely as ancient as the dugong–manatee split, and it appears to have been much earlier (22 million years ago) than had been previously estimated (4–8 million years ago) by immunological comparison. Received: 8 August 1996 / Accepted: 30 September 1996     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

Molecular Clock Calibrations and Metazoan Divergence Dates

It has recently been argued that living metazoans diverged over 800 million years ago, based on evidence from 22 nuclear genes for such a deep divergence between vertebrates and arthropods (Gu 1998). Two ``internal' calibration points were used. However, only one fossil divergence date (the mammal–bird split) was directly used to calibrate the molecular clock. The second calibration point (the primate–rodent split) was based on molecular estimates that were ultimately also calibrated by the same mammal–bird split. However, the first tetrapods that can be assigned with confidence to either the mammal (synapsid) lineage or the bird (diapsid) lineage are approximately 288 million years old, while the first mammals that can be assigned with confidence to either the primate or the rodent lineages are 65 million years old, or 85 million years old if ferungulates are part of the primate lineage and zhelestids are accepted as ferungulate relatives. Recalibration of the protein data using these fossil dates indicates that metazoans diverged between 791 and 528 million years ago, a result broadly consistent with the palaeontological documentation of the ``Cambrian explosion.' The third, ``external' calibration point (the metazoan–fungal divergence) was similarly problematic, since it was based on a controversial molecular study (which in turn used fossil dates including the mammal–bird split); direct use of fossils for this calibration point gives the absurd dating of 455 million years for metazoan divergences. Similar calibration problems affect another recent study (Wang et al. 1999), which proposes divergences for metazoans of 1000 million years or more: recalibrations of their clock again yields much more recent dates, some consistent with a ``Cambrian explosion' scenario. Molecular clock studies have persuasively argued for the imperfection of the fossil record but have rarely acknowledged that their inferences are also directly based on this same record. Received: 26 January 1999 / Accepted: 14 April 1999     

Heterogeneity in Codon Usage in the Flatworm Schistosoma mansoni

Synonymous codon choices vary considerably among Schistosoma mansoni genes. Principal components analysis detects a single major trend among genes, which highly correlates with GC content in third codon positions and exons, but does not discriminate among putatively highly and lowly expressed genes. The effective number of codons used in each gene, and its distribution when plotted against GC₃, suggests that codon usage is shaped mainly by mutational biases. The GC content of exons, GC₃, 5′, 3′, and flanking (5′+ 3′+ introns) regions are all correlated among them, suggesting that variations in GC content may exist among different regions of the S. mansoni genome. We propose that this genome structure might be among the most important factors shaping codon usage in this species, although the action of selection on certain sequences cannot be excluded. Received: 10 March 1997 / Accepted: 27 June 1997     

Rapid Evolution of the Ribonuclease A Superfamily: Adaptive Expansion of Independent Gene Clusters in Rats and Mice

The two eosinophil ribonucleases, eosinophil-derived neurotoxin (EDN/RNase 2) and eosinophil cationic protein (ECP/RNase 3), are among the most rapidly evolving coding sequences known among primates. The eight mouse genes identified as orthologs of EDN and ECP form a highly divergent, species-limited cluster. We present here the rat ribonuclease cluster, a group of eight distinct ribonuclease A superfamily genes that are more closely related to one another than they are to their murine counterparts. The existence of independent gene clusters suggests that numerous duplications and diversification events have occurred at these loci recently, sometime after the divergence of these two rodent species (∼10–15 million years ago). Nonsynonymous substitutions per site (d _N) calculated for the 64 mouse/rat gene pairs indicate that these ribonucleases are incorporating nonsilent mutations at accelerated rates, and comparisons of nonsynonymous to synonymous substitution (d _N / d _S) suggest that diversity in the mouse ribonuclease cluster is promoted by positive (Darwinian) selection. Although the pressures promoting similar but clearly independent styles of rapid diversification among these primate and rodent genes remain uncertain, our recent findings regarding the function of human EDN suggest a role for these ribonucleases in antiviral host defense. Received: 8 April 1999 / Accepted: 22 June 1999     

The Evolution of MHC Diversity by Segmental Duplication and Transposition of Retroelements

Sequence analysis of a 237 kb genomic fragment from the central region of the MHC has revealed that the HLA-B and HLA-C genes are contained within duplicated segments peri-B (53 kb) and peri-C (48 kb), respectively, and separated by an intervening sequence (IF) of 30 kb. The peri-B and peri-C segments share at least 90% sequence homology except when interrupted by insertions/deletions including Alu, L1, an endogenous retrovirus, and pseudogenes. The sequences of peri-B, IF, and peri-C were searched for the presence of Alu elements to use as markers of evolution, chromosomal rearrangements, and polymorphism. Of 29 Alu elements, 14 were identified in peri-B, 11 in peri-C, and 4 in IF. The Alu elements in peri-B and peri-C clustered phylogenetically into two clades which were classified as ``preduplication' and ``postduplication' clades. Four Alu J elements that are shared by peri-B and peri-C and are flanked by homologous sequences in their paralogous locations, respectively, clustered into a ``preduplication' clade. By contrast, the majority of Alu elements, which are unique to either peri-B or peri-C, clustered into a postduplication clade together with the Alu consensus subfamily members ranging from platyrrhine-specific (Spqxcg) to catarrhine-specific Alu sequences (Y). The insertion of platyrrhine-specific Alu elements in postduplication locations of peri-B and peri-C implies that these two segments are the products of a duplication which occurred in primates prior to the divergence of the New World primate from the human lineage (35–44 mya). Examination of the paralogous Alu integration sites revealed that 9 of 14 postduplication Alu sequences have produced microsatellites of different length and sequence within the Alu 3′-poly A tail. The present analysis supports the hypothesis that HLA-B and HLA-C genes are products of an extended segmental duplication between 44 and 81 million years ago (mya), and that subsequent diversification of both genomic segments occurred because of the mobility and mutation of retroelements such as Alu repeats. Received: 21 May 1997 / Accepted: 9 July 1997     

Activation of the Cardiac Ryanodine Receptor by Sulfhydryl Oxidation is Modified by Mg2+ and ATP

The reactive disulfide 4,4′-dithiodipyridine (4,4′DTDP) was added to single cardiac ryanodine receptors (RyRs) in lipid bilayers. The activity of native RyRs, with cytoplasmic (cis) [Ca²⁺] of 10⁻⁷ m (in the absence of Mg²⁺ and ATP), increased within ∼1 min of addition of 1 mm 4,4′-DTDP, and then irreversibly ceased 5 to 6 min after the addition. Channels, inhibited by either 1 mm cis Mg²⁺ (10⁻⁷ m cis Ca²⁺) or by 10 mm cis Mg²⁺ (10⁻³ m cis Ca²⁺), or activated by 4 mm ATP (10⁻⁷ m cis Ca²⁺), also responded to 1 mm cis 4,4′-DTDP with activation and then loss of activity. P _o and mean open time (T _o ) of the maximally activated channels were lower in the presence of Mg²⁺ than in its absence, and the number of openings within the long time constant components of the open time distribution was reduced. In contrast to the reduced activation by 1 mm 4,4′-DTDP in channels inhibited by Mg²⁺, and the previously reported enhanced activation by 4,4′-DTDP in channels activated by Ca²⁺ or caffeine (Eager et al., 1997), the activation produced by 1 mm cis 4,4′-DTDP was the same in the presence and absence of ATP. These results suggest that there is a physical interaction between the ATP binding domain of the cardiac RyR and the SH groups whose oxidation leads to channel activation. Received: 8 September 1997/Revised: 20 January 1998     

Phylogenetic Affinity of Mitochondria of Euglena gracilis and Kinetoplastids Using Cytochrome Oxidase I and hsp60

The mitochondrial DNA-encoded cytochrome oxidase subunit I (COI) gene and the nuclear DNA-encoded hsp60 gene from the euglenoid protozoan Euglena gracilis were cloned and sequenced. The COI sequence represents the first example of a mitochondrial genome-encoded gene from this organism. This gene contains seven TGG tryptophan codons and no TGA tryptophan codons, suggesting the use of the universal genetic code. This differs from the situation in the mitochondrion of the related kinetoplastid protozoa, in which TGA codes for tryptophan. In addition, a complete absence of CGN triplets may imply the lack of the corresponding tRNA species. COI cDNAs from E. gracilis possess short 5′ and 3′ untranslated transcribed sequences and lack a 3′ poly[A] tail. The COI gene does not require uridine insertion/deletion RNA editing, as occurs in kinetoplastid mitochondria, to be functional, and no short guide RNA-like molecules could be visualized by labeling total mitochondrial RNA with [α-³²P]GTP and guanylyl transferase. In spite of the differences in codon usage and the 3′ end structures of mRNAs, phylogenetic analysis using the COI and hsp60 protein sequences suggests a monophyletic relationship between the mitochondrial genomes of E. gracilis and of the kinetoplastids, which is consistent with the phylogenetic relationship of these groups previously obtained using nuclear ribosomal RNA sequences. Received: 5 March 1996 / Accepted: 31 July 1996     

Modulation of the Calmodulin-induced Inhibition of Sarcoplasmic Reticulum Calcium Release Channel (Ryanodine Receptor) by Sulfhydryl Oxidation in Single Channel Current Recordings and [3H]Ryanodine Binding

The modulation of the calmodulin-induced inhibition of the calcium release channel (ryanodine receptor) by two sulfhydryl oxidizing compounds, 4-(chloromercuri)phenyl–sulfonic acid (4-CMPS) and 4,4′-dithiodipyridine (4,4′-DTDP) was determined by single channel current recordings with the purified and reconstituted calcium release channel from rabbit skeletal muscle sarcoplasmic reticulum (HSR) and [³H]ryanodine binding to HSR vesicles. 0.1 μm CaM reduced the open probability (P _o ) of the calcium release channel at maximally activating calcium concentrations (50–100 μm) from 0.502 ± 0.02 to 0.137 ± 0.022 (n= 28), with no effect on unitary conductance. 4-CMPS (10–40 μm) and 4,4′-DTDP (0.1–0.3 mm) induced a concentration dependent increase in P _o (> 0.9) and caused the appearance of longer open states. CaM shifted the activation of the calcium release channel by 4-CMPS or 4,4′-DTDP to higher concentrations in single channel recordings and [³H]ryanodine binding. 40 μm 4-CMPS induced a near maximal (P _o > 0.9) and 0.3 mm 4,4′-DTDP a submaximal (P _o = 0.74) channel opening in the presence of CaM, which was reversed by the specific sulfhydryl reducing agent DTT. Neither 4-CMPS nor 4,4′-DTDP affected Ca-[¹²⁵I]calmodulin binding to HSR. 1 mm MgCl₂ reduced P _o from 0.53 to 0.075 and 20–40 μm 4-CMPS induced a near maximal channel activation (P _o > 0.9). These results demonstrate that the inhibitory effect of CaM or magnesium in a physiological concentration is diminished or abolished at high concentrations of 4-CMPS or 4,4′-DTDP through oxidation of activating sulfhydryls on cysteine residues of the calcium release channel. Received: 22 July 1999/Revised: 15 November 1999     

Erratic Evolution of Glycerol-3-Phosphate Dehydrogenase in Drosophila,Chymomyza, and Ceratitis

We have studied the evolution of Gpdh in 18 fruitfly species by sequencing 1,077 nucleotides per species on average. The region sequenced includes four exons coding for 277 amino acids and three variable-length introns. Phylogenies derived by a variety of methods confirm that the nominal genus Zaprionus belongs within the genus Drosophila, whereas Scaptodrosophila and Chymomyza are outside. The rate of GPDH evolution is erratic. The rate of amino acid replacements in a lineage appears to be 1.0 × 10⁻¹⁰/site/year when Drosophila species are considered (diverged up to 55 million years ago), but becomes 2.3 × 10⁻¹⁰ when they are compared to Chymomyza species (divergence around 60 My ago), and 4.6 × 10⁻¹⁰ when species of those two genera are compared with the medfly Ceratitis capitata (divergence around 100 My ago). In order to account for these observations, the rate of amino acid replacement must have been 15 or more times greater in some lineages and at some times than in others. At the nucleotide level, however, Gpdh evolves in a fairly clockwise fashion. Received: 13 June 1996 / Accepted: 16 August 1996     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

The complete mitochondrial DNA (mtDNA) of the donkey and mtDNA comparisons among four closely related mammalian species-pairs

The nucleotide sequence of the complete mitochondrial genome of the donkey, Equus asinus, was determined. The length of the molecule is 16,670 bp. The length, however, is not absolute due to pronounced heteroplasmy caused by variable numbers of two types of repetitive motifs in the control region. The sequence of the repeats is (a) 5′-CACACCCA and (b) 5′-TGCGCGCA, respectively. The order of (a) and (b) can be expressed as {n[2(a)+(b)]+m(a)}. In 32 different clones analyzed the number of n and m ranged from 0 to 9 and 1 to 7. The two rRNA genes, the 13 peptide-coding genes, and the 22 tRNA genes of the donkey and the horse, Equus caballus, were compared in detail. Total nucleotide difference outside the control region was 6.9%. Nucleotide difference between peptide-coding genes ranged from 6.4% to 9.4% with a mean of 8.0%. In the inferred protein sequences of the 13 peptide-coding genes the amino acid difference was 0.2–8.8%, and the mean for the 13 concatenated amino acid sequences was 1.9%. In the 22 tRNA genes, the mean difference was 3.5%, and that in the two rRNA genes was 4.1%. The mtDNA differences between the donkey and the horse suggest that the evolutionary separation of the two species occurred ≈9 million years ago. Analyses of differences among the mtDNAs of three other species-pairs, harbor seal/grey seal, fin whale/blue whale, and Homo/common chimpanzee, showed that the relative evolutionary rate of individual peptide-coding genes varies among different species-pairs and modes of comparison. The findings show that the superimposition of sequence data of one lineage for resolving and dating evolutionary divergences of other lineages should be performed with caution unless based on comprehensive data. Received: 15 October 1995 / Accepted: 15 April 1996     

Partial Sequence of a Sponge Mitochondrial Genome Reveals Sequence Similarity to Cnidaria in Cytochrome Oxidase Subunit II and the Large Ribosomal RNA Subunit

A 2550-bp portion of the mitochondrial genome of a Demosponge, genus Tetilla, was amplified from whole genomic DNA extract and sequenced. The sequence was found to code for the 3′ end of the 16S rRNA gene, cytochrome c oxidase subunit II, a lysine tRNA, ATPase subunit 8, and a 5′ portion of ATPase subunit 6. The Porifera cluster distinctly within the eumetazoan radiation, as a sister group to the Cnidaria. Also, the mitochondrial genetic code of this sponge is likely identical to that found in the Cnidaria. Both the full COII DNA and protein sequences and a portion of the 16S rRNA gene were found to possess a striking similarity to published Cnidarian mtDNA sequences, allying the Porifera more closely to the Cnidaria than to any other metazoan phylum. The gene arrangement, COII—tRNA^Lys—ATP8—ATP6, is observed in many Eumetazoan phyla and is apparently ancestral in the metazoa. Received: 24 November 1997 / Accepted: 14 September 1998     

Phylogenetic Analysis of Vertebrate Lactate Dehydrogenase (LDH) Multigene Families

In this paper we analyzed 49 lactate dehydrogenase (LDH) sequences, mostly from vertebrates. The amino acid sequence differences were found to be larger for a human–killifish pair than a human–lamprey pair. This indicates that some protein sequence convergence may occur and reduce the sequence differences in distantly related species. We also examined transitions and transversions separately for several species pairs and found that the transitions tend to be saturated in the distantly related species pair, while transversions are increasing. We conclude that transversions maintain a conservative rate through the evolutionary time. Kimura's two-parameter model for multiple-hit correction on transversions only was used to derive a distance measure and then construct a neighbor-joining (NJ) tree. Three findings were revealed from the NJ tree: (i) the branching order of the tree is consistent with the common branch pattern of major vertebrates; (ii) Ldh-A and Ldh-B genes were duplicated near the origin of vertebrates; and (iii) Ldh-C and Ldh-A in mammals were produced by an independent gene duplication in early mammalian history. Furthermore, a relative rate test showed that mammalian Ldh-C evolved more rapidly than mammalian Ldh-A. Under a two-rate model, this duplication event was calibrated to be approximately 247 million years ago (mya), dating back to the Triassic period. Other gene duplication events were also discovered in Xenopus, the first duplication occurring approximately 60–70 mya in both Ldh-A and Ldh-B, followed by another recent gene duplication event, approximately 20 mya, in Ldh-B. Received: 5 October 2001 / Accepted: 24 October 2001     

Evolutionary History of B1 Retroposons in the Genus Mus

Short interspersed DNA elements (SINEs) amplify by retroposition either by (i) successive waves of amplification from one or a few evolving master genes or by (ii) the generation of new master genes that coexist with their progenitors. Individual, highly conserved, elements of the B1 SINE family were identified from the GenBank nucleotide database using various B1 subfamily consensus query sequences to determine their integration times into the mouse genome. A comparison of orthologous loci in various species of the genus Mus demonstrated that four subfamilies of B1 elements have been amplifying within the last 1–3 million years. Therefore, B1 sequences are generated by coexisting source genes. Additionally, three B1 subfamilies have been concurrently propagated during subspecies divergence and strain formation in Mus, indicating very recent activity of this retroposon family. The patterns of intra- and interspecies variations of orthologous loci demonstrate the usefulness of B1 integrations as a phylogenetic tool. A single inconsistency in the phylogenetic trends was depicted by the presence of a B1 insert in an orthologous locus exclusively in M. musculus and M. pahari. However, DNA sequence analysis revealed that these were independent integrations at the same genomic site. One highly conserved B1 element that integrated at least 4–6 million years ago suggests the possibility of occasional function for B1 integrations. Received: 25 February 2000 / Accepted: 5 June 2000     

Organization, Structure, and Evolution of the Nonadult Rat β-Globin Gene Cluster

The β-globin gene cluster of Wistar rat was extensively cloned and the embryonic genes were mapped and sequenced. Four overlapping λ Dash recombinant clones cover about 31 kb and contain four nonadult β-globin genes, 5′–ε1–γ1–γ2–ψγ3–3′. The ε1 and γ2 are active genes, since their protein products were detected in the fetal stage of the rat (Iwahara et al., J Biochem 119:360–366, 1996). The γ1 locus might be a pseudogene, since the ATA box in the promoter region is mutated to GTA; however, no other defect is observed. The ψγ3 locus is a truncated pseudogene because a 19-base deletion, which causes a shift of the reading frame, is observed between the second nucleotide of the putative codon 68 and codon 76. A sequence comparison suggests that the ψγ3 might be produced by a gene conversion event of the proto-γ-globin gene set. Possible histories of the evolution of rat nonadult β-globin genes are discussed. Received: 6 August 1998 / Accepted: 12 February 1999