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     

Using Alu J Elements as Molecular Clocks to Trace the Evolutionary Relationships Between Duplicated HLA Class I Genomic Segments

The class I region of the major histocompatibility complex contains two subgenomic blocks (250–350 kb each), known as the alpha and beta blocks. These blocks contain members of multicopy gene families including HLA class I, HERV-16 (previously called P5 sequences), and PERB11 (MIC). We have previously shown that each block consists of imperfect duplicated segments (duplicons) containing linked members of different gene families, retroelements and transposons that have coevolved as part of two separate evolutionary events. Another region provisionally designated here as the kappa block is located between the alpha and the beta blocks and contains HLA-E, -30, and -92, HERV-16 (P5.3), and PERB11.3 (MICC) within about 250 kb of sequence. Using Alu elements to trace the evolutionary relationships between different class I duplicons, we have found that (a) the kappa block contains paralogous (duplicated) Alu J sequences and other retroelement patterns more in common with the beta than the alpha block; (b) the retroelement pattern associated with the HLA-E duplicon is different from all other HLA class I duplicons, indicating a more complex evolution; (c) the HLA-92 duplicon, although substantially shorter, is closely related in sequence to the HLA-B and -C duplicons; (d) two of the six paralogous Alu J elements within the HLA-B and -C duplicons are associated with the HLA-X duplicon, confirming their evolutionary relationships within the beta block; and (e) the paralogous Alu J elements within the alpha block are distinctly different from those identified within the beta and kappa blocks. The sequence conservation and location of duplicated (paralogous) Alu J elements in the MHC class I region show that the beta and kappa blocks have evolved separately from the alpha block beginning at a time before or during the evolution of Alu J elements in primates. Received: 22 September 1999 / Accepted: 24 January 2000     

Genomic Organization Around the Centromeric End of the HLA Class I Region: Large-Scale Sequence Analysis

We previously sequenced two regions around the centromeric end of HLA class I and the boundary between class I and class III. In this paper we analyze the two regions of about 385 kb and confirm, giving a new line of evidence, that the following two pairs of the genomic segments were duplicated in evolution: (i) a 43-kb genomic segment including the HLA-B gene showing the highest polymorphism among the classical HLA class I loci (class Ia) and a 40-kb segment including the HLA-C locus showing the lowest polymorphism and (ii) a 52-kb segment including the MIC (MHC class I chain related gene) B and a 35-kb segment including MICA. We also found that repetitive elements such as SINEs, LINEs, and LTRs occupy as much as 47% of nucleotides in this 385-kb region. This unusually high content of repetitive elements indicates that repeat-mediated rearrangements have frequently occurred in the evolutionary history of the HLA class Ia region. Analysis of LINE compositions within the two pairs of duplicated segments revealed that (i) LINEs in these regions had been dispersed prior to both the duplication of the HLA-B and -C loci and the duplication of the MICB and MICA loci, and (ii) the divergence of the HLA-B and -C loci occurred prior to the duplication of the MICA and MICB loci. To find novel genes responsible for HLA class I-associated or other diseases, we performed computer analysis applying GenScan and GRAIL to GenBank's dbEST. As a result, at least five as yet uncharacterized genes were newly mapped on the HLA class I centromeric region studied. These novel genes should be analyzed further to determine their relationships to diseases associated with this region. Received: 16 June 1998 / Accepted: 18 August 1998     

Comparison Between Two Human Endogenous Retrovirus (HERV)-Rich Regions Within the Major Histocompatibility Complex

Sixteen human endogenous retrovirus (HERV) sequences were detected within 656 kb of genomic sequence obtained from the alpha- and beta-block of the class I region of the major histocompatibility complex (MHC). The HERVs were identified and characterized as family members of HERV-16 (11 copies), HERV-L (1 copy), HERV-I (2 copies), HERV-K91 (1 copy), and HARLEQUIN (1 copy) by sequence comparison using CENSOR or Repeat Masker, BLAST searches, and dot plots. The 11 copies of HERV-16 arose as products of duplication of genomic segments containing HLA class I (HLAcI) and PERB11 (MIC) genes inter alia, whereas the other five HERVs arose after duplication probably as a consequence of single insertion events or translocations. HERV-L and HERV-I are located between the duplicated genes PERB11.2 (MICB) and PERB11.1 (MICA), and HLA-B and HLA-C, respectively, whereas HERV-K91 and HARLEQUIN are located telomeric of HLA-C. A highly fragmented copy of HERV-I was also found telomeric of PERB11.4. Structural analysis of open reading frames (ORFs) revealed the absence of intact coding sequence within the putative gag, pol, and env gene regions of all the HERVs with the exception of HERV-K91, which had two large ORFs within the region of the putative protease and pol genes. In addition, the 5′-LTR of HERV-L contained a 2.5-kb element that was AT-rich and large ORFs with putative amino acid sequences rich in tyrosines and isoleucines. HERV-I, HARLEQUIN, and at least four copies of HERV-16 appear to have been receptors for the insertion of other retrotransposons including Alu elements and fragments of L1 and THE1. Examination of flanking sequences suggests that HERV-I and HERV-L had occurred by insertion into ancient L1 fragments. This study has revealed that the alpha- and beta-block region within the MHC is rich in HERV sequences occurring at a much higher ratio (10 to 1) than normally observed in the human genome. These HERV sequences will therefore enhance further studies on disease associations and differences between human haplotypes and primates and their role in the evolution of class I genes in the MHC. Received: 17 September 1998 / Accepted: 8 January 1999     

Coevolution of PERB11 (MIC) and HLA Class I Genes with HERV-16 and Retroelements by Extended Genomic Duplication

The recent availability of genomic sequence information for the class I region of the MHC has provided an opportunity to examine the genomic organization of HLA class I (HLAcI) and PERB11/MIC genes with a view to explaining their evolution from the perspective of extended genomic duplications rather than by simple gene duplications and/or gene conversion events. Analysis of genomic sequence from two regions of the MHC (the alpha- and beta-blocks) revealed that at least 6 PERB11 and 14 HLAcI genes, pseudogenes, and gene fragments are contained within extended duplicated segments. Each segment was searched for the presence of shared (paralogous) retroelements by RepeatMasker in order to use them as markers of evolution, genetic rearrangements, and evidence of segmental duplications. Shared Alu elements and other retroelements allowed the duplicated segments to be classified into five distinct groups (A to E) that could be further distilled down to an ancient preduplication segment containing a HLA and PERB11 gene, an endogenous retrovirus (HERV-16), and distinctive retroelements. The breakpoints within and between the different HLAcI segments were found mainly within the PERB11 and HLA genes, HERV-16, and other retroelements, suggesting that the latter have played a major role in duplication and indel events leading to the present organization of PERB11 and HLAcI genes. On the basis of the features contained within the segments, a coevolutionary model premised on tandem duplication of single and multipartite genomic segments is proposed. The model is used to explain the origins and genomic organization of retroelements, HERV-16, DNA transposons, PERB11, and HLAcI genes as distinct segmental combinations within the alpha- and beta-blocks of the human MHC. Received: 5 December 1998 / Accepted: 27 January 1999     

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     

Monophyletic Origin of Alu Elements in Primates

To get insight into the early evolution of the primate Alu elements, we characterized sequences of these repeats from the Malagasy prosimians, lemurs (Lemuridae) and sifakas (Indriidae), as well as from galagos (Lorisidae). These sequences were compared with the oldest Alu species known from the human genome: dimeric Alu J and S and free Alu monomers. Our analysis indicates that about 60 Myr ago, before the prosimian divergence, free left and right monomers formed an Alu heterodimer connected by a 19-nucleotide-long A-rich linker. The resulting elements successfully propagated in diverging primate lineages until about ∼20 Myr ago, conserving similar sequence features and essentially the same Alu RNA secondary structure. We suggest that until that time the same ``retropositional niche', molecular machinery making possible the proliferation by retroposition, constrained the evolution of Alu elements in extant primate species. These constraints became subsequently relaxed. In the Malagasy prosimians the dimeric Alu continued to amplify after acquiring a 34- to 36-nucleotide extension of their linker segment, whereas in the galago genome the ``retropositional niche' was occupied by novel short elements. Received: 1 December 1997 / Accepted: 30 January 1998     

Coevolution of HLA-B and PERB11.1 (MICA): Significance of Independent Triplet Expansion Within the Transmembrane Region of PERB11.1 (MICA)

Several highly polymorphic sequences are present in the beta block of the MHC, especially HLA-B, HLA-C, PERB11.1 (MICA), and PERB11.2 (MICB). It is now apparent that the polymorphism of PERB11.1 is of the same order as that of HLA-A, -B, and -C and it has been suggested that PERB11 could explain some of the disease associations previously attributed to HLA-B. Phylogenetic analysis of PERB11 α-domain sequences demonstrates relationships with HLA-B cross-reactive serogroups. In contrast, the transmembrane polymorphisms do not appear to be associated with either PERB11 or HLA-B. These data indicate that PERB11 and HLA-B have evolved in concert from their common ancestors and that the transmembrane polymorphisms have arisen independently and more recently. MHC disease associations will need to be reviewed in the light of mechanisms such as receptor binding and signaling. Received: 22 September 1999 / Accepted: 11 January 2000 (Revised)     

Molecular Evolution of Nitrogen Fixation: The Evolutionary History of the nifD, nifK, nifE, and nifN Genes

The pairs of nitrogen fixation genes nifDK and nifEN encode for the α and β subunits of nitrogenase and for the two subunits of the NifNE protein complex, involved in the biosynthesis of the FeMo cofactor, respectively. Comparative analysis of the amino acid sequences of the four NifD, NifK, NifE, and NifN in several archaeal and bacterial diazotrophs showed extensive sequence similarity between them, suggesting that their encoding genes constitute a novel paralogous gene family. We propose a two-step model to reconstruct the possible evolutionary history of the four genes. Accordingly, an ancestor gene gave rise, by an in-tandem paralogous duplication event followed by divergence, to an ancestral bicistronic operon; the latter, in turn, underwent a paralogous operon duplication event followed by evolutionary divergence leading to the ancestors of the present-day nifDK and nifEN operons. Both these paralogous duplication events very likely predated the appearance of the last universal common ancestor. The possible role of the ancestral gene and operon in nitrogen fixation is also discussed. Received: 21 June 1999 / Accepted: 1 March 2000     

Genomic and Phylogenetic Analysis of the Human CD1 and HLA Class I Multicopy Genes

The human CD1 proteins belong to a lipid-glycolipid antigen-presenting gene family and are related in structure and function to the MHC class I molecules. Previous mapping and DNA hybridization studies have shown that five linked genes located within a cluster on human chromosome 1q22-23 encode the CD1 protein family. We have analyzed the complete genomic sequence of the human CD1 gene cluster and found that the five active genes are distributed over 175,600 nucleotides and separated by four expanded intervening genomic regions (IGRs) ranging in length between 20 and 68 kb. The IGRs are composed mostly of retroelements including five full-length L1 PA sequences and various pseudogenes. Some L1 sequences have acted as receptors for other subtypes or families of retroelements. Alu molecular clocks that have evolved during primate history are found distributed within the HLA class I duplicated segments (duplicons) but not within the duplicons of CD1. Phylogeny of the alpha3 domain of the class I-like superfamily of proteins shows that the CD1 cluster is well separated from HLA class I by a number of superfamily members including MIC (PERB11), HFE, Zn-alpha2-GP, FcRn, and MR1. Phylogenetically, the human CD1 sequences are interspersed by CD1 sequences from other mammalian species, whereas the human HLA class I sequences cluster together and are separated from the other mammalian sequences. Genomic and phylogenetic analyses support the view that the human CD1 gene copies were duplicated prior to the evolution of primates and the bulk of the HLA class I genes found in humans. In contrast to the HLA class I genomic structure, the human CD1 duplicons are smaller in size, they lack Alu clocks, and they are interrupted by IGRs at least 4 to 14 times longer than the CD1 genes themselves. The IGRs seem to have been created as "buffer zones" to protect the CD1 genes from disruption by transposable elements.     

Rapid evolution of a young L1 (LINE-1) clade in recently speciated rattus taxa

L1 elements are retrotransposons that have been replicating and evolving in mammalian genomes since before the mammalian radiation. Rattus norvegicus shares the young L1_mlvi2 clade only with its sister taxon, Rattus cf moluccarius. Here we compared the L1_mlvi2 clade in these recently diverged species and found that it evolved rapidly into closely related but distinct clades: the L1_mlvi2-rm clade (or subfamily), characterized here from R. cf moluccarius, and the L1_mlvi2-rn clade, originally described in R. norvegicus. In addition to other differences, these clades are distinguished by a cluster of amino acid replacement substitutions in ORF I. Both rat species contain the L1_mlvi2-rm clade, but the L1_mlvi2-rn clade is restricted to R. norvegicus. Therefore, the L1_mlvi2-rm clade arose prior to the divergence of R. norvegicus and R. cf moluccarius, and the L1_mlvi2-rn clade amplified after their divergence. The total number of L1_mlvi2-rm elements in R. cf moluccarius is about the same as the sum of the L1_mlvi2-rm and L1_mlvi2-rn elements in R. norvegicus. The possibility that L1 amplification is in some way limited so that the two clades compete for replicative supremacy as well as the implications of the other distinguishing characteristic of the L1_mlvi2-rn and L1_mlvi2-rm clades are discussed. Received: 6 November 1996 / Accepted: 12 April 1997     

Six Opsins from the Butterfly Papilio glaucus: Molecular Phylogenetic Evidence for Paralogous Origins of Red-Sensitive Visual Pigments in Insects

It has been hypothesized that the UV-, blue-, and green-sensitive visual pigments of insects were present in the common ancestor of crustaceans and insects, whereas red-sensitive visual pigments evolved later as a result of convergent evolution. This hypothesis is examined with respect to the placement of six opsins from the swallowtail butterfly Papilio glaucus (PglRh1–6) in relationship to 46 other insect, crustacean, and chelicerate opsin sequences. All basal relationships established with maximum parsimony analysis except two are present in the distance and maximum likelihood analyses. In all analyses, the six P. glaucus opsins fall into three well-supported clades, comprised, respectively, of ultraviolet (UV), blue, and long-wavelength (LW) pigments, which appear to predate the radiation of the insects. Lepidopteran green- and red-sensitive visual pigments form a monophyletic clade, which lends support to the hypothesis from comparative physiological studies that red-sensitive visual pigments in insects have paralogous origins. Polymorphic amino acid sites (180, 197, 277, 285, 308), which are essential for generating the spectral diversity among the vertebrate red- and green-sensitive pigments are notably invariant in the Papilio red- and green-sensitive pigments. Other major tuning sites must be sought to explain the spectral diversification among these and other insect visual pigments. Received: 6 December 1999 / Accepted: 3 April 2000     

Diversity, Distribution, and Ancient Taxonomic Relationships Within the TIR and Non-TIR NBS-LRR Resistance Gene Subfamilies

Phylogenetic relationships among the NBS-LRR (nucleotide binding site–leucine-rich repeat) resistance gene homologues (RGHs) from 30 genera and nine families were evaluated relative to phylogenies for these taxa. More than 800 NBS-LRR RGHs were analyzed, primarily from Fabaceae, Brassicaceae, Poaceae, and Solanaceae species, but also from representatives of other angiosperm and gymnosperm families. Parsimony, maximum likelihood, and distance methods were used to classify these RGHs relative to previously observed gene subfamilies as well as within more closely related sequence clades. Grouping sequences using a distance cutoff of 250 PAM units (point accepted mutations per 100 residues) identified at least five ancient sequence clades with representatives from several plant families: the previously observed TIR gene subfamily and a minimum of four deep splits within the non-TIR gene subfamily. The deep splits in the non-TIR subfamily are also reflected in comparisons of amino acid substitution rates in various species and in ratios of nonsynonymous-to-synonymous nucleotide substitution rates (K _A/K _S values) in Arabidopsis thaliana. Lower K _A/K _S values in the TIR than the non-TIR sequences suggest greater functional constraints in the TIR subfamily. At least three of the five identified ancient clades appear to predate the angiosperm–gymnosperm radiation. Monocot sequences are absent from the TIR subfamily, as observed in previous studies. In both subfamilies, clades with sequences separated by approximately 150 PAM units are family but not genus specific, providing a rough measure of minimum dates for the first diversification event within these clades. Within any one clade, particular taxa may be dramatically over- or underrepresented, suggesting preferential expansions or losses of certain RGH types within particular taxa and suggesting that no one species will provide models for all major sequence types in other taxa. Received: 13 June 2001 / Accepted: 22 October 2001     

The Archaea Monophyly Issue: A Phylogeny of Translational Elongation Factor G(2) Sequences Inferred from an Optimized Selection of Alignment Positions

A global alignment of EF-G(2) sequences was corrected by reference to protein structure. The selection of characters eligible for construction of phylogenetic trees was optimized by searching for regions arising from the artifactual matching of sequence segments unique to different phylogenetic domains. The spurious matchings were identified by comparing all sections of the global alignment with a comprehensive inventory of significant binary alignments obtained by BLAST probing of the DNA and protein databases with representative EF-G(2) sequences. In three discrete alignment blocks (one in domain II and two in domain IV), the alignment of the bacterial sequences with those of Archaea–Eucarya was not retrieved by database probing with EF-G(2) sequences, and no EF-G homologue of the EF-2 sequence segments was detected by using partial EF-G(2) sequences as probes in BLAST/FASTA searches. The two domain IV regions (one of which comprises the ADP-ribosylatable site of EF-2) are almost certainly due to the artifactual alignment of insertion segments that are unique to Bacteria and to Archaea–Eucarya. Phylogenetic trees have been constructed from the global alignment after deselecting positions encompassing the unretrieved, spuriously aligned regions, as well as positions arising from misalignment of the G′ and G″ subdomain insertion segments flanking the ``fifth' consensus motif of the G domain (?varsson, 1995). The results show inconsistencies between trees inferred by alternative methods and alternative (DNA and protein) data sets with regard to Archaea being a monophyletic or paraphyletic grouping. Both maximum-likelihood and maximum-parsimony methods do not allow discrimination (by log-likelihood difference and difference in number of inferred substitutions) between the conflicting (monophyletic vs. paraphyletic Archaea) topologies. No specific EF-2 insertions (or terminal accretions) supporting a crenarchaeal–eucaryal clade are detectable in the new EF-G(2) sequence alignment.     

Intraspecies Analysis: Comparison of ITS Sequence Data and Gene Intron Sequence Data with Breeding Data for a Worldwide Collection of Gonium pectorale

The morphologically uniform species Gonium pectorale is a colonial green flagellate of worldwide distribution. The affinities of 25 isolates from 18 sites on five continents were assessed by both DNA sequence comparisons and sexual compatibility. Complete sequences were obtained (i) for the internal transcribed spacer ITS-1 and ITS-2 regions of ribosomal DNA and (ii) for each of three single-copy spliceosomal introns, two in a small G protein and one in the actin gene. ITS sequences appeared to homogenize sufficiently rapidly to behave as a single copy gene. Intron sequence differences between isolates in this species reached nucleotide substitution saturation, while ITS sequences did not. Parsimony and evolutionary distance analysis of the two types of DNA data gave essentially the same tree conformation. By all these criteria, the group of G. pectorale isolates fell into two main clades, A and B. Clade A, with isolates from four continents, was comprised of four subclades of quite closely related isolates, plus one strain of ambiguous affinity. Clade B was comprised of two subclades represented by South African and South American isolates, respectively; thus, only subclades of clade B showed geographical localization. With respect to mating, all isolates except one homothallic strain and one apparently sterile strain fell into either one or the other of two mating types. Pairings in all possible combinations revealed that isolates from the same site formed abundant zygotes, which germinated to produce new, sexually active organisms. Zygotes were also formed in many pairings of other combinations, including crosses of clade A with clade B organisms, but none of the latter produced viable germlings. The ability to mate and produce viable progeny that were themselves capable of sexual reproduction was restricted to members of subclades established on the basis of DNA sequence similarities. Thus, the grades of difference in both nuclear intron sequences and rDNA ITS sequences paralleled those observed in the sexual analysis. Received: 9 March 1998 / Accepted: 1 June 1998     

Gene Conversion Drives Within Genic Sequences: Concerted Evolution of Ribosomal RNA Genes in Bacteria and Archaea

Multiple copies of a given ribosomal RNA gene family undergo concerted evolution such that sequences of all gene copies are virtually identical within a species although they diverge normally between species. In eukaryotes, gene conversion and unequal crossing over are the proposed mechanisms for concerted evolution of tandemly repeated sequences, whereas dispersed genes are homogenized by gene conversion. However, the homogenization mechanisms for multiple-copy, normally dispersed, prokaryotic rRNA genes are not well understood. Here we compared the sequences of multiple paralogous rRNA genes within a genome in 12 prokaryotic organisms that have multiple copies of the rRNA genes. Within a genome, putative sequence conversion tracts were found throughout the entire length of each individual rRNA genes and their immediate flanks. Individual conversion events convert only a short sequence tract, and the conversion partners can be any paralogous genes within the genome. Interestingly, the genic sequences undergo much slower divergence than their flanking sequences. Moreover, genomic context and operon organization do not affect rRNA gene homogenization. Thus, gene conversion underlies concerted evolution of bacterial rRNA genes, which normally occurs within genic sequences, and homogenization of flanking regions may result from co-conversion with the genic sequence. Received: 31 March 2000 / Accepted: 15 June 2000     

The Evolution of Hexamerins and the Phylogeny of Insects

The evolutionary relationships among arthropod hemocyanins and insect hexamerins were investigated. A multiple sequence alignment of 12 hemocyanin and 31 hexamerin subunits was constructed and used for studying sequence conservation and protein phylogeny. Although hexamerins and hemocyanins belong to a highly divergent protein superfamily and only 18 amino acid positions are identical in all the sequences, the core structures of the three protein domains are well conserved. Under the assumption of maximum parsimony, a phylogenetic tree was obtained that matches perfectly the assumed phylogeny of the insect orders. An interesting common clade of the hymenopteran and coleopteran hexamerins was observed. In most insect orders, several paralogous hexamerin subclasses were identified that diversified after the splitting of the major insect orders. The dipteran arylphorin/LSP-1-like hexamerins were subject to closer examination, demonstrating hexamerin gene amplification and gene loss in the brachyceran Diptera. The hexamerin receptors, which belong to the hexamerin/hemocyanin superfamily, diverged early in insect evolution, before the radiation of the winged insects. After the elimination of some rapidly or slowly evolving sequences, a linearized phylogenetic tree of the hexamerins was constructed under the assumption of a molecular clock. The inferred time scale of hexamerin evolution, which dates back to the Carboniferous, agrees with the available paleontological data and reveals some previously unknown divergence times among and within the insect orders. Received: 4 August 1997 / Accepted: 29 October 1997     

The genomic organization and evolution of the natural killer immunoglobulin-like receptor (KIR) gene cluster

Natural killer (NK) immunoglobulin-like receptors (KIRs) are a family of polymorphic receptors which interact with specific motifs on HLA class I molecules and modulate NK cytolytic activity. In this study, we analyzed a recently sequenced subgenomic region on chromosome 19q13.4 containing eight members of the KIR receptor repertoire. Six members are clustered within a 100-kb continuous sequence. These genes include a previously unpublished member of the KIR gene family 2DS6, as well as 2DL1, 2DL4, 3DL1, 2DS4, 3DL2, from centromere to telomere. Two additional KIR genes, KIRCI and 2DL3, which may be located centromeric of this cluster were also analyzed. We show that the KIR genes have undergone repeated gene duplications. Diversification between the genes has occurred postduplication primarily as a result of retroelement indels and gene truncation. Using pre- and postduplication Alu sequences identified within these genes as evolutionary molecular clocks, the evolution and duplication of this gene cluster is estimated to have occurred 30–45 million years ago, during primate evolution. A proposed model of the duplication history of the KIR gene family leading to their present organization is presented. Received: 25 November 1999 / Revised: 10 January 2000     

Structure and Origin of a Novel Dimeric Retroposon B1-dID

Here we describe a new short retroposon family of rodents. Like the primate Alu element consisting of two similar monomers, it is dimeric, but the left and right monomers are different and descend from B1 and ID short retroposons, respectively. Such elements (B1-dID) were found in the genomes of Gliridae, Sciuridae, Castoridae, Caviidae, and Hystricidae. Nucleotide sequences of this retroposon can be assigned to several structural variants. Phylogenetic analysis of B1-dID and related sequences suggests a possible scenario of B1-dID evolution in the context of rodent evolution. Received: 30 August 1999 / Accepted: 20 March 2000     

Duplication and Quadruplication of Arabidopsis thaliana Cysteinyl- and Asparaginyl-tRNA Synthetase Genes of Organellar Origin

Two cysteinyl-tRNA synthetases (CysRS) and four asparaginyl-tRNA synthetases (AsnRS) from Arabidopsis thaliana were characterized from genome sequence data, EST sequences, and RACE sequences. For one CysRS and one AsnRS, sequence alignments and prediction programs suggested the presence of an N-terminal organellar targeting peptide. Transient expression of these putative targeting sequences joined to jellyfish green fluorescent protein (GFP) demonstrated that both presequences can efficiently dual-target GFP to mitochondria and plastids. The other CysRS and AsnRSs lack targeting sequences and presumably aminoacylate cytosolic tRNAs. Phylogenetic analysis suggests that the four AsnRSs evolved by repeated duplication of a gene transferred from an ancestral plastid and that the CysRSs also arose by duplication of a transferred organelle gene (possibly mitochondrial). These case histories are the best examples to date of capture of organellar aminoacyl-tRNA synthetases by the cytosolic protein synthesis machinery. Received: 8 October 1999 / Accepted: 23 January 2000