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     

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     

Evolution of DNA Polymerase Families: Evidences for Multiple Gene Exchange Between Cellular and Viral Proteins

A phylogenetic analysis of the five major families of DNA polymerase is presented. Viral and plasmid sequences are included in this compilation along with cellular enzymes. The classification by Ito and Braithwaite (Ito and Braithwaite 1991) of the A, B, C, D, and X families has been extended to accommodate the ``Y family' of DNA polymerases that are related to the eukaryotic RAD30 and the bacterial UmuC gene products. After analysis, our data suggest that no DNA polymerase family was universally conserved among the three biological domains and no simple evolutionary scenario could explain that observation. Furthermore, viruses and plasmids carry a remarkably diverse set of DNA polymerase genes, suggesting that lateral gene transfer is frequent and includes non-orthologous gene displacements between cells and viruses. The relationships between viral and host genes appear very complex. We propose that the gamma DNA polymerase of the mitochondrion replication apparatus is of phage origin and that this gene replaced the one in the bacterial ancestor. Often there was no obvious relation between the viral and the host DNA polymerase, but an interesting exception concerned the family B enzymes: in which ancient gene exchange can be detected between the viruses and their hosts. Additional evidence for horizontal gene transfers between cells and viruses comes from an analysis of the small damage-inducible DNA polymerases. Taken together, these findings suggest a complex evolutionary history of the DNA replication apparatus that involved significant exchanges between viruses, plasmids, and their hosts.     

Identification of Multiple Gypsy LTR-Retrotransposon Lineages in Vertebrate Genomes

Gypsy LTR-retrotransposons have been identified in the genomes of many organisms, but only a small number of vertebrate examples have been reported to date. Here we show that members of this family are likely to be widespread in many vertebrate classes with the possible exceptions of mammals and birds. Phylogenetic analyses demonstrate that although there are several distinct lineages of vertebrate gypsy LTR-retrotransposons, the majority clusters into one monophyletic clade. Groups of fungal, plant, and insect elements were also observed, suggesting horizontal transfer between phyla may be infrequent. However, in contrast to this, there was little evidence to support sister relationships between elements derived from vertebrate and insect hosts. In fact, the majority of the vertebrate elements appeared to be most closely related to a group of gypsy LTR-retrotransposons present within fungi. This implies either that at least one horizontal transmission between these two phyla has occurred previously or that a gypsy LTR-retrotransposon lineage has been lost from insect taxa. Received: 22 December 1998 / Accepted: 6 April 1999     

A Phylogenetic Perspective on Sequence Evolution in Microsatellite Loci

We examined the evolution of the repeat regions of three noncoding microsatellite loci in 58 species of the Polistinae, a subfamily of wasps that diverged over 140 million years ago. A phylogenetic approach allows two new kinds of approaches to studying microsatellite evolution: character mapping and comparative analysis. The basic repeat structure of the loci was highly conserved, but was often punctuated with imperfections that appear to be phylogenetically informative. Repeat numbers evolved more rapidly than other changes in the repeat region. Changes in number of repeats among species seem consistent with the stepwise mutation model, which is based on slippage during replication as the main source of mutations. Changes in repeat numbers can occur even when there are very few tandem repeats but longer repeats, especially perfect repeats led to greater rates of evolutionary change. Species phylogenetically closer to the one from which we identified the loci had longer stretches of uninterrupted repeats and more different motifs, but not longer total repeat regions. The number of perfect repeats increased more often than it decreased. However, there was no evidence that some species have consistently greater numbers of repeats across loci than other species have, once ascertainment bias is eliminated. We also found no evidence for a population size effect posited by one form of the directionality hypothesis. Overall, phylogenetic variation in repeat regions can be explained by adding neutral evolution to what is already known about the mutation process. The life cycle of microsatellites appears to reflect a balance between growth by slippage and degradation by an essentially irreversible accumulation of imperfections. Received: 13 April 1999 / Accepted: 8 September 1999     

Are RNA Viruses Adapting or Merely Changing?

RNA viruses and retroviruses fix substitutions approximately 1 million-fold faster than their hosts. This diversification could represent an inevitable drift under purifying selection, the majority of substitutions being phenotypically neutral. The alternative is to suppose that most fixed mutations are beneficial to the virus, allowing it to keep ahead of the host and/or host population. Here, relative sequence diversification of different proteins encoded by viral genomes is found to be linear. The examples encompass a wide variety of retroviruses and RNA viruses. The smoothness of relative divergence spans quasispeciation following clonal infection, to variation among different isolates of the same virus, to viruses from different species or those associated with different diseases, indicating that the majority of fixed mutations likely reflects drift. This held for both mammalian and plant viruses, indicating that adaptive immunity doesn't necessarily shape the relative accumulation of amino acid substitutions. When compared to their hosts RNA viruses evolution appears conservative. Received: 16 November 1999 / Accepted: 10 March 2000     

Molecular Evolution in a Multisite Nearly Neutral Mutation Model

A simple nearly neutral mutation model of protein evolution was studied using computer simulation assuming a constant population size. In this model, a gene consists of a finite number of codons and there is no recombination within a gene. Each codon has two replacement and one silent sites. The fitness of a gene was determined multiplicatively by amino acids specified by codons (the independent multicodon model). Nucleotide diversity at replacement sites decreases as selection becomes stronger. A reduction of nucleotide diversity at silent sites also occurs as selection intensifies but the magnitude of the reduction is not a monotone function of the intensity of selection. The dispersion index is close to one. The average value of Tajima's and Fu and Li's statistics are negative and their absolute values increases as selection intensifies. However, their powers of detecting selection under the present model were not high unless the number of sites is large or mutation rate is high. The MK test was shown to detect intermediate selection fairly well. For comparison, the house-of-cards model was also investigated and its behavior was shown to be more sensitive to changes of population size than that of the independent multicodon model. The relevance of the present model for explaining protein evolution was discussed comparing its prediction and recent DNA data. Received: 24 May 1999 / Accepted: 17 August 1999     

Characteristics of the Quenching of 9-Aminoacridine Fluorescence by Liposomes Made from Plant Lipids

Several laboratories have determined the surface charge density of membranes utilizing methods based on vesicle-induced quenching of the fluorescence of 9-aminoacridine and its relief by other cations. However, the computational methods by which surface charge density were calculated have not been verified in a model system. In this study, the quenching of 9-aminoacridine fluorescence by liposomes made from varying amounts of digalactosyldiacylglyceride and phosphatidic acid and relief of quenching by salts was examined. Quenching of 9-aminoacridine fluorescence increased with increasing amounts of phosphatidic acid added, independent of the composition of the added liposomes. In certain instances, the computational methods did not yield the surface charge density of the liposomes expected from their composition. However, when the effects of background ionic strength on surface potential were considered, there was a positive correlation between expected and calculated values. Therefore, the data support the contention that changes in the fluorescence of 9-aminoacridine can be used to calculate surface charge density of membranes. Received: 29 November 1999/Revised: 31 July 2000     

The association between HLA-A alleles and an Alu dimorphism near HLA-G

The AluYb8 sequences are a subfamily of short interspersed Alu retroelements that have been amplified within the human genome during recent evolutionary time and are useful polymorphic markers for studies on the origin of human populations. We have identified a new member of the Yb8 subfamily, AluyHG, located between the HLA-H and -G genes and 88-kb telomeric of the highly polymorphic HLA-A gene within the alpha block of the major histocompatibility complex (MHC). The AluyHG element was characterised with a view to examining the association between AluyHG and HLA-A polymorphism and reconstructing the history of the MHC alpha block. A specific primer pair was designed for a simple PCR assay to detect the absence or presence (dimorphism) of the AluyHG element within the DNA samples prepared from a panel of 46 homozygous cell-lines containing complete or recombinant ancestral haplotypes (AH) of diverse ethnic origin and 92 Caucasoid and Asian subjects on which HLA-A typing was available. The AluyHG insertion was most strongly associated with HLA-A2 and, to a lesser degree with HLA-A1, -A3, -A11, and A-19. The gene frequency of the AluyHG insertion for 146 Caucasians and 94 Chinese-Han was 0.30 and 0.32 and there was no significant difference between the observed and expected frequencies. The results of the association studies and the phylogenetic analysis of HLA-A alleles suggest that the AluyHG sequence was integrated within the progenitor of HLA-A2, but has been transferred by recombination to other human ancestral populations. In this regard, the dimorphic AluyHG element is an important diagnostic marker for HLA association studies and could help in elucidating the evolution and functions of the MHC alpha block and polymorphism within and between ancestral haplotypes. Received: 7 December 2000 / Accepted: 28 February 2001     

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     

Speciation Versus Phenotypic Plasticity in Coral Inhabiting Barnacles: Darwin's Observations in an Ecological Context

Speciation and phenotypic plasticity are two extreme strategic modes enabling a given taxon to populate a broad ecological niche. One of the organismal models which stimulated Darwin's ideas on speciation was the Cirripedia (barnacles), to which he dedicated a large monograph. In several cases, including the coral-inhabiting barnacle genera Savignium and Cantellius (formerly Pyrgoma and Creusia, respectively), Darwin assigned barnacle specimens to morphological ``varieties' (as opposed to species) within a genus. Despite having been the subject of taxonomic investigations and revisions ever since, the significance of these varieties has never been examined with respect to host-associated speciation processes. Here we provide evidence from molecular (12S mt rDNA sequences) and micromorphological (SEM) studies, suggesting that these closely related barnacle genera utilize opposite strategies for populating a suite of live-coral substrates. Cantellius demonstrates a relatively low genetic variability, despite inhabiting a wide range of corals. The species C. pallidus alone was found on three coral families, belonging to distinct higher-order classification units. In contrast, Savignium barnacles exhibit large between- and within-species variations with respect to both micromorphology and DNA sequences, with S. dentatum ``varieties' clustering phylogenetically according to their coral host species (all of which are members of a single family). Thus, whereas Savignium seems to have undergone intense host-associated speciation over a relatively narrow taxonomic range of hosts, Cantellius shows phenotypic plasticity over a much larger range. This dichotomy correlates with differences in life-history parameters between these barnacle taxa, including host-infestation characteristics, reproductive strategies, and larval trophic type. Received: 18 January 1999 / Accepted: 9 May 1999     

Testing the Extent of Sequence Similarity Among Viroids, Satellite RNAs, and Hepatitis Delta Virus

A Monte Carlo method was used to test the extent of sequence similarity among viroids, satellite RNAs, and hepatitis delta virus. This analysis revealed that there is insufficient sequence similarity among these pathogens to support the hypothesis that they have a common evolutionary origin. Furthermore, while definite patterns of sequence similarity were observed among some viroids, there was a clear lack of overall similarity, indicating that a monophyletic origin for even this group cannot be reliably supported from sequence data alone. Received: 30 April 1999 / Accepted: 24 August 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     

Spontaneously Opening GABAA Channels in CA1 Pyramidal Neurones of Rat Hippocampus

Spontaneous, single channel, chloride currents were recorded in 48% of cell-attached patches on neurones in the CA1 region of rat hippocampal slices. In some patches, there was more than 1 channel active. They showed outward rectification: both channel conductance and open probability were greater at depolarized than at hyperpolarized potentials. Channels activated by γ-aminobutyric acid (GABA) in silent patches on the same neurones had similar conductance and outward rectification. The spontaneous currents were inhibited by bicuculline and potentiated by diazepam. It was concluded that the spontaneously opening channels were constitutively active, nonsynaptic GABA_A channels. Such spontaneously opening GABA_A channels may provide a tonic inhibitory mechanism in these cells and perhaps in other cells that have GABA_A receptors although not having a GABA_A synaptic input. They may also be a target for clinically useful drugs such as the benzodiazepines. Received: 31 August 1999/Revised: 2 November 1999     

Defining the Core of Nontransferable Prokaryotic Genes: The Euryarchaeal Core

If lateral gene transfer (LGT) has affected all genes over the course of prokaryotic evolution, reconstruction of organismal phylogeny is compromised. However, if a core of genes is immune to transfer, then the evolutionary history of that core might be our most reliable guide to the evolution of organisms. Such a core should be preferentially included in the subset of genes shared by all organisms, but where universally conserved genes have been analyzed, there is too little phylogenetic signal to allow determination of whether or not they indeed have the same history (Hansmann and Martin 2000; Teichmann and Mitchison 1999). Here we look at a more restricted set, 521 homologous genes (COGs) simultaneously present in four sequenced euryarchaeal genomes. Although there is overall little robust phylogenetic signal in this data set, there is, among well-supported trees, strong representation of all three possible four-taxon topologies. ``Informational' genes seem no less subject to LGT than are ``operational genes,' within the euryarchaeotes. We conclude that (i) even in this collection of conserved genes there has been extensive LGT (orthologous gene replacement) and (ii) the notion that there is a core of nontransferable genes (the ``core hypothesis') has not been proven and may be unprovable. Received: 7 November 2000 / Accepted: 20 February 2001     

Electric Field Pulses Can Induce Apoptosis

Injection of electric field pulses of high intensity (kV/cm) and short duration (microsecond range) into a cell suspension results in a temporary increase of the membrane permeability due to a reversible electric breakdown of the cell membrane. Here we demonstrate that application of supercritical field pulses between 4.5 and 8.1 kV/cm strength and 40 μsec duration induce typical features of apoptosis in Jurkat T-lymphoblasts and in HL-60 cells including DNA fragmentation and cleavage of the poly(ADP ribose) polymerase. Apoptosis induction did not depend on the presence of any particular electrolyte in the extracellular medium. However, no apoptosis was observed in solutions without a minimum amount of salt. Apoptotic DNA fragmentation was prevented by the caspase inhibitor zVAD. Received: 16 December 1998/Revised: 24 February 1999     

Molecular Parasites That Evolve Longer Genomes

Molecular parasites that utilize the replication machinery of cells or of in vitro amplification reactions have previously been characterized. By and large, these parasites have been smaller than the viruses or amplicons that gave rise to them. This is likely because shorter genomes can be replicated more quickly. In contrast, we have identified and characterized parasites of an isothermal amplification reaction that are longer than their parental molecules yet replicate much more efficiently. These results raise interesting questions regarding whether the optimal size of replicators reflects a trade-off between the information encoded in a parasite and the information encoded in the machinery replicating that parasite. Received: 9 February 1999 / Accepted: 7 June 1999     

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.     

Evolutionary Analysis of the ErbB Receptor and Ligand Families

We have compared all available deduced protein sequences of the ErbB family of receptors and their ligands. Analysis of the aligned sequences of the receptors indicates that there are some differences in the receptors that are specific to invertebrates. In addition, comparison of the vertebrate ErbB receptors suggest that a gene duplication event generated two ancestral receptors, the ErbB3/ErbB4 precursor and the ErbB1/ErbB2 precursor. Subsequent gene duplications of these precursors generated the four receptors present in mammals. Analysis of the sequences for the known ligands of the ErbB receptors suggests that the vertebrate ligands segregate into the ErbB1 ligands and the ErbB3/ErbB4 ligands, paralleling the evolution of the receptors; however, it is difficult to ascertain any correlation between the invertebrate and the vertebrate ligands. Even though ErbB3 is kinase-impaired, there is significant conservation of the kinase domain within the vertebrate lineage (human, rat, and F. rubripes), suggesting some function for this domain other than kinase activity, such as mediating protein–protein interactions that are involved in receptor dimerization and/or activation of the kinase domain of the heterodimerization partner. To date, no ligand for ErbB2 has been identified, and comparison of the extracellular domains of ErbB2 reveals two regions that are not conserved across the mammalian species. These two regions of divergence align with sequences in ErbB1 that have been shown to be proximal to the amino-terminus and to the carboxyl-terminal region, respectively, of bound EGF. Further, one of these regions contains an insertion, relative to the other members of the mammalian ErbB family, which might affect the ligand binding site and provide a structural basis for this receptor's apparent inability to bind ligand independently. Received: 8 September 1999 / Accepted: 17 January 2000     

An Attempt to Measure the Patchwork Pattern Observed Among Alleles at Major Histocompatibility Complex Loci

By means of simulations and DNA sequence analyses, standardized identity excess (a measure of linkage disequilibrium) between segregating nucleotide sites was studied as an effort to quantify the patchwork pattern among alleles of the major histocompatibility complex loci. It was found that the pattern under selective neutrality, and/or no intralocus recombination does not fit the observed pattern based on DNA sequences. However, the intensity and type of selection and the rate of recombination are difficult to estimate by comparing simulation results with the observed pattern. Received: 10 December 1999 / Accepted: 2 March 2000