Population Genetics of an Expanding Family of Mobile Genetic Elements

A model of an expanding family of dispersed repetitive DNA was studied. Based on the previous result of the model of duplicative transposition, an approximate solution to give allelism and identify coefficients as functions of time was obtained, and theoretical predictions were verified by Monte Carlo experiments. The results show that, even if the copy number per genome increases very rapidly, allelism and identity coefficients may take a long time to reach equilibrium. The changes of allelism and allelic identity are similar to that of homozygosity at an ordinary single locus, whereas that of nonallelic identity can be much slower, particularly when the copy number per genome is large. Thus, many existing families of highly repetitive sequences may represent nonequilibrium states for nonallelic identity. The present model may be extended to include other evolutionary forces such as gene conversion or the recurrent insertion from normal gene copies.     

A strong deletion bias in nonallelic gene conversion

Gene conversion is the unidirectional transfer of genetic information between orthologous (allelic) or paralogous (nonallelic) genomic segments. Though a number of studies have examined nucleotide replacements, little is known about length difference mutations produced by gene conversion. Here, we investigate insertions and deletions produced by nonallelic gene conversion in 338 Drosophila and 10,149 primate paralogs. Using a direct phylogenetic approach, we identify 179 insertions and 614 deletions in Drosophila paralogs, and 132 insertions and 455 deletions in primate paralogs. Thus, nonallelic gene conversion is strongly deletion-biased in both lineages, with almost 3.5 times as many conversion-induced deletions as insertions. In primates, the deletion bias is considerably stronger for long indels and, in both lineages, the per-site rate of gene conversion is orders of magnitudes higher than that of ordinary mutation. Due to this high rate, deletion-biased nonallelic gene conversion plays a key role in genome size evolution, leading to the cooperative shrinkage and eventual disappearance of selectively neutral paralogs.     

Frequent independent duplicative transpositions activate a single VSG gene.

The expression of several surface antigen genes in Trypanosoma brucei is mediated by the duplicative transposition of a basic-copy variant surface glycoprotein (VSG) gene into an expression site. We determined that the appearance of variant 118, in a parasitemia, resulted from at least four independent duplicative transpositions of the same VSG 118 gene. Variants 117 and 118 both appeared at specific periods but resulted from multiple independent activations. Antigenic variants thus occur in an ordered manner. We show that in the duplicative transpositions of VSG genes, the ends of the transposed segments were homologous between the basic copy and the expression site. Sequences other than the previously reported 70-base-pair (bp) repeats could be involved. In one variant, 118 clone 1, the homology was between a sequence previously transposed into the expression site and a sequence located 6 kilobases upstream of the VSG 118 gene. In variant 118b the homology was presumably in 70-bp repeat arrays, while in a third 118 variant yet another sequence was involved. The possibility that the 70-bp repeats are important in the initial steps of the recombinational events was illustrated by a rearrangement involving a 70-bp repeat array. The data provide strong evidence for the notion that gene conversion mediates the duplicative transposition of VSG genes. We discuss a model that explains how the process of duplicative transposition can occur at random and still produce an ordered appearance of variants.     

Gene conversions mediating antigenic variation in Trypanosoma brucei can occur in variant surface glycoprotein expression sites lacking 70-base-pair repeat sequences.

African trypanosomes undergo antigenic variation of their variant surface glycoprotein (VSG) coat to avoid immune system-mediated killing by their mammalian host. An important mechanism for switching the expressed VSG gene is the duplicative transposition of a silent VSG gene into one of the telomeric VSG expression sites of the trypanosome, resulting in the replacement of the previously expressed VSG gene. This process appears to be a gene conversion reaction, and it has been postulated that sequences within the expression site may act to initiate and direct the reaction. All bloodstream form expression sites contain huge arrays (many kilobase pairs) of 70-bp repeat sequences that act as the 5' boundary of gene conversion reactions involving most silent VSG genes. For this reason, the 70-bp repeats seemed a likely candidate to be involved in the initiation of switching. Here, we show that deletion of the 70-bp repeats from the active expression site does not affect duplicative transposition of VSG genes from silent expression sites. We conclude that the 70-bp repeats do not appear to function as indispensable initiation sites for duplicative transposition and are unlikely to be the recognition sequence for a sequence-specific enzyme which initiates recombination-based VSG switching.     

Nonallelic gene conversion is not GC-biased in Drosophila or primates

Gene conversion is the unidirectional transfer of genetic information between allelic (orthologous) or nonallelic (paralogous) DNA segments. Recently, there has been much interest in understanding how gene conversion shapes the nucleotide composition of the genomic landscape. A widely held hypothesis is that gene conversion is universally GC-biased. However, direct experimental evidence of this hypothesis is limited to a single study of meiotic crossovers in yeast. Although there have been a number of indirect studies of gene conversion, evidence of GC-biased replacements gathered from such studies can also be attributed to positive selection, which has the same evolutionary dynamics as biased gene conversion. Here, we apply a direct phylogenetic approach to examine nucleotide replacements produced by nonallelic gene conversion in Drosophila and primate genomes. We find no evidence for GC-biased gene conversion in either lineage, suggesting that previously observed GC biases may be due to positive selection rather than to biased gene conversion.     

Factors affecting ectopic gene conversion in mice

Duplicated genes and repetitive sequences are distributed throughout the genomes of complex organisms. The homology between related sequences can promote nonallelic (ectopic) recombination, including gene conversion and reciprocal exchange. Resolution of these events can result in translocations, deletions, or other harmful rearrangements. In yeast, ectopic recombination between sequences on nonhomologous chromosomes occurs at high frequency. Because the mammalian genome is replete with duplicated sequences and repetitive elements, high levels of ectopic exchange would cause aneuploidy and genome instability. To understand the factors regulating ectopic recombination in mice, we evaluated the effects of homology length on gene conversion between unlinked sequences in the male germline. Previously, we found high levels of gene conversion between lacZ transgenes containing 2557 bp of homology. We report here that genetic background can play a major role in ectopic recombination; frequency of gene conversion was reduced by more than an order of magnitude by transferring the transgenes from a CF1 strain background to C57BL/6J. Additionally, conversion rates decreased as the homology length decreased. Sequences sharing 1214 bp of sequence identity underwent ectopic conversion less frequently than a pair sharing 2557 bp of identity, while 624 bp was insufficient to catalyze gene conversion at significant levels. These results suggest that the germline recombination machinery in mammals has evolved in a way that prevents high levels of ectopic recombination between smaller classes of repetitive sequences, such as the Alu family. Additionally, genomic location appeared to influence the availability of sequences for ectopic recombination. Received: 12 September 1997 / Accepted: 29 December 1997     

Mass-murder deletion of 19 ORFs from Saccharomyces cerevisiae chromosome XI

In eukaryote genomes, there are many kinds of gene families. Gene duplication and conversion are sources of the evolution of gene families, including those with uniform members and those with diverse functions. Population genetics theory on identity coefficients among gene members of a gene family shows that the balance between diversification by mutation, and homogenization by unequal crossing over and gene conversion, is important. Also, evolution of new functions is due to gene duplication followed by differentiation. Positive selection is necessary for the evolution of novel functions. However, many examples of current gene families suggest that both drift and selection are at work on their evolution.     

Genetic Differentiation of Transposable Elements under Mutation and Unbiased Gene Conversion

A model is developed to predict the extent of genetic differentiation in a family of transposable elements under the combined effects of genetic drift, transposition, mutation and unbiased gene conversion. The model is based on simplifying assumptions that are valid when transposition is always to new sites and copy number per site is low. In the absence of gene conversion, the degree of differentiation as measured by the probability of identity of different elements is the same as at a single locus with the same mutation rate but in a population of effective size Nc/2, where N is the population size and c is the number of copies per individual. The inclusion of unbiased gene conversion does not significantly change this result. If, as seems to be the case, families of transposable elements are relatively homogeneous, then the model implies either that mutation rates for transposable elements are much lower than at comparable single-copy loci or that some other force, such as natural selection or biased gene conversion, is at work. Transposition is a very ineffective force for homogenizing a family of transposable elements.     

Gene conversion as a mechanism for antigenic variation in trypanosomes

Expression of the gene coding for the trypanosome AnTat 1.1 surface antigen is linked to the duplicative transposition of a basic copy (BC) of this gene to an expression site. In two trypanosome clones successively derived from AnTat 1.1 (AnTat 1.10 and AnTat 1.1B) we found evidence that gene conversions are involved in the transformation of the AnTat 1.1 transposed element into the two new surface antigen coding sequences. Although the three resultant mRNAs--AnTat 1.1, 1.10, and 1.1B--are different, they still share large homologies. Two of them, AnTat 1.1 and 1.1B, code for surface coats that are indistinguishable by conventional serological techniques, whereas AnTat 1.10 has been found different by the same methods. The three genomic rearrangements involve two of the five members of the AnTat 1.1 gene family. These two members are both located in unstable telomeric regions similar to the expression site, each in a different orientation with respect to the DNA terminus. We have concluded that the duplicative transposition is achieved by a gene conversion that may affect variable lengths of the same silent genes, and that different members of the same surface antigen gene family can contribute to the diversification of the antigen repertoire.     

Trypanosome variant surface glycoprotein genes expressed early in infection

We have studied further the genes for trypanosomal variant surface glycoproteins expressed during a chronic infection of rabbits with Trypanosoma brucei, strain 427. We show that there are three closely related chromosomal-internal isogenes for VSG 121; expression of one of these genes is accompanied by the duplicate transposition of the gene to a telomeric expression site, also used by other chromosome-internal VSG genes. The 3' end of the 121 gene is replaced during transposition with another sequence, also found in the VSG mRNAs of two other variants. We infer that an incoming VSG gene duplicate recombines with the resident gene in the expression site and may exchange ends in this process. The extra expression-linked copy of the 121 gene is lost when another gene enters the expression site. However, when the telomeric VSG gene 221 is activated without duplication the extra 121 gene copy is inactivated without detectable alterations in or around the gene. We have also analysed the VSG genes expressed very early when trypanosomes are introduced into rats or tissue culture. The five genes identified in 24 independent switching events were all found to be telomeric genes and we calculate that the telomeric 1.8 gene has a 50% chance of being activated in this trypanosome strain when the trypanosome switches the VSG that is synthesized. We argue that the preferential expression of telomeric VSG genes is due to two factors: first, some telomeric genes reside in an inactive expression site, that can be reactivated; second, telomeric genes can enter an active expression site by a duplicative telomere conversion and this process occurs more frequently than the duplicative transposition of chromosome-internal genes to an expression site.     

The evolution of mobile DNAs: when will transposons create phylogenies that look as if there is a master gene?

Some families of mammalian interspersed repetitive DNA, such as the Alu SINE sequence, appear to have evolved by the serial replacement of one active sequence with another, consistent with there being a single source of transposition: the "master gene." Alternative models, in which multiple source sequences are simultaneously active, have been called "transposon models." Transposon models differ in the proportion of elements that are active and in whether inactivation occurs at the moment of transposition or later. Here we examine the predictions of various types of transposon model regarding the patterns of sequence variation expected at an equilibrium between transposition, inactivation, and deletion. Under the master gene model, all bifurcations in the true tree of elements occur in a single lineage. We show that this property will also hold approximately for transposon models in which most elements are inactive and where at least some of the inactivation events occur after transposition. Such tree shapes are therefore not conclusive evidence for a single source of transposition.     

Transposition, amplification, and divergence in the origin of the DNF15 loci, a polymorphic repetitive sequence family on chromosomes 1 and 3

The loci DNF15S1 and DNF15S2 are members of a small repetitive sequence family at discrete chromosomal locations, namely, 1p36 and 3p21, respectively. Studies of the structure, arrangement, and interrelations of the family suggest that the single copy on chromosome 3 is the original member and that this gave rise to the several members on chromosome 1 by transposition, partial duplication, and amplification. Several restriction fragment length polymorphisms have been discovered at the DNF15S1 locus and these have been assigned to the different subfamilies of the repeat at this locus. The existence of these RFLPs, and the nonallelic restriction site variation also found in this sequence family, suggests that transposition and amplification occurred as discrete events. We sequenced across the ancient junction between chromosomes 1 and 3 and noted features which might explain the mechanics of the transposition and amplification events.     

A double-strand break in a chromosomal LINE element can be repaired by gene conversion with various endogenous LINE elements in mouse cells

A double-strand break (DSB) in the mammalian genome has been shown to be a very potent signal for the cell to activate repair processes. Two different types of repair have been identified in mammalian cells. Broken ends can be rejoined with or without loss or addition of DNA or, alternatively, a homologous template can be used to repair the break. For most genomic sequences the latter event would involve allelic sequences present on the sister chromatid or homologous chromosome. However, since more than 30% of our genome consists of repetitive sequences, these would have the option of using nonallelic sequences for homologous repair. This could have an impact on the evolution of these sequences and of the genome itself. We have designed an assay to look at the repair of DSBs in LINE-1 (L1) elements which number 10(5) copies distributed throughout the genome of all mammals. We introduced into the genome of mouse epithelial cells an L1 element with an I-SceI endonuclease site. We induced DSBs at the I-SceI site and determined their mechanism of repair. We found that in over 95% of cases, the DSBs were repaired by an end-joining process. However, in almost 1% of cases, we found strong evidence for repair involving gene conversion with various endogenous L1 elements, with some being used preferentially. In particular, the T(F) family and the L1Md-A2 subfamily, which are the most active in retrotransposition, appeared to be contributing the most in this process. The degree of homology did not seem to be a determining factor in the selection of the endogenous elements used for repair but may be based instead on accessibility. Considering their abundance and dispersion, gene conversion between repetitive elements may be occurring frequently enough to be playing a role in their evolution.     

Gene Conversion Among Paralogs Results in Moderate False Detection of Positive Selection Using Likelihood Methods

Previous studies have shown that recombination between allelic sequences can cause likelihood-based methods for detecting positive selection to produce many false-positive results. In this article, we use simulations to study the impact of nonallelic gene conversion on the specificity of PAML to detect positive selection among gene duplicates. Our results show that, as expected, gene conversion leads to higher rates of false-positive results, although only moderately. These rates increase with the genetic distance between sequences, the length of converted tracts, and when no outgroup sequences are included in the analysis. We also find that branch-site models will incorrectly identify unconverted sequences as the targets of positive selection when their close paralogs are converted. Bayesian prediction of sites undergoing adaptive evolution implemented in PAML is affected by conversion, albeit in a less straightforward way. Our work suggests that particular attention should be devoted to the evolutionary analysis of recent duplicates that may have experienced gene conversion because they may provide false signals of positive selection. Fortunately, these results also imply that those cases most susceptible to false-positive results—i.e., high divergence between paralogs, long conversion tracts—are also the cases where detecting gene conversion is the easiest. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Polymorphism and Gene Conversion in Mouse α-Globin Haplotypes

We have cloned and characterized three distinct alpha-globin haplotypes obtained from inbred strains of the mouse, Mus domesticus. We report here the complete nucleotide sequence of the six alpha-globin genes that the haplotypes contain. Our analysis of these genes and those from one other previously described haplotype indicates that recurrent gene conversion events have played a major role in their history. The pattern of nucleotide substitutions suggests that conversions have occurred both within and between haplotypes. Limited segments of coding and noncoding DNA have been involved in these gene conversion events. In two of the haplotypes, the nonallelic genes of each maintain DNA sequence identity over discrete intervals and encode the same alpha-globin polypeptide. On the other hand, the coding regions of some genes have accumulated replacement changes that result in distinct alpha-globins. In one instance, these changes appear to reflect positive selection of advantageous mutations.     

Duplicative and conservative transpositions of larval serum protein 1 genes in the genus Drosophila

Interspecific comparative molecular analyses of transposed genes and their flanking regions can help to elucidate the time, direction, and mechanism of gene transposition. In the Drosophila melanogaster genome, three Larval serum protein 1 (Lsp1) genes (alpha, beta and gamma) are present and each of them is located on a different chromosome, suggesting multiple transposition events. We have characterized the molecular organization of Lsp1 genes in D. buzzatii, a species of the Drosophila subgenus and in D. pseudoobscura, a species of the Sophophora subgenus. Our results show that only two Lsp1 genes (beta and gamma) exist in these two species. The same chromosomal localization and genomic organization, different from that of D. melanogaster, is found in both species for the Lsp1beta and Lsp1gamma genes. Overall, at least two duplicative and two conservative transpositions are necessary to explain the present chromosomal distribution of Lsp1 genes in the three Drosophila species. Clear evidence for implication of snRNA genes in the transposition of Lsp1beta in Drosophila has been found. We suggest that an ectopic exchange between highly similar snRNA sequences was responsible for the transposition of this gene. We have also identified the putative cis-acting regulatory regions of these genes, which seemingly transposed along with the coding sequences.     

Telomere conversion in trypanosomes.

Activation of the gene coding for variant surface glycoprotein (VSG) 118 in Trypanosoma brucei proceeds via a duplicative transposition to a telomeric expression site. The resulting active expression-linked extra copy (ELC) is usually flanked by DNA that lacks sites for most restriction enzymes and that is thought to interfere with the cloning of the ELC as recombinant DNA in Escherichia coli. We have circumvented this problem by cloning an aberrant 118 ELC gene, flanked at the 3'-side by at least 1 kb DNA, that contains restriction enzyme sites. Our analysis shows that this DNA and the 3'-end of the 118 ELC gene are derived from another VSG gene (1.1006) that is permanently located at a telomeric position. We propose that the 3'-end of the 1.1006 gene and (all of) its 3' flanking sequence moved to the expression site by a telomere conversion. Such a telomere conversion can also account for the appearance of an extra copy of the 1.1006 gene detected in a sub-population of our trypanosome strain.