The Ghost of Selection Past: Rates of Evolution and Functional Divergence of Anciently Duplicated Genes

The duplication of genes and even complete genomes may be a prerequisite for major evolutionary transitions and the origin of evolutionary novelties. However, the evolutionary mechanisms of gene evolution and the origin of novel gene functions after gene duplication have been a subject of many debates. Recently, we compiled 26 groups of orthologous genes, which included one gene from human, mouse, and chicken, one or two genes from the tetraploid Xenopus and two genes from zebrafish. Comparative analysis and mapping data showed that these pairs of zebrafish genes were probably produced during a fish-specific genome duplication that occurred between 300 and 450 Mya, before the teleost radiation (Taylor et al. 2001). As discussed here, many of these retained duplicated genes code for DNA binding proteins. Different models have been developed to explain the retention of duplicated genes and in particular the subfunctionalization model of Force et al. (1999) could explain why so many developmental control genes have been retained. Other models are harder to reconcile with this particular set of duplicated genes. Most genes seem to have been subjected to strong purifying selection, keeping properties such as charge and polarity the same in both duplicates, although some evidence was found for positive Darwinian selection, in particular for Hox genes. However, since only the cumulative pattern of nucleotide substitutions can be studied, clear indications of positive Darwinian selection or neutrality may be hard to find for such anciently duplicated genes. Nevertheless, an increase in evolutionary rate in about half of the duplicated genes seems to suggest that either positive Darwinian selection has occurred or that functional constraints have been relaxed at one point in time during functional divergence. Received: 4 January 2001 / Accepted: 29 March 2001     

Phylogenetic Timing of the Fish-Specific Genome Duplication Correlates with the Diversification of Teleost Fish

For many genes, ray-finned fish (Actinopterygii) have two paralogous copies, where only one ortholog is present in tetrapods. The discovery of an additional, almost-complete set of Hox clusters in teleosts (zebrafish, pufferfish, medaka, and cichlid) but not in basal actinopterygian lineages (Polypterus) led to the formulation of the fish-specific genome duplication hypothesis. The phylogenetic timing of this genome duplication during the evolution of ray-finned fish is unknown, since only a few species of basal fish lineages have been investigated so far. In this study, three nuclear genes (fzd8, sox11, tyrosinase) were sequenced from sturgeons (Acipenseriformes), gars (Semionotiformes), bony tongues (Osteoglossomorpha), and a tenpounder (Elopomorpha). For these three genes, two copies have been described previously teleosts (e.g., zebrafish, pufferfish), but only one orthologous copy is found in tetrapods. Individual gene trees for these three genes and a concatenated dataset support the hypothesis that the fish-specific genome duplication event took place after the split of the Acipenseriformes and the Semionotiformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. If these three genes were duplicated during the proposed fish-specific genome duplication event, then this event separates the species-poor early-branching lineages from the species-rich teleost lineage. The additional number of genes resulting from this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish.[Reviewing Editor: Martin Kreitman]     

Evolutionary Patterns of Gene Families Generated in the Early Stage of Vertebrates

In this paper we have analyzed 49 vertebrate gene families that were generated in the early stage of vertebrates and/or shortly before the origin of vertebrates, each of which consists of three or four member genes. We have dated the first (T₁) and second (T₂) gene duplications of 26 gene families with 3 member genes. The means of T₁ (594 mya) and T₂ (488 mya) are largely consistent to a well-cited version of two-round (2R) genome duplication theory. Moreover, in most cases, the time interval between two successive gene duplications is large enough that the fate of duplicate genes generated by the first gene duplication was likely to be determined before the second one took place. However, the phylogenetic pattern of 23 gene families with 4 members is complicated; only 5 of them are predicted by 2R model, but 11 families require an additional gene (or genome) duplication. For the rest (7 families), at least one gene duplication event had occurred before the divergence between vertebrate and Drosophila, indicating a possible misleading of the 4:1 rule (member gene ratio between vertebrates and invertebrates). Our results show that Ohno's 2R conjecture is valid as a working hypothesis for providing a most parsimonious explanation. Although for some gene families, additional gene duplication is needed, the credibility of the third genome duplication (3R) remains to be investigated. Received: 13 December 1999 / Accepted: 7 April 2000     

A Phylogenomic Study of the OCTase Genes in Pseudomonas syringae Pathovars: The Horizontal Transfer of the argK–tox Cluster and the Evolutionary History of OCTase Genes on Their Genomes

Phytopathogenic Pseudomonas syringae is subdivided into about 50 pathovars due to their conspicuous differentiation with regard to pathogenicity. Based on the results of a phylogenetic analysis of four genes (gyrB, rpoD, hrpL, and hrpS), Sawada et al. (1999) showed that the ancestor of P. syringae had diverged into at least three monophyletic groups during its evolution. Physical maps of the genomes of representative strains of these three groups were constructed, which revealed that each strain had five rrn operons which existed on one circular genome. The fact that the structure and size of genomes vary greatly depending on the pathovar shows that P. syringae genomes are quite rich in plasticity and that they have undergone large-scale genomic rearrangements. Analyses of the codon usage and the GC content at the codon third position, in conjunction with phylogenomic analyses, showed that the gene cluster involved in phaseolotoxin synthesis (argK–tox cluster) expanded its distribution by conducting horizontal transfer onto the genomes of two P. syringae pathovars (pv. actinidiae and pv. phaseolicola) from bacterial species distantly related to P. syringae and that its acquisition was quite recent (i.e., after the ancestor of P. syringae diverged into the respective pathovars). Furthermore, the results of a detailed analysis of argK [an anabolic ornithine carbamoyltransferase (anabolic OCTase) gene], which is present within the argK–tox cluster, revealed the plausible process of generation of an unusual composition of the OCTase genes on the genomes of these two phaseolotoxin-producing pathovars: a catabolic OCTase gene (equivalent to the orthologue of arcB of P. aeruginosa) and an anabolic OCTase gene (argF), which must have been formed by gene duplication, have first been present on the genome of the ancestor of P. syringae; the catabolic OCTase gene has been deleted; the ancestor has diverged into the respective pathovars; the foreign-originated argK–tox cluster has horizontally transferred onto the genomes of pv. actinidiae and pv. phaseolicola; and hence two copies of only the anabolic OCTase genes (argK and argF) came to exist on the genomes of these two pathovars. Thus, the horizontal gene transfer and the genomic rearrangement were proven to have played an important role in the pathogenic differentiation and diversification of P. syringae. Received: 22 May 2001 / Accepted: 26 September 2001     

Histone H1 Genes and Histone Gene Clusters in the Genus Drosophila

Whereas the genomes of many organisms contain several nonallelic types of linker histone genes, one single histone H1 type is known in Drosophila melanogaster that occurs in about 100 copies per genome. Amplification of H1 gene sequences from genomic DNA of wild type strains of D. melanogaster from Oregon, Australia, and central Africa yielded numerous clones that all exhibited restriction patterns identical to each other and to those of the known H1 gene sequence. Nucleotide sequences encoding the evolutionarily variable domains of H1 were determined in two gene copies of strain Niamey from central Africa and were found to be identical to the known H1 sequence. Most likely therefore, the translated sequences of D. melanogaster H1 genes do not exhibit intragenomic or intergenomic variations. In contrast, three different histone H1 genes were isolated from D. virilis and found to encode proteins that differ remarkably from each other and from the H1 of D. melanogaster and D. hydei. About 40 copies of H1 genes are organized in the D. virilis genome with copies of core histone genes in gene quintets that were found to be located in band 25F of chromosome 2. Another type of histone gene cluster is present in about 15 copies per genome and contains a variable intergenic sequence instead of an H1 gene. The H1 heterogeneity in D. virilis may have arisen from higher recombination rates than occur near the H1 locus in D. melanogaster and might provide a basis for formation of different chromatin subtypes. Received: 2 March 2000 / Accepted: 1 June 2000     

Characterization of the Hydra Lamin and Its Gene: A Molecular Phylogeny of Metazoan Lamins

We report sequences for nuclear lamins from the teleost fish Danio and six invertebrates. These include two cnidarians (Hydra and Tealia), one priapulid, two echinoderms, and the cephalochordate Branchiostoma. Combining these results with earlier data on Drosophila, Caenorhabditis elegans, and various vertebrates, the following conclusions on lamin evolution can be drawn. First, all invertebrate lamins resemble in size the vertebrate B-type lamin. Second, all lamins described previously for amphibia, birds and mammals as well as the first lamin of a fish, characterized here, show a cluster of 7 to 12 acidic residues in the tail domain. Since this acidic cluster is absent from all invertebrate lamins including that of the cephalochordate Branchiostoma, it was acquired with the vertebrate lineage. The larger A-type lamin of differentiated cells must have arisen subsequently by gene duplication and insertion of an extra exon. This extra exon of the vertebrate A-lamins is the only major change in domain organization in metazoan lamin evolution. Third, the three introns of the Hydra and Priapulus genes correspond in position to the last three introns of vertebrate B-type lamin genes. Thus the entirely different gene organization of the C. elegans and Drosophila Dmo genes seems to reflect evolutionary drift, which probably also accounts for the fact that C. elegans has the most diverse lamin sequence. Finally we discuss the possibility that two lamin types, a constitutively expressed one and a developmentally regulated one, arose independently on the arthropod and vertebrate lineages. Received: 4 February 1999 / Accepted: 1 April 1999     

Evolution of MADS-box gene induction by FLO/LFY genes

Some MADS-box genes function as floral homeotic genes. The Arabidopsis LFY gene is a positive regulator of floral homeotic genes, and homologs of the FLO/LFY gene family in other angiosperms and gymnosperms are likely to have a similar function. To investigate the origin of the floral homeotic gene regulatory cascade involving the FLO/LFY gene, FLO/LFY homologs were cloned from a leptosporangiate fern (Ceratopteris richardii), two eusporangiate ferns (Angiopteris lygodiifolia and Botrychium multifidum var. robustum), three fern allies (Psilotum nudum, Equisetum arvense, and Isoetes asiatica), and a moss (Physcomitrella patens). The FLO/LFY gene phylogenetic tree indicates that both duplication and loss of FLO/LFY homologs occurred during the course of vascular plant evolution. The expression patterns of the Ceratopteris LFY genes (CrLFY1 and 2) were assessed. CrLFY1 expression was prominent in tissues including shoot tips and circinate reproductive leaves, but very weak in other tissues examined. Expression of CrLFY2 was also prominent in tissues, including shoot tips and circinate reproductive leaves. These patterns of expression are dissimilar to that of any Ceratopteris MADS-box gene previously reported, suggesting that the induction of MADS-box genes by FLO/LFY is not established at the stage of ferns. Received: 4 January 2001 / Accepted: 28 February 2001     

Phylogenomic Analysis of the α Proteasome Gene Family from Early-Diverging Eukaryotes

We employed a phylogenomic approach to study the evolution of α subunits of the proteasome gene family from early diverging eukaryotes. BLAST similarity searches of the Giardia lamblia genome identified all seven α proteasome genes characteristic of eukaryotes from the crown group. In addition, a PCR strategy for the amplification of multiple α subunit sequences generated single α proteasome products for representatives of the Kinetoplastida (Leishmania major), the Parabasalia (Trichomonas vaginalis), and the Microsporidia (Vairimorpha sp., Nosema sp., Endoreticulata sp., and Spraguea lophii). The kinetoplastid Trypanosoma cruzi and the eukaryote crown group Acanthamoeba castellanii yielded two distinct α proteasome genes each. The presence of seven distinct α proteasome genes in G. lamblia, one of the earliest-diverging eukaryotes, indicates that the α proteasome gene family evolved rapidly from a minimum of one gene in Archaea to seven or more in Eukarya. Results from the phylogenomic analysis are consistent with the idea that the Diplomonida (as represented by G. lamblia), the Kinetoplastida, the Parabasalia, and the Microsporidia diverged after the duplication events that originated the α proteasome gene family. A model for the early origin and evolution of the proteasome gene family is presented. Received: 14 February 2000 / Accepted: 14 August 2000     

The evolution of teleost pigmentation and the fish‐specific genome duplication

Teleost fishes have evolved a unique complexity and diversity of pigmentation and colour patterning that is unmatched among vertebrates. Teleost colouration is mediated by five different major types of neural‐crest derived pigment cells, while tetrapods have a smaller repertoire of such chromatophores. The genetic basis of teleost colouration has been mainly uncovered by the cloning of pigmentation genes in mutants of zebrafish Danio rerio and medaka Oryzias latipes. Many of these teleost pigmentation genes were already known as key players in mammalian pigmentation, suggesting partial conservation of the corresponding developmental programme among vertebrates. Strikingly, teleost fishes have additional copies of many pigmentation genes compared with tetrapods, mainly as a result of a whole‐genome duplication that occurred 320–350 million years ago at the base of the teleost lineage, the so‐called fish‐specific genome duplication. Furthermore, teleosts have retained several duplicated pigmentation genes from earlier rounds of genome duplication in the vertebrate lineage, which were lost in other vertebrate groups. It was hypothesized that divergent evolution of such duplicated genes may have played an important role in pigmentation diversity and complexity in teleost fishes, which therefore not only provide important insights into the evolution of the vertebrate pigmentary system but also allow us to study the significance of genome duplications for vertebrate biodiversity.     

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     

The Cephalopod Loligo bleekeri Mitochondrial Genome: Multiplied Noncoding Regions and Transposition of tRNA Genes

We previously reported the sequence of a 9260-bp fragment of mitochondrial (mt) DNA of the cephalopod Loligo bleekeri [J. Sasuga et al. (1999) J. Mol. Evol. 48:692–702]. To clarify further the characteristics of Loligo mtDNA, we have sequenced an 8148-bp fragment to reveal the complete mt genome sequence. Loligo mtDNA is 17,211 bp long and possesses a standard set of metazoan mt genes. Its gene arrangement is not identical to any other metazoan mt gene arrangement reported so far. Three of the 19 noncoding regions longer than 10 bp are 515, 507, and 509 bp long, and their sequences are nearly identical, suggesting that multiplication of these noncoding regions occurred in an ancestral Loligo mt genome. Comparison of the gene arrangements of Loligo, Katharina tunicata, and Littorina saxatilis mt genomes revealed that 17 tRNA genes of the Loligo mt genome are adjacent to noncoding regions. A majority (15 tRNA genes) of their counterparts is found in two tRNA gene clusters of the Katharina mt genome. Therefore, the Loligo mt genome (17 tRNA genes) may have spread over the genome, and this may have been coupled with the multiplication of the noncoding regions. Maximum likelihood analysis of mt protein genes supports the clade Mollusca + Annelida + Brachiopoda but fails to infer the relationships among Katharina, Loligo, and three gastropod species. Received: 9 May 2001 / Accepted: 3 October 2001     

Rapid Concerted Evolution via Gene Conversion at the <Emphasis Type="BoldItalic">Drosophila hsp70</Emphasis> Genes

We analyzed nucleotide variation in the hsp70 genes of Drosophila melanogaster (five genes) and D. simulans (four genes) to characterize the homogenizing and diversifying roles of gene conversion in their evolution. Gene conversion within and between the 87A7 and 87C1 gene clusters homogenize the hsp70 coding regions; in both D. melanogaster and D. simulans, same-cluster paralogues are virtually identical, and large intercluster conversion tracts diminish 87A7/87C1 divergence. Same-cluster paralogues share many polymorphisms, consistent with frequent intracluster conversion. Shared polymorphism is highly biased toward silent variation; homogenizing conversion interacts with purifying selection. In contrast to the coding regions, some hsp70 flanking regions show conversion-mediated diversification. Strong reductions of nucleotide variability and linkage disequilibria among conversion-mediated sites in hsp70Ab and hsp70Bb alleles sampled from a single natural population are consistent with a selective sweep. Comparison of the D. melanogaster and D. simulans hsp70 genes reveals whole-family fixed differences, consistent with rapid propagation of novel mutations among duplicate genes. These results suggest that the homogenizing and diversifying roles of conversion interact to drive dynamic concerted evolution of the hsp70 genes. Received: 25 June 2001 / Accepted: 10 October 2001     

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     

Sucrose Phosphate Synthase Genes in Plants Belong to Three Different Families

We present phylogenetic analyses to demonstrate that there are three families of sucrose phosphate synthase (SPS) genes present in higher plants. Two data sets were examined, one consisting of full-length proteins and a second larger set that covered a highly conserved region including the 14-3-3 binding region and the UDPGlu active site. Analysis of both datasets showed a well supported separation of known genes into three families, designated A, B, and C. The genomic sequences of Arabidopsis thaliana include a member in each family: two genes on chromosome 5 belong to Family A, one gene on chromosome 1 to Family B, and one gene on chromosome 4 to Family C. Each of three Citrus genes belong to one of the three families. Intron/exon organization of the four Arabidopsis genes differed according to phylogenetic analysis, with members of the same family from different species having similar genomic organization of their SPS genes. The two Family A genes on Arabidopsis chromosome 5 appear to be due to a recent duplication. Analysis of published literature and ESTs indicated that functional differentiation of the families was not obvious, although B family members appear not to be expressed in roots. B family genes were cloned from two Actinidia species and southern analysis indicated the presence of a single gene family, which contrasts to the multiple members of Family A in Actinidia. Only two family C genes have been reported to date. Received: 17 April 2001 / Accepted: 27 August 2001     

Limitations of compositional approach to identifying horizontally transferred genes

Genes with atypical G+C content and pattern of codon usage in a certain genome are possibly of exotic origin, and this idea has been applied to identify horizontal events. In this way, it was postulated that a total of 755 genes in the E. coli genome are relics of horizontal events after the divergence of E. coli from the Salmonella lineage 100 million years ago (Lawrence and Ochman, 1998). In this paper we propose a new way to study sequence composition more thoroughly. We found that although the 755 genes differ in composition from other genes in the E. coli genome, the difference is minor. If we accepted that these genes are horizontally transferred, then (1) it would be more likely that they were transferred from genomes evolutionarily closely related to E. coli; but (2) the dating method used by Lawrence and Ochman (1997, 1998) largely underestimated the average age of introduced sequences in the E. coli genome, in particular, most of the 755 genes should be introduced into E. coli before, instead of after, the divergence of E. coli from the Salmonella lineage. Our study reveals that atypical G+C content and pattern of codon usage are not reliable indicators of horizontal gene transfer events. Received: 27 September 2000 / Accepted: 9 April 2001     

Extensive Gene Duplication in the Early Evolution of Animals Before the Parazoan–Eumetazoan Split Demonstrated by G Proteins and Protein Tyrosine Kinases from Sponge and Hydra

To know whether genes involved in cell–cell communication typical of multicellular animals dramatically increased in concert with the Cambrian explosion, the rapid evolutionary burst in the major groups of animals, and whether these genes exist in the sponge lacking cell cohesiveness and coordination typical of eumetazoans, we have carried out cloning of the G-protein α subunit (Gα) and the protein tyrosine kinase (PTK) cDNAs from Ephydatia fluviatilis (freshwater sponge) and Hydra magnipapillata strain 105 (hydra). We obtained 13 Gα and 20 PTK cDNAs. Generally animal gene families diverged first by gene duplication (subtype duplication) that gave rise to diverse subtypes with different primary functions, followed by further gene duplication in the same subtype (isoform duplication) that gave rise to isoform genes with virtually identical function. Phylogenetic trees of Gα and PTK families including cDNAs from sponge and hydra revealed that most of the present-day subtypes had been established in the very early evolution of animals before the parazoan–eumetazoan split, the earliest branching among the extant animal phyla, by extensive subtype duplication: for PTK and Gα families, 23 and 9 subtype duplications were observed in the early stage before the parazoan–eumetazoan split, respectively, and after that split, only 2 and 1 subtype duplications were found, respectively. After the separation from arthropods, vertebrates underwent frequent isoform duplications before the fish–tetrapod split. Furthermore, rapid amino acid changes appear to have occurred in concert with the extensive subtype duplication and isoform duplication. Thus the pattern of gene diversification during animal evolution might be characterized by bursts of gene duplication interrupted by considerably long periods of silence, instead of proceeding gradually, and there might be no direct link between the Cambrian explosion and the extensive gene duplication that generated diverse functions (subtypes) of these families. Received: 4 November 1998 / Accepted: 17 November 1998     

Gene Duplication and Gene Conversion in the Caenorhabditis elegans Genome

A comprehensive analysis of duplication and gene conversion for 7394 Caenorhabditis elegans genes (about half the expected total for the genome) is presented. Of the genes examined, 40% are involved in duplicated gene pairs. Intrachromosomal or cis gene duplications occur approximately two times more often than expected. In general the closer the members of duplicated gene pairs are, the more likely it is that gene orientation is conserved. Gene conversion events are detectable between only 2% of the duplicated pairs. Even given the excesses of cis duplications, there is an excess of gene conversion events between cis duplicated pairs on every chromosome except the X chromosome. The relative rates of cis and trans gene conversion and the negative correlation between conversion frequency and DNA sequence divergence for unconverted regions of converted pairs are consistent with previous experimental studies in yeast. Three recent, regional duplications, each spanning three genes are described. All three have already undergone substantial deletions spanning hundreds of base pairs. The relative rates of duplication and deletion may contribute to the compactness of the C. elegans genome. Received: 30 July 1998 / Accepted: 12 October 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     

Biglycan-Like Extracellular Matrix Genes of Agnathans and Teleosts

Biglycan and decorin are two members of a family of small extracellular matrix proteoglycans characterized by the presence of 10 leucine-rich repeats and one or two attachment sites for glucosaminoglycans. Both have thus far been described only from tetrapod species, mainly mammals. Because the extracellular matrix has played an important part in the evolution of Metazoa, the phylogeny of its components is of considerable interest. In this study, biglycan-like (BGL) cDNA sequences have been obtained from two teleost (Oreochromis cichlid and zebrafish) and two lamprey species. The analysis of the sequences suggests that, like tetrapods, the lampreys possess two types of proteoglycans, both of which are biglycan-like; decorin-like proteoglycans could not be identified in these species. The genes specifying these two types apparently arose by duplication in the lamprey lineage after its divergence from gnathostomes. The two teleost species possess a BGL proteoglycan and a bona fide decorin. The BGL proteoglycan is highly divergent from the tetrapod biglycan and related to the BGL proteoglycans of the lamprey. Hence, although the duplication generating the ancestors of biglycan and decorin genes occurred after the divergence of agnathans but before the emergence of teleosts, only decorin acquired its characteristic properties in the bony fishes. The BGL gene presumably turned into a typical biglycan only in the tetrapod lineages. The presumed acquisitions of new functions appear to have been accompanied by changes in the evolutionary rate. Received: 13 April 2000 / Accepted: 4 July 2000