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     

Interspecific Comparison in the Frequency of Concerted Evolution at the Polyubiquitin Gene Locus

The polyubiquitin gene, encoding tandemly repeated multiple ubiquitins, constitutes a uniquitin gene subfamily. It has been demonstrated that polyubiquitin genes are subject to concerted evolution; namely, the individual ubiquitin coding units contained within a polyubiquitin gene are more similar to one another than they are to the ubiquitin coding units in the orthologous gene from other species. However there has been no comprehensive study on the concerted evolution of polyubiquitin genes in a wide range of species, because the relationships (orthologous or paralogous) among multiple polyubiquitin genes from different species have not been extensively analyzed yet. In this report, we present the results of analyzing the nucleotide sequence of polyubiquitin genes of mammals, available in the DDBJ/EMBL/GenBank nucleotide sequence databases, in which we found that there are two groups of polyubiquitin genes in an orthologous relationship. Based on this result, we analyzed the concerted evolution of the polyubiquitin gene in various species and compared the frequency of concerted evolutionary events interspecifically by taking into consideration that the rate of synonymous substitution at the polyubiquitin gene locus may vary depending on species. We found that the concerted evolutionary events in polyubiquitin genes have been more frequent in rats and Chinese hamsters than those in humans, cows, and sheep. The guinea pig polyubiquitin gene was an intermediate example. The frequency of concerted evolution in the mouse gene was unexpectedly low compared to that of other rodent genes. Received: 18 January 2000 / Accepted: 26 April 2000     

Gene conversion and the evolution of euryarchaeal chaperonins: a maximum likelihood-based method for detecting conflicting phylogenetic signals

Recombination is well known as a complicating factor in the interpretation of molecular phylogenies. Here we describe a maximum likelihood sliding window method based on a likelihood ratio test for scanning DNA sequence alignments for regions of incongruent phylogenetic signals, such as those influenced by recombination. Using this method, we identify several instances of gene conversion between paralogous chaperonin genes in euryarchaeote Archaea, many of which are not detected by two other widely used methods. In the Thermococcus/Pyrococcus lineage, where a gene duplication producing a and b paralogues predates the divergence of Thermococcus strains KS-1 and KS-8, gene conversion has homogenized portions of the a and b genes in KS-8 since the divergence of these two strains. A region near the 3′ end of the a and b paralogues in the methanogen Methanobacterium thermoautotrophicum also appears to have undergone gene conversion. We apply the method to two additional test data sets, the argF gene of Neisseria and a set of actin paralogues in maize, and show that it successfully identifies all the recombinant regions that were previously detected with other methods. Our approach is relatively insensitive to the presence of divergent sequences in the alignment, making it ideal for detecting recombination between both closely and distantly related genes.     

Molecular Evolutionary Analysis of the Thiamine-Diphosphate-Dependent Enzyme,Transketolase

Members of the transketolase group of thiamine-diphosphate-dependent enzymes from 17 different organisms including mammals, yeast, bacteria, and plants have been used for phylogenetic reconstruction. Alignment of the amino acid and DNA sequences for 21 transketolase enzymes and one putative transketolase reveals a number of highly conserved regions and invariant residues that are of predicted importance for enzyme activity, based on the crystal structure of yeast transketolase. One particular sequence of 36 residues has some similarities to the nucleotide-binding motif and we designate it as the transketolase motif. We report further evidence that the recP protein from Streptococcus pneumoniae might be a transketolase and we list a number of invariant residues which might be involved in substrate binding. Phylogenies derived from the nucleotide and the amino acid sequences by various methods show a conventional clustering for mammalian, plant, and gram-negative bacterial transketolases. The branching order of the gram-positive bacteria could not be inferred reliably. The formaldehyde transketolase (sometimes known as dihydroxyacetone synthase) of the yeast Hansenula polymorpha appears to be orthologous to the mammalian enzymes but paralogous to the other yeast transketolases. The occurrence of more than one transketolase gene in some organisms is consistent with several gene duplications. The high degree of similarity in functionally important residues and the fact that the same kinetic mechanism is applicable to all characterized transketolase enzymes is consistent with the proposition that they are all derived from one common ancestral gene. Transketolase appears to be an ancient enzyme that has evolved slowly and might serve as a model for a molecular clock, at least within the mammalian clade. Received: 13 September 1995 / Accepted: 14 November 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     

Long-range and targeted ectopic recombination between the two homeologous chromosomes 11 and 12 in Oryza species

Whole genome duplication (WGD) and subsequent evolution of gene pairs have been shown to have shaped the present day genomes of most, if not all, plants and to have played an essential role in the evolution of many eukaryotic genomes. Analysis of the rice (Oryza sativa ssp. japonica) genome sequence suggested an ancestral WGD ～50-70 Ma common to all cereals and a segmental duplication between chromosomes 11 and 12 as recently as 5 Ma. More recent studies based on coding sequences have demonstrated that gene conversion is responsible for the high sequence conservation which suggested such a recent duplication. We previously showed that gene conversion has been a recurrent process throughout the Oryza genus and in closely related species and that orthologous duplicated regions are also highly conserved in other cereal genomes. We have extended these studies to compare megabase regions of genomic (coding and noncoding) sequences between two cultivated (O. sativa, Oryza glaberrima) and one wild (Oryza brachyantha) rice species using a novel approach of topological incongruency. The high levels of intraspecies conservation of both gene and nongene sequences, particularly in O. brachyantha, indicate long-range conversion events less than 4 Ma in all three species. These observations demonstrate megabase-scale conversion initiated within a highly rearranged region located at ～2.1 Mb from the chromosome termini and emphasize the importance of gene conversion in cereal genome evolution.     

Multiple Substitutions Affect the Phylogenetic Utility of Cytochrome b and 12S rDNA Data: Examining a Rapid Radiation in Leporid (Lagomorpha) Evolution

Partial sequences of two mitochondrial genes, the 12S ribosomal gene (739 bp) and the cytochrome b gene (672 bp), were analyzed in hopes of reconstructing the evolutionary relationships of 11 leporid species, representative of seven genera. However, partial cytochrome b sequences were of little phylogenetic value in this study. A suite of pairwise comparisons between taxa revealed that at the intergeneric level, the cytochrome b gene is saturated at synonymous coding positions due to multiple substitution events. Furthermore, variation at the nonsynonymous positions is limited, rendering the cytochrome b gene of little phylogenetic value for assessing the relationships between leporid genera. If the cytochrome b data are analyzed without accounting for these two classes of nucleotides (i.e., synonymous and nonsynonymous sites), one may incorrectly conclude that signal exists in the cytochrome b data. The mitochondrial 12S rRNA gene, on the other hand, has not experienced excessive saturation at either stem or loop positions. Phylogenies reconstructed from the 12S rDNA data support hypotheses based on fossil evidence that African rock rabbits (Pronolagus) are outside of the main leporid stock and that leporids experienced a rapid radiation. However, the molecular data suggest that this radiation event occurred in the mid-Miocene several millions of years earlier than the Pleistocene dates suggested by paleontological evidence. Received: 23 April 1998 / Accepted: 14 May 1998     

Evidence for multiple functional copies of the male sex-determining locus, sry, in African murine rodents

Southern hybridization data suggest that the male sex-determining locus, Sry, is often duplicated in rodents. Here we explore DNA sequence evolution of orthologous and paralogous copies of Sry isolated from six species of African murines. PCR amplification followed by direct sequencing revealed from two to four copies of Sry per species. All copies include a long open reading frame, with a stop codon that coincides closely with the stop codon of the house mouse, Mus musculus, a species known to have a single copy of Sry. A phylogenetic analysis suggests that there are at least seven paralogous copies of Sry in this group of rodents. Putative orthologues are identical; sequence divergence among putative paralogues ranges from 1 to 8% (excluding the CAG repeat), with much lower levels of divergence in the high-mobility group (HMG-box) region than in the C-terminal region. A high proportion of nucleotide substitutions in both regions result in amino-acid replacement. The long open reading frame, conserved HMG-box, and pattern of evolution of the putative paralogues suggest that they are functional. Received: 4 October 1996 / Accepted: 17 January 1997     

Are Noncoding Sequences of Rickettsia prowazekii Remnants of ``Neutralized' Genes?

It has been hypothesized that a large fraction of 24% noncoding DNA in R. prowazekii consists of degraded genes. This hypothesis has been based on the relatively high G+C content of noncoding DNA. However, a comparison with other genomes also having a low overall G+C content shows that this argument would also apply to other bacteria. To test this hypothesis, we study the coding potential in sets of genes, pseudogenes, and intergenic regions. We find that the correlation function and the χ²-measure are clearly indicative of the coding function of genes and pseudogenes. However, both coding potentials make almost no indication of a preexisting reading frame in the remaining 23% of noncoding DNA. We simulate the degradation of genes due to single-nucleotide substitutions and insertions/deletions and quantify the number of mutations required to remove indications of the reading frame. We discuss a reduced selection pressure as another possible origin of this comparatively large fraction of noncoding sequences. Received: 27 December 1999 / Accepted: 5 July 2000     

Comparison of Bombyx mori and Helicoverpa armigera cytoplasmic actin genes provides clues to the evolution of actin genes in insects

The cytoplasmic actin genes BmA3 and BmA4 of Bombyx mori were found clustered in a single genomic clone in the same orientation. As a similar clustering of the two cytoplasmic actin genes Ha3a and Ha3b also occurs in another lepidopteran, Helicoverpa armigera, we analyzed the sequence of the pair of genes from each species. Due to the high conservation of cytoplasmic actins, the coding sequence of the four genes was easily aligned, allowing the detection of similarities in noncoding exon and intron sequences as well as in flanking sequences. All four genes exhibited a conserved intron inserted in codon 117, an original position not encountered in other species. It can thus be postulated that all of these genes derived from a common ancestral gene carrying this intron after a single event of insertion. The comparison of the four genes revealed that the genes of B. mori and H. armigera are related in two different ways: the coding sequence and the intron that interrupts it are more similar between paralogous genes within each species than between orthologous genes of the two species. In contrast, the other (noncoding) regions exhibited the greatest similarity between a gene of one species and a gene of the other species, defining two pairs of orthologous genes, BmA3 and HaA3a on one hand and BmA4 and HaA3b on the other. However, in each species, the very high similarities of the coding sequence and of the single intron that interrupts it strongly suggest that gene conversion events have homogenized this part of the sequence. As the divergence of the B. mori genes was higher than that of the H. armigera genes, we postulated that the gene conversion occurred earlier in the B. mori lineage. This leads us to hypothesize that gene conversion could also be responsible for the original transfer of the common intron to the second gene copy before the divergence of the B. mori and H. armigera lineages.     

Dynamic Insertion–Deletion of Introns in Deuterostome EF-1α Genes

To test the validity of intron–exon structure as a phylogenetic marker, the intron–exon structure of EF-1α genes was investigated for starfish, acornworms, ascidians, larvaceans, and amphioxus and compared with that of vertebrates. Of the 11 distinct intron insertion sites found within the coding regions of the deuterostome EF-1α genes, 7 are shared by several taxa, while the remainder are unique to certain taxa. Examination of the shared introns of the deuterostome EF-1α gene revealed that independent intron loss or intron insertion must have occurred in separate lineages of the deuterostome taxa. Maximum parsimony analysis of the intron–exon data matrix recovered five parsimonious trees (consistency index = 0.867). From this result, we concluded that the intron–exon structure of deuterostome EF-1α has evolved more dynamically than previously thought, rendering it unsuitable as a phylogenetic marker. We also reconstructed an evolutionary history of intron insertion–deletion events on the deuterostome phylogeny, based on several molecular phylogenetic studies. These analyses revealed that the deuterostome EF-1α gene has lost individual introns more frequently than all introns simultaneously.     

The Effect of Recombination on the Accuracy of Phylogeny Estimation

Phylogenetic studies based on DNA sequences typically ignore the potential occurrence of recombination, which may produce different alignment regions with different evolutionary histories. Traditional phylogenetic methods assume that a single history underlies the data. If recombination is present, can we expect the inferred phylogeny to represent any of the underlying evolutionary histories? We examined this question by applying traditional phylogenetic reconstruction methods to simulated recombinant sequence alignments. The effect of recombination on phylogeny estimation depended on the relatedness of the sequences involved in the recombinational event and on the extent of the different regions with different phylogenetic histories. Given the topologies examined here, when the recombinational event was ancient, or when recombination occurred between closely related taxa, one of the two phylogenies underlying the data was generally inferred. In this scenario, the evolutionary history corresponding to the majority of the positions in the alignment was generally recovered. Very different results were obtained when recombination occurred recently among divergent taxa. In this case, when the recombinational breakpoint divided the alignment in two regions of similar length, a phylogeny that was different from any of the true phylogenies underlying the data was inferred.     

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions. Received: 14 November 1996 / Accepted: 7 May 1997     

DNA Repair Systems in Archaea: Mementos from the Last Universal Common Ancestor?

DNA repair in the Archaea is relevant to the consideration of genome maintenance and replication fidelity in the last universal common ancestor (LUCA) from two perspectives. First, these prokaryotes embody a mix of bacterial and eukaryal molecular features. Second, DNA repair proteins would have been essential in LUCA to maintain genome integrity, regardless of the environmental temperature. Yet we know very little of the basic molecular mechanisms of DNA damage and repair in the Archaea in general. Many studies on DNA repair in archaea have been conducted with hyperthermophiles because of the additional stress imposed on their macromolecules by high temperatures. In addition, of the six complete archaeal genome sequences published so far, five are thermophilic archaea. We have recently shown that the hyperthermophile Pyrococcus furiosus has an extraordinarily high capacity for repair of radiation-induced double-strand breaks and we have identified and sequenced several genes involved in DNA repair in P. furiosus. At the sequence level, only a few genes share homology with known bacterial repair genes. For instance, our phylogenetic analysis indicates that archaeal recombinases occur in two paralogous gene families, one of which is very deeply branched, and both recombinases are more closely related to the eukaryotic RAD51 and Dmc1 gene families than to the Escherichia coli recA gene. We have also identified a gene encoding a repair endo/exonuclease in the genomes of several Archaea. The archaeal sequences are highly homologous to those of the eukaryotic Rad2 family and they cluster with genes of the FEN-1 subfamily, which are known to be involved in DNA replication and repair in eukaryotes. We argue that there is a commonality of mechanisms and protein sequences, shared between prokaryotes and eukaryotes for several modes of DNA repair, reflecting diversification from a minimal set of genes thought to represent the genome of the LUCA.     

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     

Cyclostome Hemoglobins Are Possibly Paralogous to Gnathostome Hemoglobins

Dimeric hemoglobins in lampreys are thought to be orthologous to gnathostome hemoglobins, comprising a familiar tetrameric assembly, despite their different subunit interface. To elucidate this contradictory problem, a phylogenetic analysis of vertebrate globins was conducted. The inferred maximum-likelihood trees revealed that the cyclostome hemoglobins are closely related to STAPs, recently identified stellate cell activation-associated proteins, and are paralogous to gnathostome hemoglobins. The quaternary structural difference between cyclostome and gnathostome hemoglobins is well understandable in light of their paralogous relationship.     

Evolution of the nad3-rps12 Gene Cluster in Angiosperm Mitochondria: Comparison of Edited and Unedited Sequences

We have analyzed the nad3-rps12 locus for eight angiosperms in order to compare the utility of mitochondrial DNA and edited mRNA sequences in phylogenetic reconstruction. The two coding regions, containing from 25 to 35 editing sites in the various plants, have been concatenated in order to increase the significance of the analysis. Differing from the corresponding chloroplast sequences, unedited mitochondrial DNA sequences seem to evolve under a quasi-neutral substitution process which undifferentiates the nucleotide substitution rates for the three codon positions. By using complete gene sequences (all codon positions) we found that genomic sequences provide a classical angiosperm phylogenetic tree with a clear-cut grouping of monocotyledons and dicotyledons with Magnoliidae at the basal branch of the tree. Conversely, owing to their low nucleotide substitution rates, edited mRNA sequences were found not to be suitable for studying phylogenetic relationships among angiosperms. Received: 24 January 1996 / Accepted: 5 June 1996