Ectopic gene conversions in bacterial genomes.

We characterized the gene conversions found between the duplicated genes of 75 bacterial genomes from five species groups (archaea, nonpathogenic and pathogenic firmicutes, and nonpathogenic and pathogenic proteobacteria). The number of gene conversions is positively correlated with the size of multigene families and the size of multigene families is not significantly different between pathogenic and nonpathogenic taxa. However, gene conversions occur twice as frequently in pathogenic species as in nonpathogenic species. Comparisons between closely related species also indicate a trend towards increased gene conversion in pathogenic species. Whereas the length of the conversions is positively correlated with flanking sequence similarity in all five groups, these correlations are smaller for pathogenic firmicutes and proteobacteria than for nonpathogenic firmicutes and proteobacteria. These results are consistent with our previous work on E. coli genomes and suggest that pathogenic bacteria allow recombination between more divergent gene sequences. This higher permissiveness is likely adaptive because it allows them to generate more genetic variability.     

Distributions of exons and introns in the human genome

The human genome is revisited using exon and intron distribution profiles. The 26,564 annotated genes in the human genome (build October, 2003) contain 233,785 exons and 207,344 introns. On average, there are 8.8 exons and 7.8 introns per gene. About 80% of the exons on each chromosome are < 200 bp in length. < 0.01% of the introns are < 20 bp in length and < 10% of introns are more than 11,000 bp in length. These results suggest constraints on the splicing machinery to splice out very long or very short introns and provide insight to optimal intron length selection. Interestingly, the total length in introns and intergenic DNA on each chromosome is significantly correlated to the determined chromosome size with a coefficient of correlation r = 0.95 and r = 0.97, respectively. These results suggest their implication in genome design.     

Ectopic Gene Conversions in Four <Emphasis Type="BoldItalic">Escherichia coli</Emphasis> Genomes: Increased Recombination in Pathogenic Strains

We characterized the ectopic gene conversions in the genomes of the K-12 MG1655, O157:H7 Sakai, O157:H7 EDL933, and CFT073 strains of E coli. Compared to the three pathogenic strains, the K-12 strain has a much smaller number of gene families, its gene families contain fewer genes, and gene conversions are less frequent. Whereas the three pathogenic strains have gene conversions covering hundreds of nucleotides when their flanking regions have as little as 50% similarity, flanking region similarity of at least 94% on both sides of the converted region is required to observe conversions of more than 87 nucleotides in the K-12 strain. Recombination is therefore more frequent and requires less sequence similarity in the three pathogenic strains than in K-12. This higher recombination level might be due to mutations in some of their mismatch-repair genes. In contrast with the gene conversions present in the yeast genome, the gene conversions found in the E. coli genomes do not occur more frequently between duplicated genes that are close to one another than between duplicated genes that are far apart and are randomly distributed along the length of the genes. In E. coli, gene conversions are not more frequent near the origin of replication. However, they do occur more frequently near the terminus of replication of the Sakai genome, where multigene family members are more abundant. This suggests that, in E. coli, gene conversions occur randomly between genes located in different chromosomal locations or located on different copies of the multiple chromosomes found in E. coli cells.     

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     

Characterization of the gene conversions between the multigene family members of the yeast genome

Stanley Sawyer's gene conversion detection method, implemented in his GENECONV computer program, was used to detect and characterize the gene conversions between the multigene family members of the yeast genome. This method gave different gene conversion frequencies and size distribution for gene families with two members and multigene families with more than two members. The 69 gene conversions detected in multigene families with more than two members occur at a frequency of 7.8% gene conversion/pair of genes compared and have an average size of 173+/-220 nucleotides. Larger gene conversions are found only between more similar genes, the genes involved in gene conversions are distributed almost randomly among the 16 yeast chromosomes, and the frequency of gene conversions increases as the distance between repeated genes decreases. In contrast to previous studies, no relationship was observed between the level of expression of a gene and its involvement in gene conversions. These analyses also suggest that gene conversions might occur by different mechanisms in closely linked genes and unlinked genes. The excess of converted regions at the 3? end of unlinked genes suggests that recombination with incomplete cDNA molecules is the main mechanism responsible for gene conversions between such genes.     

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.     

Intron-exon structures of eukaryotic model organisms.

To investigate the distribution of intron-exon structures of eukaryotic genes, we have constructed a general exon database comprising all available intron-containing genes and exon databases from 10 eukaryotic model organisms: Homo sapiens, Mus musculus, Gallus gallus, Rattus norvegicus, Arabidopsis thaliana, Zea mays, Schizosaccharomyces pombe, Aspergillus, Caenorhabditis elegans and Drosophila. We purged redundant genes to avoid the possible bias brought about by redundancy in the databases. After discarding those questionable introns that do not contain correct splice sites, the final database contained 17 102 introns, 21 019 exons and 2903 independent or quasi-independent genes. On average, a eukaryotic gene contains 3.7 introns per kb protein coding region. The exon distribution peaks around 30-40 residues and most introns are 40-125 nt long. The variable intron-exon structures of the 10 model organisms reveal two interesting statistical phenomena, which cast light on some previous speculations. (i) Genome size seems to be correlated with total intron length per gene. For example, invertebrate introns are smaller than those of human genes, while yeast introns are shorter than invertebrate introns. However, this correlation is weak, suggesting that other factors besides genome size may also affect intron size. (ii) Introns smaller than 50 nt are significantly less frequent than longer introns, possibly resulting from a minimum intron size requirement for intron splicing.     

Organization of the human hepatocyte growth factor-encoding gene.

Human genomic phage libraries were screened for the human hepatocyte growth factor (HGF)-encoding gene (HGF) using a cDNA encoding the human protein as a probe. Characterization of the clones revealed that this gene is composed of 18 exons interrupted by 17 introns spanning approx. 70 kb. The first exon contains the 5'-untranslated region and the signal peptide. The next ten exons encode the alpha-chain which contains four kringle structures. Each kringle domain is encoded by two exons as observed in other kringle-containing proteins. The twelfth exon contains the short spacer region between the alpha- and beta-chains and the remaining six exons comprise the beta-chain. The beta-chain is structurally similar to the catalytic domains of serine proteases; amino acid substitutions in the active site were found. The organization of the HGF gene is highly homologous to those of the serine proteases involved in blood coagulation and fibrinolysis, especially with that of plasminogen. This suggests that the human HGF gene is evolutionally related to these genes.     

Sequence-Dependent Gene Conversion: Can Duplicated Genes Diverge Fast Enough to Escape Conversion?

Conversion between duplicated genes limits their independent evolution. Models in which conversion frequencies decrease as genes diverge are examined to determine conditions under which genes can "escape" further conversion and hence escape from a gene family. A review of results from various recombination systems suggests two classes of sequence-dependence models: (1) the "k-hit" model in which conversion is completely inactivated by a few (k) mutational events, such as the insertion of a mobile element, and (2) more general models where conversion frequency gradually declines as genes diverge through the accumulation of point mutants. Exact analysis of the k-hit model is given and an approximate analysis of a more general sequence-dependent model is developed and verified by computer simulation. If mu is the per nucleotide mutation rate, then neutral duplicated genes diverging through point mutants are likely to escape conversion provided 2 mu/lambda much greater than 0.1, where lambda is the conversion rate between identical genes. If 2 mu/lambda much less than 0.1, the expected number of conversions before escape increases exponentially so that, for biological purposes, the genes never escape conversion. For single mutational events sufficient to block further conversions, occurring at rate nu per copy per generation, many conversions are expected if 2 nu/lambda much less than 1, while the genes essentially evolve independently if 2 nu/lambda much greater than 1. Implications of these results for both models of concerted evolution and the evolution of new gene functions via gene duplication are discussed.     

Concerted evolution of the primate immunoglobulin alpha-gene through gene conversion.

We determined four nucleotide sequences of the hominoid immunoglobulin alpha (C alpha) genes (chimpanzee C alpha 2, gorilla C alpha 2, and gibbon C alpha 1 and C alpha 2 genes), which made possible the examination of gene conversions in all hominoid C alpha genes. The following three methods were used to detect gene conversions: 1) phenetic tree construction; 2) detection of a DNA segment with extremely low variability between duplicated C alpha genes; and 3) a site by site search of shared nucleotide changes between duplicated C alpha genes. Results obtained from method 1 indicated a concerted evolution of the duplicated C alpha genes in the human, chimpanzee, gorilla, and gibbon lineages, while results obtained from method 2 suggested gene conversions in the human, gorilla, and gibbon C alpha genes. With method 3 we identified clusters of shared nucleotide changes between duplicated C alpha genes in human, chimpanzee, gorilla, and gibbon lineages, and in their hypothetical ancestors. In the present study converted regions were identified over the entire C alpha gene region excluding a few sites in the coding region which have escaped from gene conversion. This indicates that gene conversion is a general phenomenon in evolution, that can be clearly observed in non-functional regions.