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.     

Detecting and characterizing gene conversions between multigene family members

We used a variety of methods to detect known gene conversions in the actin gene families of five angiosperm species, the beta-globin gene families of two primate species, and the Zfx/Zfy gene families of seven mammalian species. Our goal was to devise a working strategy which would allow the analysis of the members of a multigene family in order to determine whether there had been gene conversions between its members, identify the genes involved in the gene conversions, establish the lengths of the converted regions, and determine the polarities of the gene conversions. We show that three phylogenetic methods and the homoplasy test of Maynard Smith and Smith perform relatively poorly on our data sets because the sequences we analyzed had large levels of multiple substitutions. The method of Sawyer, the compatibility method of Jakobsen and Easteal, the partition matrix method of Jakobsen, Wilson, and Easteal, and the co-double method of Balding, Nichols, and Hunt can be used to identify the genes which have been involved in gene conversions. The co-double method is more powerful than other methods but requires orthologous sequences from related species. Compatibility, phylogenetic, and nucleotide substitution distribution statistics methods can be used to identify the location of the converted region(s). Site-by-site compatibility analyses can also be used to identify the direction of the conversion event(s). Combinations of these methods can therefore be used to establish the presence, locations, and polarities of gene conversions between multigene family members.     

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.     

Phylogenetic analysis reveals gene conversions in multigene families of rhizobia

Gene families are an important and intrinsic trait of rhizobial species. These gene copies can participate in non-reciprocal recombination events, also called gene conversions. Gene conversion has diverse roles, but it is usually implicated in the evolution of multigene families. Here, we searched for gene conversions in multigene families of six representative rhizobial genomes. We identified 11 gene families with different numbers of copies, genome location and function in CFN42 and CIAT652 strains of Rhizobium etli, Rhizobium sp NGR234, Mesorhizobium loti MAFF303099, Sinorhizobium meliloti 1021, and Bradyrhizobium japonicum USDA110. Gene conversions were detected by phylogenetic inference in the nifD and nifK gene families in R. etli. Sequence analysis confirmed multiple gene conversions in these two gene families. We suggest that gene conversion events have an important role in homogenizing multigene families in rhizobia.     

Trypanosoma brucei: the extent of conversion in antigen genes may be related to the DNA coding specificity

The boundaries of gene conversion in variant-specific antigen genes have been determined in six clones of Trypanosoma brucei. In each clone, antigenic switching involved interaction between two telomeric members of the AnTat 1.1 multigene family, which share extensive homology throughout their coding regions. All conversion events occurred by substitution of faithful copies of donor sequences. Conversion endpoints were nonrandomly distributed. In four clones, the 5' conversion limit was near the antigen translation initiation codon, while in three clones, the 3' conversion limit was located at the "hinge" between the two major antigen domains. In one case, two segmental conversions were involved in antigen switching. These observations reveal that antigen gene conversion can occur without generating point mutations, and suggest that postrecombinational selection may impose a limit on the number of possible rearrangements within antigen genes.     

Tarsius delta- and beta-globin genes: conversions, evolution, and systematic implications

Comparisons between duplicated genes have shown that gene conversions play an important role in the evolution of multigene families. Previous comparisons have documented in the recently duplicated gamma-fetal globin genes of catarrhine primates, over 15 separate conversions affecting extensive stretches of coding and noncoding sequences. In the present study, delta- and beta- globin genes from a lower primate Tarsius syrichta, and the delta-globin gene of the Asian great ape, Pongo pygmaeus, have been isolated and sequenced. Comparisons of these sequences with other primate delta and beta sequences confirmed a previously reported conversion in an anthropoid ancestor and revealed additional conversions in basal primate, stem haplorhine, tarsier, and early lemur lineages. Conversions found between primate delta- and beta-globin genes contrast with those found in the gamma-genes in that delta-beta conversions appear much less frequently and are more restricted to regions conserved by selection (i.e. coding and 5'-regulatory sequences). These differences indicate that soon after a duplication occurs, conversions can be quite frequent and encompass extensive portions of the duplicated region. With time, sequence differences accumulate, particularly in noncoding regions, and limit both the frequency and size of the conversions. Sequences conserved by selection accumulate differences more slowly and are therefore subject to gene conversions for a longer period of time. Both unconverted and converted sequences were consistent in supporting the placement of tarsier with anthropoids.     

On the evolution of multigene families

Multigene families are classified into three groups: small families as exemplified by hemoglobin genes of mammals; middlesize multigene families, by genes of mammalian histocompatibility antigens; and large multigene families, by variable region genes of immunoglobulins. Facts and theories on these evolving multigene families are reviewed, with special reference to the population genetics of their concerted evolution. It is shown that multigene families are evolving under continued occurrence of unequal (but homologous) crossing-over and gene conversion, and that mechanisms for maintaining genetic variability are totally different from the conventional models of population genetics. Thus, in view of widespread occurrence of multigene families in genomes of higher organisms, the evolutionary theory based mainly on change of gene frequency at each locus would appear to need considerable revision.     

On the divergence of genes in multigene families

Statistical properties of the amount of divergence of genes in multigene families are studied. The model considered is an infinite-site neutral model with unbiased intrachromosomal conversion, unbiased interchromosomal conversion, and recombination. By considering the time back to the most recent common ancestor of two genes, both the probability of identity and the moments of S, the number of sites that differ between two sampled genes, are obtained. We find that if recombination rates are large or conversion is always interchromosomal, then the expectation of S is 4N mu n where N is the population size, mu is the rate of mutation per generation per gene and n is the number of genes in the gene family, as the conversion rates approach zero, the moments of divergence do not approach the moments of divergence with conversion rates equal to zero, and it is possible for a decrease in the rate of intrachromosomal conversion to result in a higher probability of identity, but a greater mean divergence of the two genes.     

Sequence Identity in an Early Chorion Multigene Family Is the Result of Localized Gene Conversion

The multigene families that encode the chorion (eggshell) of the silk moth, Bombyx mori, are closely linked on one chromosome. We report here the isolation and characterization of two segments, totaling 102 kb of genomic DNA, containing the genes expressed during the early period of choriogenesis. Most of these early genes can be divided into two multigene families, ErA and ErB, organized into five divergently transcribed ErA/ErB gene pairs. Nucleotide sequence identity in the major coding regions of the ErA genes was 96%, while nucleotide sequence identity for the ErB major coding regions was only 63%. Selection pressure on the encoded proteins cannot explain this difference in the level of sequence conservation between the ErA and ErB gene families, since when only fourfold redundant codon positions are considered, the divergence within the ErA genes is 8%, while the divergence within the ErB genes (corrected for multiple substitutions at the same site) is 110%. The high sequence identity of the ErA major exons can be explained by sequence exchange events similar to gene conversion localized to the major exon of the ErA genes. These gene conversions are correlated with the presence of clustered copies of the nucleotide sequence GGXGGX, encoding paired glycine residues. This sequence has previously been correlated with gradients of gene conversion that extend throughout the coding and noncoding regions of the High-cysteine (Hc) chorion genes of B. mori. We suggest that the difference in the extent of the conversion tracts in these gene families reflects a tendency for these recombination events to become localized over time to the protein encoding regions of the major exons.     

GABI‐DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock‐out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non‐standard conditions or at the molecular level. Knock‐out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.     

Role of gene duplication in evolution

It is now known that many multigene and supergene families exist in eukaryote genomes: multigene families with uniform copy members like genes for ribosomal RNA, those with variable members like immunoglobulin genes, and supergene families such as those for various growth factor and hormone receptors. Many such examples indicate that gene duplication and subsequent differentiation are extremely important for organismal evolution. In particular, gene duplication could well have been the primary mechanism for the evolution of complexity in higher organisms. Population genetic models for the origin of gene families with diverse functions are presented, in which natural selection favors those genomes with more useful mutants in duplicated genes. Since any gene has a certain probability of degenerating by mutation, success versus failure in acquiring a new gene by duplication may be expressed as the ratio of probabilities of spreading of useful versus detrimental mutations in redundant gene copies. Also examined are the effects of gene duplication on evolution by compensatory advantageous mutations. Results of the analyses show that both natural selection and random drift are important for the origin of gene families. In addition, interaction between molecular mechanisms such as unequal crossing-over and gene conversion, and selection or drift is found to have a large effect on evolution by gene duplication.     

Molecular history of gene conversions in the primate fetal gamma-globin genes. Nucleotide sequences from the common gibbon, Hylobates lar

Comparative and phylogenetic analyses of homologous sequences from closely related species reveal genetic events which have happened in the past and thus provide considerable insight into molecular genetic processes. One such process which has been especially important in the evolution of multigene families is gene conversion. The fetal gamma 1 and gamma 2-globin genes of catarrhine primates (humans, apes, and Old World monkeys) underwent numerous gene conversion events after they arose from a gene duplication event 25-35 million years ago. By including the gamma 1- and gamma 2-globin gene sequences from the common gibbon, Hylobates lar, the present work expands the gamma-globin data set to represent all major groups of hominoid primates. A computer-assisted algorithm is introduced which reveals converted DNA segments and provides results very similar to those obtained by site-by-site evolutionary reconstruction. Both methods provide strong evidence for at least 14 different converted stretches in catarrhine primates as well as five conversions in ancestral lineages. Features of gene conversions generalized from this molecular history are 1) conversions are restricted to regions maintaining high degrees of sequence similarity, 2) one gene may dominate in converting another gene, 3) sequences involved in conversions may accumulate changes more rapidly than expected, and 4) certain elements, such as polypurine/polypyrimidine [Y)n) and (TG)n elements, appear to be hotspots for initiating or terminating conversion events.     

The mutational load of a multigene family with uniform members

The mutational load of a multigene family with uniform members was studied by computer simulations. Two models of selection, truncation and exponential fitness, were examined, by using a simple model of gene conversion. It was found that the load is much smaller than the Haldane-Muller prediction under the truncation selection, and that it becomes approximately equal to the value calculated by the formula, nv(1-q)/(m-nq), where n is the copy number, v is the rate of detrimental mutation per gene copy, m is the truncation point in terms of the number of detrimental genes eliminated, and q is the equilibrium frequency of detrimental mutation. However the equilibrium frequency cannot be analytically obtained. For the exponential fitness model, the load is close to the Haldane-Muller value. When there is no gene conversion, the load becomes larger than the cases with conversion both for the truncation and the exponential fitness models. Thus, gene conversion or other mechanisms that are responsible for contraction-expansion of mutants on chromosomes helps eliminating deleterious mutations occurring in multigene families.     

Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Simulation Study of a Multigene Family, with Special Reference to the Evolution of Compensatory Advantageous Mutations

We investigate the evolution of a multigene family incorporating the forces of drift, mutation, gene conversion, unequal crossing over and selection. The use of simulation studies is required due to the complexity of the model. Selection is modeled in two modes: positive selection as a function of the number of different beneficial alleles and negative selection against deleterious alleles. We assume that gene conversion is unbiased, and that all mutations are initially deleterious. Compensation between mutants creates beneficial and neutral alleles, and allowances are made for compensatory mutations either within or between the members of a multigene family. We find that gene conversion can enhance the rate of acquisition of compensatory advantageous mutations when genes are redundant.     

Saccharomyces cerevisiae PAU genes are induced by anaerobiosis

Saccharomyces cerevisiae PAU genes constitute the largest multigene family in yeast, with 23 members located mainly in subtelomeric regions. The role and regulation of these genes were previously unknown. We detected PAU gene expression during alcoholic fermentation. An analysis of PAU gene regulation using PAU-lacZ fusions and Northern analyses revealed that they were regulated by anaerobiosis. PAU genes display, however, different abilities to be induced by anaerobiosis and this appears to be related to their chromosomal localization; two subtelomeric copies are more weakly inducible than an interstitial one. We show that PAU genes are negatively regulated by oxygen and repressed by haem. Examination of PAU gene expression in rox1Delta and tup1Delta strains indicates that PAU repression by oxygen is mediated by an unknown, haem-dependent pathway, which does not involve the Rox1p anaerobic repressor but requires Tup1p. Given the size of the gene family, PAU genes could be expected to be important during yeast life and some of them probably help the yeast to cope with anaerobiosis.