Identification of Nucleotide Substitutions Necessary for Trans-Activation of Mariner Transposable Elements in Drosophila: Analysis of Naturally Occurring Elements

Six copies of the mariner element from the genomes of Drosophila mauritiana and Drosophila simulans were chosen at random for DNA sequencing and functional analysis and compared with the highly active element Mos1 and the inactive element peach. All elements were 1286 base pairs in length, but among them there were 18 nucleotide differences. As assayed in Drosophila melanogaster, three of the elements were apparently nonfunctional, two were marginally functional, and one had moderate activity that could be greatly increased depending on the position of the element in the genome. Both molecular (site-directed mutagenesis) and evolutionary (cladistic analysis) techniques were used to analyze the functional effects of nucleotide substitutions. The nucleotide sequence of the element is the primary determinant of function, though the activity level of elements is profoundly influenced by position effects. Cladistic analysis of the sequences has identified a T----A transversion at position 1203 (resulting in a Phe----Leu amino acid replacement in the putative transposase) as being primarily responsible for the low activity of the barely functional elements. Use of the sequences from the more distantly related species, Drosophila yakuba and Drosophila teissieri, as outside reference species, indicates that functional mariner elements are ancestral and argues against their origination by a novel mutation or by recombination among nonfunctional elements.     

Pseudogene evolution and natural selection for a compact genome

Pseudogenes are nonfunctional copies of protein-coding genes that are presumed to evolve without selective constraints on their coding function. They are of considerable utility in evolutionary genetics because, in the absence of selection, different types of mutations in pseudogenes should have equal probabilities of fixation. This theoretical inference justifies the estimation of patterns of spontaneous mutation from the analysis of patterns of substitutions in pseudogenes. Although it is possible to test whether pseudogene sequences evolve without constraints for their protein-coding function, it is much more difficult to ascertain whether pseudogenes may affect fitness in ways unrelated to their nucleotide sequence. Consider the possibility that a pseudogene affects fitness merely by increasing genome size. If a larger genome is deleterious--for example, because of increased energetic costs associated with genome replication and maintenance--then deletions, which decrease the length of a pseudogene, should be selectively advantageous relative to insertions or nucleotide substitutions. In this article we examine the implications of selection for genome size relative to small (1-400 bp) deletions, in light of empirical evidence pertaining to the size distribution of deletions observed in Drosophila and mammalian pseudogenes. There is a large difference in the deletion spectra between these organisms. We argue that this difference cannot easily be attributed to selection for overall genome size, since the magnitude of selection is unlikely to be strong enough to significantly affect the probability of fixation of small deletions in Drosophila.     

Periodic extinctions of transposable elements in bacterial lineages: evidence from intragenomic variation in multiple genomes

Most previous work on the evolution of mobile DNA was limited by incomplete sequence information. Whole genome sequences allow us to overcome this limitation. I study the nucleotide diversity of prominent members of five insertion sequence families whose transposition activity is encoded by a single transposase gene. Eighteen among 376 completely sequenced bacterial genomes and plasmids carry between 3 and 20 copies of a given insertion sequence. I show that these copies generally show very low DNA divergence. Specifically, more than 68% of the transposase genes are identical within a genome. The average number of amino acid replacement substitutions at amino acid replacement sites is Ka = 0.013, that at silent sites is Ks = 0.1. This low intragenomic diversity stands in stark contrast to a much higher divergence of the same insertion sequences among distantly related genomes. Gene conversion among protein-coding genes is unlikely to account for this lack of diversity. The relation between transposition frequencies and silent substitution rates suggests that most insertion sequences in a typical genome are evolutionarily young and have been recently acquired. They may undergo periodic extinction in bacterial lineages. By implication, they are detrimental to their host in the long run. This is also suggested by the highly skewed and patchy distribution of insertion sequences among genomes. In sum, one can think of insertion sequences as slow-acting infectious diseases of cell lineages.     

Genomic heterogeneity of background substitutional patterns in Drosophila melanogaster

Mutation is the underlying force that provides the variation upon which evolutionary forces can act. It is important to understand how mutation rates vary within genomes and how the probabilities of fixation of new mutations vary as well. If substitutional processes across the genome are heterogeneous, then examining patterns of coding sequence evolution without taking these underlying variations into account may be misleading. Here we present the first rigorous test of substitution rate heterogeneity in the Drosophila melanogaster genome using almost 1500 nonfunctional fragments of the transposable element DNAREP1_DM. Not only do our analyses suggest that substitutional patterns in heterochromatic and euchromatic sequences are different, but also they provide support in favor of a recombination-associated substitutional bias toward G and C in this species. The magnitude of this bias is entirely sufficient to explain recombination-associated patterns of codon usage on the autosomes of the D. melanogaster genome. We also document a bias toward lower GC content in the pattern of small insertions and deletions (indels). In addition, the GC content of noncoding DNA in Drosophila is higher than would be predicted on the basis of the pattern of nucleotide substitutions and small indels. However, we argue that the fast turnover of noncoding sequences in Drosophila makes it difficult to assess the importance of the GC biases in nucleotide substitutions and small indels in shaping the base composition of noncoding sequences.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Rapid sequence turnover at an intergenic locus in Drosophila

Closely related species of Drosophila tend to have similar genome sizes. The strong imbalance in favor of small deletions relative to insertions implies that the unconstrained DNA in Drosophila is unlikely to be passively inherited from even closely related ancestors, and yet most DNA in Drosophila genomes is intergenic and potentially unconstrained. In an attempt to investigate the maintenance of this intergenic DNA, we studied the evolution of an intergenic locus on the fourth chromosome of the Drosophila melanogaster genome. This 1.2-kb locus is marked by two distinct, large insertion events: a nuclear transposition of a mitochondrial sequence and a transposition of a nonautonomous DNA transposon DNAREP1_DM. Because we could trace the evolutionary histories of these sequences, we were able to reconstruct the length evolution of this region in some detail. We sequenced this locus in all four species of the D. melanogaster species complex: D. melanogaster, D. simulans, D. sechellia, and D. mauritiana. Although this locus is similar in size in these four species, less than 10% of the sequence from the most recent common ancestor remains in D. melanogaster and all of its sister species. This region appears to have increased in size through several distinct insertions in the ancestor of the D. melanogaster species complex and has been shrinking since the split of these lineages. In addition, we found no evidence suggesting that the size of this locus has been maintained over evolutionary time; these results are consistent with the model of a dynamic equilibrium between persistent DNA loss through small deletions and more sporadic DNA gain through less frequent but longer insertions. The apparent stability of genome size in Drosophila may belie very rapid sequence turnover at intergenic loci.     

Long interspersed repeated DNA (LINE) causes polymorphism at the rat insulin 1 locus.

The insulin 1, but not the insulin 2, locus is polymorphic (i.e., exhibits allelic variation) in rats. Restriction enzyme analysis and hybridization studies showed that the polymorphic region is 2.2 kilobases upstream of the insulin 1 coding region and is due to the presence or absence of an approximately 2.7-kilobase repeated DNA element. DNA sequence determination showed that this DNA element is a member of a long interspersed repeated DNA family (LINE) that is highly repeated (greater than 50,000 copies) and highly transcribed in the rat. Although the presence or absence of LINE sequences at the insulin 1 locus occurs in both the homozygous and heterozygous states, LINE-containing insulin 1 alleles are more prevalent in the rat population than are alleles without LINEs. Restriction enzyme analysis of the LINE-containing alleles indicated that at least two versions of the LINE sequence may be present at the insulin 1 locus in different rats. Either repeated transposition of LINE sequences or gene conversion between the resident insulin 1 LINE and other sequences in the genome are possible explanations for this.     

Size matters: non-LTR retrotransposable elements and ectopic recombination in Drosophila

The Drosophila melanogaster genome contains approximately 100 distinct families of transposable elements (TEs). In the euchromatic part of the genome, each family is present in a small number of copies (5-150 copies), with individual copies of TEs often present at very low frequencies in populations. This pattern is likely to reflect a balance between the inflow of TEs by transposition and the removal of TEs by natural selection. The nature of natural selection acting against TEs remains controversial. We provide evidence that selection against chromosome abnormalities caused by ectopic recombination limits the spread of some TEs. We also demonstrate for the first time that some TE families in the Drosophila euchromatin appear to be only marginally affected by purifying selection and contain many copies at high population frequencies. We argue that TEs in these families attain high population frequencies and even reach fixation as a result of low family-wide transposition rates leading to low TE copy numbers and consequently reduced strength of selection acting on individual TE copies. Fixation of TEs in these families should provide an upward pressure on the size of intergenic sequences counterbalancing rapid DNA loss through small deletions. Copy-number-dependent selection on TE families caused by ectopic recombination may also promote diversity among TEs in the Drosophila genome.     

Human-specific subfamilies of HERV-K (HML-2) long terminal repeats: three master genes were active simultaneously during branching of hominoid lineages

Using 40 known human-specific LTR sequences, we have derived a consensus sequence for an evolutionary young HERV-K (HML-2) LTR family, which was named the HS family. In the human genome the HS family is represented by approximately 150-160 LTR sequences, 90% of them being human-specific (hs). The family can be subdivided into two subfamilies differing in five linked nucleotide substitutions: HS-a and HS-b of 5.8 and 10.3 Myr evolutionary ages, respectively. The HS-b subfamily members were transpositionally active both before the divergence of the human and chimpanzee ancestor lineages and after it in both lineages. The HS-a subfamily comprises only hs LTRs. These and other data strongly suggest that at least three "master genes" of HERV-K (HML-2) LTRs were active in the human ancestor lineage after the human-chimpanzee divergence. We also found hs HERV-K (HML-2) LTRs integrations in introns of 12 human genes and identified 13 new hs HERV-K (HML-2) LTRs.     

Gradual evolution of a specific satellite DNA family in Drosophila ambigua, D. tristis, and D. obscura

The highly repetitive satellite DNA family "ATOC180" is specific for the three closely related species Drosophila obscura, D. ambigua, and D. tristis but does not occur in their closest relatives D. subsilvestris and D. bifasciata. Approximately 10,000 copies/haploid genome of approximately 180-bp repetition units are tandemly arranged in the centromeric heterochromatin of all chromosomes of all three species. Molecular analysis of 29 cloned repeats shows much intra- and interspecific sequence homogeneity. Single nucleotide changes are the main source of variability and distinguish the sequence-, subfamily- and species-specific ATOC180 repeats from each other. Based on these nucleotide differences, phylogenetic dendrograms were constructed and compared with published trees for other traits. The data indicate that the sequences of the ATOC180 satellite DNA family probably arose in a phylogenetically "short period" during the anagenetic evolution of the common ancestor of D. obscura, D. tristis, and D. ambigua, as a consequence of a process of genome reorganization, followed by a "long period" of entirely gradual sequence evolution. For the latter period, an evolutionary rate of 3 x 10(-8) substitutions/site/year was calculated.      

Horizontal transmission, vertical inactivation, and stochastic loss of mariner-like transposable elements

Horizontal transmission has been well documented as a major mechanism for the dissemination of mariner-like elements (MLEs) among species. Less well understood are mechanisms that limit vertical transmission of MLEs resulting in the "spotty" or discontinuous distribution observed in closely related species. In this article we present evidence that the genome of the common ancestor of the melanogaster species subgroup of Drosophila contained an MLE related to the mellifera (honey bee) subfamily. Horizontal transmission, approximately 3-10 MYA, is strongly suggested by the observation that the sequence of the MLE in Drosophila erecta is 97% identical in nucleotide sequence with that of an MLE in the cat flea, Ctenocephalides felis. The D. erecta MLE has a spotty distribution among species in the melanogaster subgroup. The element has a high copy number in D. erecta and D. orena, a moderate copy number in D. teissieri and D. yakuba, and was apparently lost ("stochastic loss") in the lineage leading to D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In D. erecta, most copies are concentrated in the heterochromatin. Two copies from D. erecta, denoted De12 and De19, were cloned and sequenced, and they appear to be nonfunctional ("vertical inactivation"). It therefore appears that the predominant mode of MLE evolution is vertical inactivation and stochastic loss balanced against occasional reinvasion of lineages by horizontal transmission.      

Factors affecting transposition of the Himar1 mariner transposon in vitro.

Mariner family transposable elements are widespread in animals, but their regulation is poorly understood, partly because only two are known to be functional. These are particular copies of the Dmmar1 element from Drosophila mauritiana, for example, Mos1, and the consensus sequence of the Himar1 element from the horn fly, Haematobia irritans. An in vitro transposition system was refined to investigate several parameters that influence the transposition of Himar1. Transposition products accumulated linearly over a period of 6 hr. Transposition frequency increased with temperature and was dependent on Mg2+ concentration. Transposition frequency peaked over a narrow range of transposase concentration. The decline at higher concentrations, a phenomenon observed in vivo with Mos1, supports the suggestion that mariners may be regulated in part by "overproduction inhibition." Transposition frequency decreased exponentially with increasing transposon size and was affected by the sequence of the flanking DNA of the donor site. A noticeable bias in target site usage suggests a preference for insertion into bent or bendable DNA sequences rather than any specific nucleotide sequences beyond the TA target site.     

A cloned I-factor is fully functional in Drosophila melanogaster

Summary I-R hybrid dysgenesis in Drosophila melanogaster occurs in female progeny of crosses between reactive strain females and inducer strain males, and is controlled by transposable elements called I-factors. These are 5.3 kb elements that are structurally similar to mammalian LINE elements and other retroposons. We have tested the activity of an I-factor directly, by introducing it into the genome of a reactive strain, using P-element mediated transformation. It confers the complete inducer phenotype on the reactive strain, and can stimulate dysgenesis when transformed males are mated with reactive females. It has transposed in the transformed lines, and we have cloned one of the transposed copies. This is the first time that it has been possible to demonstrate that a particular retroposon is transposition proficient, and to compare donor and transposed elements. We propose a mechanism for I-factor transposition based on these results, and the coding capacity of these elements. We have been unable to detect either autonomous transposition of a complete I-factor from a plasmid injected into reactive strain embryos, or transposition of a marked I-factor when co-injected with a complete element.