Evolutionary dynamics of the human endogenous retrovirus family HERV-K inferred from full-length proviral genomes

Several distinct families of endogenous retroviruses exist in the genomes of primates. Most of them are remnants of ancient germ-line infections. The human endogenous retrovirus family HERV-K represents the unique known case of endogenous retrovirus that amplified in the human genome after the divergence of human and chimpanzee lineages. There are two types of HERV-K proviral genomes differing by the presence or absence of 292 bp in the pol-env boundary. Human-specific insertions exist for both types. The analyses shown in the present work reveal that several lineages of type 1 and type 2 HERV-K proviruses remained transpositionally active after the human/chimpanzee split. The data also reflect the important role of mosaic evolution (either by recombination or gene conversion) during the evolutionary history of HERV-K. Received: 5 February 2001 / Accepted: 22 March 2001     

Higher-Order Organization of Subrepeats and the Evolution of Cervid Satellite I DNA

Based on sequence analyses of 17 complete centromeric DNA monomers from ten different deer species, a model is proposed for the genesis, evolution, and genomic organization of cervid satellite I DNA. All cervid satellite I DNA arose from the initial amplification of a 31-bp DNA sequence. These 31-bp subrepeats were organized in a hierarchical fashion as 0.8-kb monomers in plesiometacarpalia deer and 1-kb monomers in telemetacarpalia deer. The higher-order repeat nature of cervid centromeric satellite DNA monomers accounts for their high intragenomic and intraspecific sequence conservation. Such high intraspecific sequence conservation validates the use of a single cervid satellite I DNA monomer from each deer species for interspecific sequence comparisons to elucidate phylogenetic relationships. Also, a specific 0.18-kb tandem duplication was observed in all 1-kb monomers, implying that 1-kb cervid satellite I DNA monomers arose from an unequal crossover event between two similar 0.8-kb ancestral DNA sequences. Received: 28 May 1996 / Accepted: 24 October 1996     

A Molecular Phylogeny of Lilium in the Internal Transcribed Spacer Region of Nuclear Ribosomal DNA

Phylogenetic relationships among 55 species of Lilium, Cardiocrinum giganteum, and Nomocharis saluenensis were inferred from nucleotide sequence variations in the internal transcribed spacer (ITS) regions of 18S–25S nuclear ribosomal DNA. The phylogeny derived from ITS sequences estimated using maximum-likelihood methods indicated that (1) most of the species construct their own clade according to the classification based on morphological features at the section level; (2) section Daurolirion is not independent of Sinomartagon, and it is appropriate to integrate two sections as Sinomartagon; (3) it is appropriate that L. henryi and L. bulbiferum are classified into subsection 6a and Sinomartagon–Daurolirion, respectively; (4) subsection 6b is much closer to Sinomartagon than subsection 6a and Archelirion, and it arose directly from Sinomartagon; and (5) Lilium is much closer to Nomocharis than Cardiocrinum. Phylogenetic estimation using sequences of the ITS region is suitable at the levels of genus, section, and most of subsection. Received: 18 December 1998 / Accepted: 14 March 1999     

The Complete Mitochondrial DNA Sequence of the Pig (Sus scrofa)

The complete mitochondrial genome sequence of the pig, Sus scrofa, was determined. The length of the sequence presented is 16,679 nucleotides. This figure is not absolute, however, due to pronounced heteroplasmy caused by variable numbers of the motif GTACACGTGC in the control region of different molecules. A phylogenetic study was performed on the concatenated amino acid and nucleotide sequences of 12 protein-coding genes of the mitochondrial genome. The analysis identified the pig (Suiformes) as a sister group of a cow/whale clade, making Artiodactyla paraphyletic. The split between pig and cow/whale was molecularly dated at 65 million years before present. Received: 2 December 1997 / Accepted: 20 February 1998     

Variable and Invariable DNA Repeat Characters Revealed by Taxonprint Approach Are Useful for Molecular Systematics

A specially optimized restriction analysis of highly repetitive DNA elements, called DNA taxonprint, was applied for phylogenetic study of primates and lizards. It was shown that electrophoretic bands of DNA repeats revealed by the taxonprint technique have valuable properties for molecular systematics. Approximately half of taxonprint bands (TB) are invariable and do not disappear from the genomes during evolution or change spontaneously. Presumably these invariable bands are restriction fragments of dispersed DNA repeats. Another group represents variable taxonprint bands that differ even between closely related species. These variable bands are probably represented by tandem DNA repeats and could be used as species-specific markers. It was shown that taxonprint bands are independent characters since the appearance of a new taxonprint band does not change the previous band pattern. Phylogenetic reconstruction carried out on taxonprint data demonstrated that this approach could be of general utility for molecular systematics and species identification. Received: 12 January 1998 / Accepted: 16 May 1998     

Concerted Evolution of a Highly Repetitive DNA Family in Eptatretidae (Cyclostomata, Agnatha) Implies Specifically Differential Homogenization and Amplification Events in Their Germ Cells

In eight hagfish species, it is known that chromosome elimination occurs during early embryogenesis, and some highly repetitive DNA families, restricted to germ cells, have been isolated. One of these families, ``EEEo2,' has been isolated as DNA fragments by restriction enzyme analyses from Eptatretus okinoseanus and E. cirrhatus. In this study, EEEo2 sequences were isolated from germline DNA in E. burgeri, Paramyxine sheni, and P. atami using PCR methods. Sequence analysis revealed that these sequences are intraspecifically homogeneous, except in E. burgeri, and are interspecifically conserved with heterogeneity. The intraspecific sequence variability tends to decrease as the copy number increases. These results indicate that EEEo2 has evolved in a concerted manner. Moreover, an ancestral repeating motif consisting of triplicate subrepeats was deduced. These results suggest that EEEo2 arose as an initial amplification of this subrepeat and has evolved by saltatory replication. Phylogenetic analyses suggested the possibility that EEEo2 in E. okinoseanus and E. cirrhatus has been subjected to strong homogenizing forces for concerted evolution, whereas the force is weak in E. burgeri. In addition, EEEo2 in P. sheni and P. atami appear to have been incompletely subjected to these forces. Chromosomal in situ hybridization experiments revealed that EEEo2 sequences were located along almost their entire length of several heterochromatic chromosomes that are restricted to germ cells. These chromosomes are disposed to form a secondary association during the first meiotic metaphases, except in P. sheni. This chromosomal distribution may promote a concerted mode of sequence evolution in both nonhomologous chromosomes and homologous chromosomes and reflect the differential driving forces between species. Received: 17 April 1999 / Accepted: 10 September 1999     

Mitochondrial DNA Phylogeny of the Family Cichlidae: Monophyly and Fast Molecular Evolution of the Neotropical Assemblage

A mitochondrial DNA (mtDNA) phylogeny of cichlid fish is presented for the most taxonomically inclusive data set compiled to date (64 taxa). 16S rDNA data establish with confidence relationships among major lineages of cichlids, with a general pattern congruent with previous morphological studies and less inclusive molecular phylogenies based on nuclear genes. Cichlids from Madagascar and India are the most basal groups of the family Cichlidae and sister to African–Neotropical cichlids. The cichlid phylogeny suggests drift-vicariance events, consistent with the fragmentation of Gondwana, to explain current biogeographic distributions. Important phylogenetic findings include the placement of the controversial genus Heterochromis basal among African cichlids, the South American genus Retroculus as the most basal taxon of the Neotropical cichlid assemblage, and the close relationship of the Neotropical genera Cichla with Astronotus rather than with the crenicichlines. Based on a large number of South American genera, the Neotropical cichlids are defined as a monophyletic assemblage and shown to harbor significantly higher levels of genetic variation than their African counterparts. Relative rate tests suggest that Neotropical cichlids have experienced accelerated rates of molecular evolution. But these high evolutionary rates were significantly higher among geophagine cichlids. Received: 18 September 1998 / Accepted: 16 December 1998     

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.     

Evolution of Histone H4 and H3 Genes in Different Ciliate Lineages

The histones H4 are known as highly conserved proteins. However, in ciliates a high degree of variation was found compared both to other eukaryotes and between the ciliate species. To date, only H4 histones of species belonging to two distantly related classes have been investigated. In order to obtain more detailed information on histone H4 variation in ciliates we undertook a comprehensive sequence analysis of PCR-amplified internal H4 fragments from 12 species belonging to seven out of the nine currently recognized ciliate classes. In addition, we used PCR primers to amplify longer fragments of H3 and H4 genes including the intergenic region. The encoded amino acid sequences reveal a high number of differences when compared with those of other eukaryotes and the ciliate species investigated. Furthermore, in some species H4 gene variants were detected, which result in amino acid differences. The greatest number of substitutions and insertions found was in the amino terminal region of the H4 histones. However, all sequences possess a conserved region corresponding to those of all other eukaryotic H4 histones. The histone gene variations were used to reconstruct phylogenetic relationships. The tree from our data matches perfectly with the ribosomal RNA data: The heterotrichs, which were considered as a late branching lineage, diverge at the base of the ciliate tree and groups formerly thought to represent ancestral lineages now appear as highly derived ciliates. Received: 4 April 1997 / Accepted: 1 August 1997     

Extensive Amplification and Transposition of a Novel Repetitive Element, Xstir, Together with Its Terminal Inverted Repeat in the Evolution of Xenopus

A DNA fragment containing short tandem repeat sequences (approximately 86-bp repeat) was isolated from a Xenopus laevis cDNA library. Southern blot and in situ hybridization analyses revealed that the repeat was highly dispersed in the genome and was present at approximately 1 million copies per haploid genome. We named this element Xstir (Xenopus short tandemly and invertedly repeating element) after its arrangement in the genome. The majority of the genomic Xstir sequences were digested to monomer and dimer sizes with several restriction enzymes. Their sequences were found to be highly homogeneous and organized into tandem arrays in the genome. Alignment analyses of several known sequences showed that some of the Xstir-like sequences were also organized into interspersed inverted repeats. The inverted repeats consisted of an inverted pair of two differently modified Xstirs separated by a short insert. In addition, these were framed by another novel inverted repeat (Xstir-TIR). The Xstir-TIR sequence was also found at the ends of tandem Xstir arrays. Furthermore, we found that Xstir-TIR was linked to a motif characterizing the T2 family which belonged to a vertebrate MITE (miniature inverted-repeat transposable element) family, suggesting the importance of Xstir-TIR for their amplification and transposition. The present study of 11 anuran and 2 urodele species revealed that Xstir or Xstir-like sequences were extensively amplified in the three Xenopus species. Genomic Xstir populations of X. borealis and X. laevis were mutually indistinguishable but significantly different from that of X. tropicalis. Received: 5 April 2000 / Accepted: 3 August 2000     

Mitochondrial DNA Phylogeny and the Evolution of Host-Plant Use in Palearctic Chrysolina (Coleoptera, Chrysomelidae) Leaf Beetles

The genus Chrysolina consists of specialized phytophagous leaf-beetles (Coleoptera, Chrysomelidae) with feed on several plant families. There is no explicit phylogenetic hypothesis available for this genus, which includes 65 subgenera and more than 400 species with a wide distribution. We obtained 839-bp sequence data from the 16S rDNA and cytochrome oxidase subunit I (COI) mitochondrial genes. Thirty Chrysolina taxa representing eight host–plant affiliations, two species of the closely related genus Oreina, and two outgroups were sampled. These data sets were used separately and combined to obtain the mitochondrial cladogram of the group using maximum-parsimony and maximum-likelihood criteria. The results were compared to current proposals for Chrysolina systematics that are based on morphological, ecological, and karyological data. The trees obtained were in the most part congruent with the proposed ancestral association of Chrysolina to Lamiaceae based on chromosome number in several lineages. A minimum of five host-plant switches from the ancestral state inferred at the family level and two at the subclass level suggests the absence of parallel evolution of beetles and their host plants. Another switch leading to oligophagy at the family level was deduced to have occurred in the lineage of the subgenus Chrysolina s.str. Received: 22 May 1998 / Accepted: 16 September 1998     

Molecular Characterization of the Principal Symbiotic Bacteria of the Weevil Sitophilus oryzae: A Peculiar G + C Content of an Endocytobiotic DNA

The principal intracellular symbiotic bacteria of the cereal weevil Sitophilus oryzae were characterized using the sequence of the 16S rDNA gene (rrs gene) and G + C content analysis. Polymerase chain reaction amplification with universal eubacterial primers of the rrs gene showed a single expected sequence of 1,501 bp. Comparison of this sequence with the available database sequences placed the intracellular bacteria of S. oryzae as members of the Enterobacteriaceae family, closely related to the free-living bacteria, Erwinia herbicola and Escherichia coli, and the endocytobiotic bacteria of the tsetse fly and aphids. Moreover, by high-performance liquid chromatography, we measured the genomic G + C content of the S. oryzae principal endocytobiotes (SOPE) as 54%, while the known genomic G + C content of most intracellular bacteria is about 39.5%. Furthermore, based on the third codon position G + C content and the rrs gene G + C content, we demonstrated that most intracellular bacteria except SOPE are A + T biased irrespective of their phylogenetic position. Finally, using the hsp60 gene sequence, the codon usage of SOPE was compared with that of two phylogenetically closely related bacteria: E. coli, a free-living bacterium, and Buchnera aphidicola, the intracellular symbiotic bacteria of aphids. Taken together, these results show a peculiar and distinctly different DNA composition of SOPE with respect to the other obligate intracellular bacteria, and, combined with biological and biochemical data, they elucidate the evolution of symbiosis in S. oryzae. Received: 8 September 1997 / Accepted: 24 October 1997     

Ribosomal External and Internal Transcribed Spacers: Combined Use in the Phylogenetic Analysis of Medicago (Leguminosae)

We have estimated the potential phylogenetic utility of the ribosomal external transcribed spacer (ETS) from the nuclear ribosomal region. The ETS was sequenced from 13 annual Medicago (Fabaceae) species upstream a highly conserved motive which was found among many different organisms. In the genus Medicago, the ETS was found to evolve 1.5 times faster than the internal transcribed spacer and to be 1.5 times more informative. Reduced ribosomal maturation process constraints on ETS are proposed to explain the different evolutionary rates between the two spacers. Maximal phylogenetic resolution and support was obtained when the two spacers were analyzed together. No incongruence between the two spacers was found and ETS appears to be a valuable source of information for solidifying ITS plant phylogeny. The phylogeny obtained in Medicago suggests that none of the three subsections included in the study is monophyletic. Received: 15 April 1997 / Accepted: 29 July 1997     

Ligand-binding domains in vitellogenin receptors and other LDL-receptor family members share a common ancestral ordering of cysteine-rich repeats

Insect vitellogenin and yolk protein receptors (VgR/YPR) are newly discovered members of the low-density lipoprotein receptor (LDLR) family, which is characterized by a highly conserved arrangement of repetitive modular elements homologous to functionally unrelated proteins. The insect VgR/YPRs are unique in having two clusters of complement-type cysteine-rich (class A) repeats or modules, with five modules in the first cluster and seven in the second cluster, unlike classical LDLRs which have a single seven-module cluster, vertebrate VgRs and very low density lipoprotein receptors (VLDLR) which have a single eight-module cluster, and LDLR-related proteins (LRPs) and megalins which have four clusters of 2–7, 8, 10, and 11 modules. Alignment of clusters across subfamilies by conventional alignment programs is problematic because of the repetitive nature of the component modules which may have undergone rearrangements, duplications, and deletions during evolution. To circumvent this problem, we ``fingerprinted' each class A module in the different clusters by identifying those amino acids that are both relatively conserved and relatively unique within the cluster. Intercluster reciprocal comparisons of fingerprints and aligned sequences allowed us to distinguish four cohorts of modules reflecting shared recent ancestry. All but two of the 57 modules examined could be assigned to one of these four cohorts designated A, B, C, and D. Alignment of clusters based on modular cohorts revealed that all clusters are derived from a single primordial cluster of at least seven modules with a consensus arrangement of CDCADBC. All extant clusters examined are consistent with this consensus, though none matches it perfectly. This analysis also revealed that the eight-module clusters in vertebrate VgRs, insect VgR/YPRs, and LRP/megalins are not directly homologous with one another. Assignment of modules to cohorts permitted us to properly align 32 class A clusters from all four LDLR subfamilies for phylogenetic analysis. The results revealed that smaller one-cluster and two-cluster members of the family did not originate from the breakup of a large two-cluster or four-cluster receptor. Similarly, the LRP/megalins did not arise from the duplication of a two-cluster insect VgR/YPR-like progenitor. Rather, it appears that the multicluster receptors were independently constructed from the same single-cluster ancestor. Received: 16 January 1997 / Accepted: 21 August 1997     

A Reexamination of the Phylogenetic Position of Callimico (Primates) Incorporating New Mitochondrial DNA Sequence Data

The New World monkeys are divided into two main groups, Callitrichidae and Cebidae. Callimico goeldii shares traits with both the Cebidae and the Callitrichidae. Recent morphological phyletic studies generally place Callimico as the most basal member of the Callitrichidae. In contrast, genetic studies (immunological, restriction fragment, and sequence data) have consistently placed Callimico somewhere within the Callitrichidae, not basal to this clade. A DNA sequence data set from the terminal 236 codons of the mitochondrial ND4 gene and the tRNA^His, tRNA^Ser, and tRNA^Leu genes was generated to clarify the position of Callimico. The sequences of 887 base pairs were analyzed by maximum-parsimony, neighbor-joining, and maximum-likelihood methods. The results of these various methods are generally congruent and place Callimico within the Callitrichidae between the marmosets (Callithrix and Cebuella) and the tamarins (Saguinus and Leontopithecus). Combined analyses of all suitable nuclear and mitochondrial gene sequences confirm the position of Callimico between the marmosets and the tamarins. As available molecular evidence indicates that Callimico is more closely related to the marmosets than to the tamarins, a reconsideration of the morphological evidence in light of the consensus tree from DNA sequence analyses is warranted. The marmosets and tamarins share four morphological characters (loss of the third molar, loss of the hypocone, reduced body size, reproductive twinning). Dwarfism may have evolved repeatedly among the Callitrichidae. It is well-known that the loss of a character can occur many times independently. The reproduction of marmosets and tamarins is extremely specialized and it is difficult to imagine that this complex and unique twinning system evolved separately in marmosets and tamarins. However, it is possible that a secondary reversal to single offspring took place in Callimico. Received: 20 March 1997 / Accepted: 17 December 1997     

DNA Polymerase C of the Thermophilic Bacterium Thermus aquaticus: Classification and Phylogenetic Analysis of the Family C DNA Polymerases

Bacterial family C DNA polymerases (DNA pol IIIs), the major chromosomal replicative enzymes, have been provisionally classified based on primary sequences and domain structures into three classes: class I (Escherichia coli DNA pol C-type), class II (Bacillus subtilis DNA pol C-type), and class III (cyanobacterial DNA pol C-type), respectively. We have sequenced the structural gene encoding the DNA pol C catalytic subunit of the thermophilic bacterium Thermus aquaticus. This gene, designated the Taq DNA pol C gene, contains a 3660-bp open reading frame which specifies a polypeptide of molecular weight of 137,388 daltons. Comparative sequence analyses revealed that Taq DNA pol C is a class I family C DNA polymerase. The Taq DNA pol C is most closely related to the Deinococcus radiodurans DNA pol C. Although a phylogenetic tree based on the class I family C DNA pols is still in the provisional stage, some important conclusion can be drawn. First, the high-G+C and the low-G+C Gram-positive bacteria are not monophyletic. Second, the low-G+C Gram-positive bacteria contain multigenes of family C DNA pols (classes I and II). Third, the cyanobacterial family C DNA pol, classified as class III because it is encoded by a split gene, forms a group with the high-G+C Gram-positive bacteria. Received: 7 October 1998 / Accepted: 12 January 1999     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Full-Length L1 Elements Have Arisen Recently in the Same 1-kb Region of the Gorilla and Human Genomes

New copies of the mammalian retrotransposon L1 arise in the germline at an undetermined rate. Each new L1 copy appears at a specific evolutionary time point that can be estimated by phylogenetic analysis. In humans, the active L1 sequence L1.2 resides at the genomic locus LRE1. Here we analyzed the region surrounding the LRE1 locus in humans and gorillas to determine the evolutionary history of the region and to estimate the age of L1.2. We found that the region was composed of an ancient L1, L1Hs-Lrg, which was significantly divergent from all other L1 sequences available in the databases. We also determined that L1.2 was absent from the gorilla genome and arose in humans after the divergence of gorilla and human lineages. In the gorilla LRE1 region, we discovered a different full-length L1 element, L1Gg-1, which was allelic and present at a high gene frequency in gorillas but absent from other primates. We determined the nucleotide sequence of L1Gg-1 and found that it was 98% identical to L1.2, suggesting a close relationship between active L1s in gorillas and humans. Received: 28 December 1997 / Accepted: 20 March 1998