The effects of weak selection on neutral diversity at linked sites

The effects of selection on variability at linked sites have an important influence on levels and patterns of within-population variation across the genome. Most theoretical models of these effects have assumed that selection is sufficiently strong that allele frequency changes at the loci concerned are largely deterministic. These models have led to the conclusion that directional selection for selectively favorable mutations, or against recurrent deleterious mutations, reduces nucleotide site diversity at linked neutral sites. Recent work has shown, however, that fixations of weakly selected mutations, accompanied by significant stochastic changes in allele frequencies, can sometimes cause higher diversity at linked sites when compared with the effects of fixations of neutral mutations. This study extends this work by deriving approximate expressions for the mean conditional times to fixation and loss of mutations subject to selection, and analyzing the conditions under which selection increases rather than reduces these times. Simulations are used to examine the relations between diversity at a neutral site and the fixation and loss times of mutations at a linked site that is subject to selection. It is shown that the long-term level of neutral diversity can be increased over the purely neutral value by recurrent fixations and losses of linked, weakly selected dominant or partially dominant favorable mutations, or linked recessive or partially recessive deleterious mutations. The results are used to examine the conditions under which associative overdominance, as opposed to background selection, is likely to operate.     

Mechanisms regulating the copy numbers of six LTR retrotransposons in the genome of Drosophila melanogaster

There has been debate over the mechanisms that control the copy number of transposable elements in the genome of Drosophila melanogaster. Target sites in D. melanogaster populations are occupied at low frequencies, suggesting that there is some form of selection acting against transposable elements. Three main theories have been proposed to explain how selection acts against transposable elements: insertions of a copy of a transposable element are selected against; chromosomal rearrangements caused by ectopic exchange between element copies are selected against; or the process of transposition itself is selected against. The three theories give different predictions for the pattern of transposable element insertions in the chromosomes of D. melanogaster. We analysed the abundance of six LTR (long terminal repeat) retrotransposons on the X and fourth chromosomes of multiple strains of D. melanogaster, which we compare with the predictions of each theory. The data suggest that no one theory can account for the insertion patterns of all six retrotransposons. Comparing our results with earlier work using these transposable element families, we find a significant correlation between studies in the particular model of copy number regulation supported by the proportion of elements on the X for the different transposable element families. This suggests that different retrotransposon families are regulated by different mechanisms.     

Genealogies and weak purifying selection

The assumption that selection alters the genealogical tree of a sample of alleles from a population relative to the neutral expectation underlies several "tests of neutrality." Two recent papers have studied the effect of purifying selection; their suggestive but incomplete results indicate that, in the single site case, the shape of a gene genealogy for a locus may differ only from the neutral expectation. We verify this finding for weak selection using the "ancestral selection graph." We consider a wider range of models, including both a four-allele single-site model and an infinite-sites model. Our results confirm the previous claim for the symmetric-mutation single site model. We emphasize, however, that a neutral-seeming genealogy is consistent with detectable effects of selection on the distribution of allele frequences within the sample. With selection operating, the information about a sample cannot be reduced to the genealogy. As a result, a distinction needs to be made between the selected sites themselves, for which the genealogy offers insufficient information, and linked neutral variation. This distinction seems to have been overlooked in previous papers, yet it has significant implications for the interpretation of data on DNA sequence variation. In particular, it predicts that under purifying selection, the frequency spectrum of neutral mutations will not reflect the skew toward rare polymorphisms at replacement sites even if there is no recombination between them. We caution, however, that the effect of weak selection on the genealogy is specific to the model; a (more realistic) model of multiple linked sites could lead to a more distorted genealogy than is observed for a single site.     

What is the impact of transposable elements on host genome variability?

The spread of a transposable element family through a wild population may be of astonishing rapidity. At least three families of transposable genetic elements have recently invaded Drosophila melanogaster worldwide, including the P element. The mechanism has been a process of effectively replicative transposition, and, for the P element, has occurred notwithstanding the sterility induced by unrestricted movement. This element's invasion into D. melanogaster has been accompanied by the development of heterogeneity between P sequences, most of which now have internal deletions. Increasing evidence suggests that some deleted elements can repress P transposition, thereby protecting the host from the harmful effects of complete elements. Such repressing elements may rise to high frequencies in populations as a result of selection at the level of the host. We here investigate selective sweeps invoked by the spread of P sequences in D. melanogaster populations. Numerous high-frequency sites have been identified on the X chromosome, which differ in frequency between populations, and which are associated with repression of P-element transposition. Unexpectedly, sequences adjacent to high-frequency P-element sites do not show reduced levels of genetic diversity, and DNA variability is in linkage equilibrium with the presence or absence of a P element at the adjacent selected site. This might be explained by multiple insertions or through a selection for recombination analogous to that seen in 'hitchhiking'.     

Does the proposed DSE motif form the active center in the Hermes transposase?

Donor cleavage and strand transfer are two functions performed by transposases during transposition of class II transposable elements. Within transposable elements, the only active center described, to date, facilitating both functions, is the so-called DDE motif. A second motif, R-K-H/K-R-H/W-Y, is found in the site-specific recombinases of the tyrosine recombinase family. While present in many bacterial insertion sequences as well as in the eukaryotic family of mariner/Tc1 elements, the DDE motif was considered absent in other classes of eukaryotic class II elements such as P, and hAT and piggyBac. Based on sequence alignments of a hobo-like element from the nematode Caenorhabditis elegans, to a variety of other hAT transposases and several members of the mariner/Tc1 group, Bigot et al. [Gene 174 (1996) 265] proposed the presence of a DSE motif in hAT transposases. In the present study we tested if each of these three residues is required for transposition of the Hermes element, a member of the hAT family commonly used for insect transformation. While D402N and E572Q mutations lead to knock-out of Hermes function, mutations S535A and S535D did not affect transposition frequency or the choice of integration sites. These data give the first experimental support that D402 and E572 are indeed required for transposition of Hermes. Furthermore, this study indicates that the active center of the Hermes transposase differs from the proposed DSE motif. It remains to be shown if other residues also form the active site of this transposase.     

A P Element Containing Suppressor of Hairy-Wing Binding Regions Has Novel Properties for Mutagenesis in Drosophila Melanogaster

P elements are widely used as insertional mutagens to tag genes, facilitating molecular cloning and analyses. We modified a P element so that it carried two copies of the suppressor of Hairy-wing [su(Hw)] binding regions isolated from the gypsy transposable element. This transposon was mobilized, and the genetic consequences of its insertion were analyzed. Gene expression can be altered by the su(Hw) protein as a result of blocking the interaction between enhancer/silencer elements and their promoter. These effects can occur over long distances and are general. Therefore, a composite transposon (SUPor-P for suppressor-P element) combines the mutagenic efficacy of the gypsy element with the controllable transposition of P elements. We show that, compared to standard P elements, this composite transposon causes an expanded repertoire of mutations and produces alleles that are suppressed by su(Hw) mutations. The large number of heterochromatic insertions obtained is unusual compared to other insertional mutagenesis procedures, indicating that the SUPor-P transposon may be useful for studying the structural and functional properties of heterochromatin.     

Ancestral inference on gene trees under selection

The extent to which natural selection shapes diversity within populations is a key question for population genetics. Thus, there is considerable interest in quantifying the strength of selection. A full likelihood approach for inference about selection at a single site within an otherwise neutral fully linked sequence of sites is described here. A coalescent model of evolution is used to model the ancestry of a sample of DNA sequences which have the selected site segregating. The mutation model, for the selected and neutral sites, is the infinitely many-sites model where there is no back or parallel mutation at sites. A unique perfect phylogeny, a gene tree, can be constructed from the configuration of mutations on the sample sequences under this model of mutation. The approach is general and can be used for any bi-allelic selection scheme. Selection is incorporated through modelling the frequency of the selected and neutral allelic classes stochastically back in time, then using a subdivided population model considering the population frequencies through time as variable population sizes. An importance sampling algorithm is then used to explore over coalescent tree space consistent with the data. The method is applied to a simulated data set and the gene tree presented in Verrelli et al. (2002).     

The Sleeping Beauty transposable element: evolution, regulation and genetic applications

Members of the Tc1/mariner superfamily of transposable elements isolated from vertebrate species are inactive due to the accumulation of mutations. A representative of a subfamily of fish elements estimated to be last active > 10 million years ago has been reconstructed, and named Sleeping Beauty(SB). This element opened up new avenues for studies on DNA transposition in vertebrates, and for the development of transposon tools for genetic manipulation in important model species and in humans. Multiple transposase binding sites within the terminal inverted repeats, a transpositional enhancer sequence, unequal affinity of the transposase to the binding sites and the activity of the cellular HMGB1 protein all contribute to a highly regulated assembly of SB synaptic complexes, which is likely a requirement for the subsequent catalytic steps. Host proteins involved in double-strand DNA break repair are limiting factors of SB transposition in mammalian cells, underscoring evolutionary, structural and functional links between DNA transposition, retroviral integration and V(D)J recombination. SB catalyzes efficient cut-and-paste transposition in a wide range of vertebrate cells in tissue culture, and in somatic tissues as well as the germline of the mouse and zebrafish in vivo, indicating its usefulness as a vector for transgenesis and insertional mutagenesis.     

Common physical properties of DNA affecting target site selection of sleeping beauty and other Tc1/mariner transposable elements

Sleeping Beauty (SB) is the most active Tc1/mariner-type transposable element in vertebrates, and is therefore a valuable vector for transposon mutagenesis in vertebrate models and for human gene therapy. We have analyzed factors affecting target site selection of SB in mammalian cells, by generating transposition events from extrachromosomal plasmids to chromosomes. In contrast to the local hopping observed when transposition is induced from a chromosomal context, mapping of 138 unique SB insertions on human chromosomes showed a fairly random genomic distribution, and a 35% occurrence of transposition into genes. Inspection of the DNA flanking the sites of element integration revealed significant differences from random DNA in both primary sequence and physical properties. The consensus sequence of SB target sites was found to be a palindromic AT-repeat, ATATATAT, in which the central TA is the canonical target site. We found however, that target site selection is determined primarily on the level of DNA structure, and not by specific base-pair interactions. Computational analyses revealed that insertion sites tend to have a bendable structure and a palindromic pattern of potential hydrogen-bonding sites in the major groove of the DNA. These features appear conserved in the Tc1/mariner family of transposons and in other, distantly related elements that share a common catalytic domain of the transposase, and integrate fairly randomly. No similar target site preference was found for non-randomly integrating elements. Our results suggest common factors influencing target site selection of a wide range of transposable elements.     

Preferential Transposition of Drosophila P Elements to Nearby Chromosomal Sites

Two different schemes were used to demonstrate that Drosophila P elements preferentially transpose into genomic regions close to their starting sites. A starting element with weak rosy(+) marker gene expression was mobilized from its location in the subtelomeric region of the 1,300-kb Dp1187 minichromosome. Among progeny lines with altered rosy(+) expression, a much higher than expected frequency contained new insertions on Dp1187. Terminal deficiencies were also recovered frequently. In a second screen, a rosy(+)-marked element causing a lethal mutation of the cactus gene was mobilized in male and female germlines, and viable revertant chromosomes were recovered that still contained a rosy(+) gene due to an intrachromosomal transposition. New transpositions recovered using both methods were mapped between 0 and 128 kb from the starting site. Our results suggested that some mechanism elevates the frequency 43-67-fold with which a P element inserts near its starting site. Local transposition is likely to be useful for enhancing the rate of insertional mutation within predetermined regions of the genome.     

The effects of deleterious mutations on evolution at linked sites

The process of evolution at a given site in the genome can be influenced by the action of selection at other sites, especially when these are closely linked to it. Such selection reduces the effective population size experienced by the site in question (the Hill-Robertson effect), reducing the level of variability and the efficacy of selection. In particular, deleterious variants are continually being produced by mutation and then eliminated by selection at sites throughout the genome. The resulting reduction in variability at linked neutral or nearly neutral sites can be predicted from the theory of background selection, which assumes that deleterious mutations have such large effects that their behavior in the population is effectively deterministic. More weakly selected mutations can accumulate by Muller's ratchet after a shutdown of recombination, as in an evolving Y chromosome. Many functionally significant sites are probably so weakly selected that Hill-Robertson interference undermines the effective strength of selection upon them, when recombination is rare or absent. This leads to large departures from deterministic equilibrium and smaller effects on linked neutral sites than under background selection or Muller's ratchet. Evidence is discussed that is consistent with the action of these processes in shaping genome-wide patterns of variation and evolution.     

Transposable Element-Induced Response to Artificial Selection in Drosophila Melanogaster: Molecular Analysis of Selection Lines

Artificial selection lines for abdominal bristle score of Drosophila melanogaster established from P-M hybrid dysgenic crosses showed increases in selection response, heritability and phenotypic variance compared to similar lines started from nondysgenic crosses. To determine whether this increased genetic variance could be due to enhanced transposition of P elements following the dysgenic cross, the cytological locations (sites) of P elements were determined by in situ hybridization for the whole genome of samples of 20 individuals from the parental P strain, 20 individuals from each of the eight dysgenic selection lines, and ten individuals from each of the eight nondysgenic selection lines. Variation among and within the selection lines and the parental P strain in P element insertion sites was exceptionally high. A total of 601 sites were identified, but there was no difference in total number of sites per line, mean number of sites per individual, mean copy number per individual, or site frequency between dysgenic and nondysgenic selection lines, or between lines selected for high and low bristle score. Transposition following nondysgenic crosses may explain additional observations of accelerated selection responses in nondysgenic selection lines. It was not possible to deduce which, if any, of the several hundred insertions in the dysgenic selection lines were responsible for their extreme bristle phenotypes.     

The Effect of Deleterious Mutations on Neutral Molecular Variation

Selection against deleterious alleles maintained by mutation may cause a reduction in the amount of genetic variability at linked neutral sites. This is because a new neutral variant can only remain in a large population for a long period of time if it is maintained in gametes that are free of deleterious alleles, and hence are not destined for rapid elimination from the population by selection. Approximate formulas are derived for the reduction below classical neutral values resulting from such background selection against deleterious mutations, for the mean times to fixation and loss of new mutations, nucleotide site diversity, and number of segregating sites. These formulas apply to random-mating populations with no genetic recombination, and to populations reproducing exclusively asexually or by self-fertilization. For a given selection regime and mating system, the reduction is an exponential function of the total mutation rate to deleterious mutations for the section of the genome involved. Simulations show that the effect decreases rapidly with increasing recombination frequency or rate of outcrossing. The mean time to loss of new neutral mutations and the total number of segregating neutral sites are less sensitive to background selection than the other statistics, unless the population size is of the order of a hundred thousand or more. The stationary distribution of allele frequencies at the neutral sites is correspondingly skewed in favor of rare alleles, compared with the classical neutral result. Observed reductions in molecular variation in low recombination genomic regions of sufficiently large size, for instance in the centromere-proximal regions of Drosophila autosomes or in highly selfing plant populations, may be partly due to background selection against deleterious mutations.     

Effect of Divergence Time and Recombination Rate on Molecular Evolution of <Emphasis Type="Italic">Drosophila</Emphasis> INE-1 Transposable Elements and Other Candidates for Neutrally Evolving Sites

Interspecies divergence of orthologous transposable element remnants is often assumed to be simply due to genetic drift of neutral mutations that occurred after the divergence of the species. However, divergence may also be affected by other factors, such as variation in the mutation rate, ancestral polymorphisms, or selection. Here we attempt to determine the impact of these forces on divergence of three classes of sites that are often assumed to be selectively unconstrained (INE-1 TE remnants, sites within short introns, and fourfold degenerate sites) in two different pairwise comparisons of Drosophila (D. melanogaster vs. D. simulans and D. simulans vs. D. sechellia). We find that divergence of these three classes of sites is strongly influenced by the recombination environment in which they are located, and this is especially true for the closer D. simulans vs. D. sechellia comparison. We suggest that this is mainly a result of the contribution of ancestral polymorphisms in different recombination regions. We also find that intergenic INE-1 elements are significantly more diverged than intronic INE-1 in both pairwise comparisons, implying the presence of either negative selection or lower mutation rates in introns. Furthermore, we show that substitution rates in INE-1 elements are not associated with the length of the noncoding sequence in which they are located, suggesting that reduced divergence in long noncoding sequences is not due to reduced mutation rates in these regions. Finally, we show that GC content for each site within INE-1 sequences has evolved toward an equilibrium value (approximately 33%) since insertion.     

A coalescent model of background selection with recombination,demography and variation in selection coefficients

There is increasing evidence that background selection, the effects of the elimination of recurring deleterious mutations by natural selection on variability at linked sites, may be a major factor shaping genome-wide patterns of genetic diversity. To accurately quantify the importance of background selection, it is vital to have computationally efficient models that include essential biological features. To this end, a structured coalescent procedure is used to construct a model of background selection that takes into account the effects of recombination, recent changes in population size and variation in selection coefficients against deleterious mutations across sites. Furthermore, this model allows a flexible organization of selected and neutral sites in the region concerned, and has the ability to generate sequence variability at both selected and neutral sites, allowing the correlation between these two types of sites to be studied. The accuracy of the model is verified by checking against the results of forward simulations. These simulations also reveal several patterns of diversity that are in qualitative agreement with observations reported in recent studies of DNA sequence polymorphisms. These results suggest that the model should be useful for data analysis.     

Population dynamics of an Ac-like transposable element in self- and cross-pollinating arabidopsis.

Theoretical models predict that the mating system should be an important factor driving the dynamics of transposable elements in natural populations due to differences in selective pressure on both element and host. We used a PCR-based approach to examine the abundance and levels of insertion polymorphism of Ac-III, a recently identified Ac-like transposon family, in natural populations of the selfing plant Arabidopsis thaliana and its close outcrossing relative, Arabidopsis lyrata. Although several insertions appeared to be ancient and shared between species, there is strong evidence for recent activity of this element family in both species. Sequences of the regions flanking insertions indicate that all Ac-III transposons segregating in natural populations are in noncoding regions and provide no evidence for local transposition events. Transposon display analysis suggests the presence of slightly higher numbers of insertion sites per individual but fewer total polymorphic insertions in the self-pollinating A. thaliana than A. lyrata. Element insertions appear to be segregating at significantly lower frequencies in A. lyrata than A. thaliana, which is consistent with a reduction in transposition rate, reduction in effective population size, or reduced efficacy of natural selection against element insertions in selfing populations.     

Transposition without transposase: a spontaneous mutation in bacteria.

Transposition mutations are typically associated with the activities of transposable elements such as transposons and insertion sequences, whose mobility is dependent upon transposase enzymes that catalyze exchanges between element ends and target sites. We describe a single transposition event in which a block of donor sequence is inserted at a target site without the involvement of any known transposase or the ends of any known transposable element. We propose that this is a new type of spontaneous mutation which may be difficult to detect in standard mutant hunts but may be of evolutionary importance.     

Molecular characterization of a mouse genomic element mobilized by advanced glycation endproduct modified-DNA (AGE-DNA).

BACKGROUND: DNA modified by advanced glycation endproducts (AGEs) undergoes a high frequency of insertional mutagenesis. In mouse lymphoid cells, these mutations are due in part to the transposition of host genomic elements that contain a DNA region homologous to the Alu family of repetitive elements. One particular 853 bp insertion, designated INS-1, was identified previously as a DNA element common to plasmids recovered from multiple, independent lymphoid cell transfections. MATERIALS AND METHODS: To characterize the genomic origin of this element, we used a 281-bp region of non-Alu-containing INS-1 sequence, designated. CORE, as a probe in Southern hybridization and for screening a bacteriophage mouse genomic DNA library. The resultant clones were sequenced and localized within the mouse genome. RESULTS: Two distinct genomic clones of 15 kB and 17 kB in size were isolated. A 522-bp unique region common to INS-1 and corresponding to the CORE sequence was identified in each clone. In both cases, CORE was found to be surrounded by repetitive DNA sequences: a 339-bp MT repeat at the 5' end, and a 150-bp B1 repeat at the 3' end. The CORE sequence was localized to mouse chromosome 1. CONCLUSIONS: These studies revealed that the CORE region of INS is present in low copy number but is associated with known repetitive DNA elements. The presence of these repetitive elements may facilitate the transposition of CORE by recombination or other, more complex rearrangement events, and explain in part the origin of AGE-induced insertional mutations.