Long microsatellite alleles in Drosophila melanogaster have a downward mutation bias and short persistence times, which cause their genome-wide underrepresentation

Microsatellites are short tandemly repeated DNA sequence motifs that are highly variable in most organisms. In contrast to mammals, long microsatellites (>15 repeats) are extremely rare in the Drosophila melanogaster genome. To investigate this paucity of long microsatellites in Drosophila, we studied 19 loci with exceptionally long microsatellite alleles. Inter- and intraspecific analysis showed that long microsatellite alleles arose in D. melanogaster only very recently. This lack of old alleles with many repeats indicated that long microsatellite alleles have short persistence times. The size distribution of microsatellite mutations in mutation-accumulation lines suggests that long alleles have a mutation bias toward a reduction in the number of repeat units. This bias causes the short persistence times of long microsatellite alleles. We propose that species-specific, size-dependent mutation spectra of microsatellite alleles may provide a general mechanism to account for the observed differences in microsatellite length between species.     

Relative rates of insertion and deletion mutations in dinucleotide repeats of various lengths in mismatch repair proficient mouse and mismatch repair deficient human cells

Microsatellites are DNA elements composed of short tandem repeats of 1-5bp. These sequences are particularly prone to frameshift mutation by insertion-deletion loop formation during replication. The mismatch repair system is responsible for correcting these replication errors, and microsatellite mutation rates are significantly elevated in the absence of mismatch repair. We have investigated the effect of varying the number of repeats in a (CA)n microsatellite on mutation rates in cultured mammalian cells proficient or deficient in mismatch repair. We have also compared the relative rates of single-repeat insertions and deletions in these cells. Two plasmid vectors were constructed for each repeat unit number (n=8, 17, and 30), such that the microsatellites, placed upstream of a bacterial neomycin resistance gene (neo), disrupted the reading frame of the gene in the (-1) or (+1) direction. Plasmids were introduced separately into the cells, where they integrated into the cellular genome. Mutation rates were determined by selection of clones with frameshift mutations in the microsatellite that restored the reading frame of the neo gene. We found that mutation rates were significantly higher for (CA)17 and (CA)30 tracts than for (CA)8 tracts in both mismatch repair proficient (mouse) and deficient (human) cells. A mutational bias favoring insertions was generally observed. In both (CA)17 and (CA)30 tracts, single-repeat insertion rates were higher than single-repeat deletion rates with or without mismatch repair; deletions of multiple repeat units (> or =8bp) were observed in these tracts, where as deletions this large were not found in the (CA)8 tract. Single-repeat mutations of both types were made at similar rates in (CA)8 tracts in human mismatch repair deficient (MMR-) cells, but single-repeat insertion rates were higher than single-repeat deletion rates in mouse mismatch repair proficient (MMR+) cells. Results of these direct studies on microsatellite mutations in cultured cells should be useful for refinement of mathematical models for microsatellite evolution.     

The mutation rates of di-, tri- and tetranucleotide repeats in Drosophila melanogaster

In a recent study, we reported that the combined average mutation rate of 10 di-, 6 tri-, and 8 tetranucleotide repeats in Drosophila melanogaster was 6.3 x 10(-6) mutations per locus per generation, a rate substantially below that of microsatellite repeat units in mammals studied to date (range = 10(-2)-10(-5) per locus per generation). To obtain a more precise estimate of mutation rate for dinucleotide repeat motifs alone, we assayed 39 new dinucleotide repeat microsatellite loci in the mutation accumulation lines from our earlier study. Our estimate of mutation rate for a total of 49 dinucleotide repeats is 9.3 x 10(-6) per locus per generation, only slightly higher than the estimate from our earlier study. We also estimated the relative difference in microsatellite mutation rate among di-, tri-, and tetranucleotide repeats in the genome of D. melanogaster using a method based on population variation, and we found that tri- and tetranucleotide repeats mutate at rates 6.4 and 8.4 times slower than that of dinucleotide repeats, respectively. The slower mutation rates of tri- and tetranucleotide repeats appear to be associated with a relatively short repeat unit length of these repeat motifs in the genome of D. melanogaster. A positive correlation between repeat unit length and allelic variation suggests that mutation rate increases as the repeat unit lengths of microsatellites increase.      

An unusually low microsatellite mutation rate in Dictyostelium discoideum, an organism with unusually abundant microsatellites

The genome of the social amoeba Dictyostelium discoideum is known to have a very high density of microsatellite repeats, including thousands of triplet microsatellite repeats in coding regions that apparently code for long runs of single amino acids. We used a mutation accumulation study to see if unusually high microsatellite mutation rates contribute to this pattern. There was a modest bias toward mutations that increase repeat number, but because upward mutations were smaller than downward ones, this did not lead to a net average increase in size. Longer microsatellites had higher mutation rates than shorter ones, but did not show greater directional bias. The most striking finding is that the overall mutation rate is the lowest reported for microsatellites: approximately 1 x 10(-6) for 10 dinucleotide loci and 6 x 10(-6) for 52 trinucleotide loci (which were longer). High microsatellite mutation rates therefore do not explain the high incidence of microsatellites. The causal relation may in fact be reversed, with low mutation rates evolving to protect against deleterious fitness effects of mutation at the numerous microsatellites.     

High mutation rate of a long microsatellite allele in Drosophila melanogaster provides evidence for allele-specific mutation rates

Within recent years, microsatellite have become one of the most powerful genetic markers in biology. For several mammalian species, microsatellite mutation rates have been estimated on the order of 10(- 3)-10(-5). A recent study, however, demonstrated mutation rates in Drosophila melanogaster of at least one order of magnitude lower than those in mammals. To further test this result, we examined mutation rates of different microsatellite loci using a larger sample size. We screened 24 microsatellite loci in 119 D. melanogaster lines maintained for approximately 250 generations and detected 9 microsatellite mutations. The average mutation rate of 6.3 x 10(-6) is identical to the mutation rate from a previous study. Most interestingly, all nine mutations occurred at the same allele of one locus (DROYANETSB). This hypermutable allele has 28 dinucleotide repeats and is among the longest microsatellite reported in D. melanogaster. The allele-specific mutation rate of 3.0 x 10(-4) per generation is within the range of mammalian mutation rates. Future microsatellite analyses will have to account for the dramatic differences in allele-specific mutation rates.      

Mobile dispersed genetic element MDG1 of Drosophila melanogaster: structural organization.

The whole-length mobile dispersed genetic element mdg1 has been cloned from D. melanogaster genome. It contains DNA fragments described earlier as Dm225 and Dm234, Mdg1 is 7.2 kb long and framed with two direct repeats of 300-400 base pairs each. Mdg1 family is represented by about 25 copies in the genome of flies and by 200 copies in the genome of cultured cell line 67J25D. Virtually all the copies in the genome of D. melanogaster have the same restriction map. Oligo(dA)-oligo(dT) regions were found within mdg1.     

Drosophila virilis has long and highly polymorphic microsatellites

Comparative genomics is a powerful approach to inference of the dynamics of genome evolution. Most information about the evolution of microsatellites in the genus Drosophila has been obtained from Drosophila melanogaster. For comparison, we collected microsatellite data for the distantly related species Drosophila virilis. Screening about 0.5 Mb of nonredundant genomic sequence from GenBank, we identified 239 dinucleotide microsatellites. On average, D. virilis dinucleotides were significantly longer than D. melanogaster microsatellites (7.69 repeats vs. 6.75 repeats). Similarly, direct cloning of microsatellites resulted in a higher mean repeat number in D. virilis than in D. melanogaster (12.7 repeats vs. 12.2 repeats). Characterization of 11 microsatellite loci mapping to division 40-49 on the fourth chromosome of D. virilis indicated that D. virilis microsatellites are more variable than those of D. melanogaster.     

Drosophila bloom helicase maintains genome integrity by inhibiting recombination between divergent DNA sequences

DNA double strand breaks (DSB) can be repaired either via a sequence independent joining of DNA ends or via homologous recombination. We established a detection system in Drosophila melanogaster to investigate the impact of sequence constraints on the usage of the homology based DSB repair via single strand annealing (SSA), which leads to recombination between direct repeats with concomitant loss of one repeat copy. First of all, we find the SSA frequency to be inversely proportional to the spacer length between the repeats, for spacers up to 2.4 kb in length. We further show that SSA between divergent repeats (homeologous SSA) is suppressed in cell cultures and in vivo in a sensitive manner, recognizing sequence divergences smaller than 0.5%. Finally, we demonstrate that the suppression of homeologous SSA depends on the Bloom helicase (Blm), encoded by the Drosophila gene mus309. Suppression of homeologous recombination is a novel function of Blm in ensuring genomic integrity, not described to date in mammalian systems. Unexpectedly, distinct from its function in Saccharomyces cerevisiae, the mismatch repair factor Msh2 encoded by spel1 does not suppress homeologous SSA in Drosophila.     

Length of CTG.CAG repeats determines the influence of mismatch repair on genetic instability

We showed previously that mutations in methyl-directed mismatch repair of Escherichia coli reduced the occurrence of large deletions in (CTG.CAG)(175) repeats contained on plasmids. By contrast, other workers reported that mutations in mismatch repair increase the frequency of small-length changes in the shorter (CTG.CAG)(64). Using plasmids with a variety of lengths and purity of (CTG.CAG) repeats, we have resolved these apparently conflicting observations. We show that all lengths of (CTG.CAG) repeats are subject to small-length changes (eight repeats) in (CTG.CAG)(n) occur more readily in cells with active mismatch repair. The frequency of large deletions is proportional to the tract length; in our assays they become prominent in tracts greater than 100 repeats. Interruptions in repeat purity enhance the occurrence of large deletions. In addition, we observed a high level of incidence of deletions in (CTG.CAG) repeats for cultures passing repeatedly through stationary phase during long-term growth experiments of all strains (i.e. with active or inactive mismatch repair). These results agree with current theories on mismatch repair acting on DNA slippage events that occur in DNA triplet-repeats.     

Characterization and isolation of novel microsatellites from the Drosophila dunni subgroup

We have isolated and characterized 77 novel microsatellites from two species, Drosophila dunni and Drosophila nigrodunni, which are closely related Caribbean-island endemics from the Drosophila cardini species group. These species are very distantly related to all other Drosophila from which microsatellites have previously been characterized. We find that the average length of microsatellites isolated in these species is quite small, with an overall mean length of 9.8 repeat units for dinucleotide microsatellites in the two study species. The nucleotide composition of dinucleotides differs between the two species: D. nigrodunni has a predominance of (AC/GT)n repeats, whereas D. dunni has equal numbers of (AC/GT)n and (AG/CT)n repeats. Tri- and tetranucleotide repeats are not abundant in either species. We assayed the variability of eight microsatellites in a closely related third species, Drosophila arawakana, using wild-caught individuals from the island of Guadeloupe. We found the microsatellites to be extremely variable in this population, with observed heterozygosities ranging from 0.541 to 0.889. DNA amplification trials suggest that these eight microsatellites are widely conserved across the D. cardini group, with five of the eight producing amplification products in every species tested. However, the loci are very poorly conserved over greater phylogenetic distances. DNA amplification of the microsatellite loci was unreliable in members of the closely related Drosophila quinaria, Drosophila calloptera, Drosophila guarani and Drosophila tripunctata species groups. Furthermore, these microsatellites could not be detected in the genome of Drosophila melanogaster, despite the conservation of microsatellite flanking regions at some loci. These data indicate that Drosophila microsatellite loci are quite short lived over evolutionary timescales relative to many other taxa.     

Germ cells microsatellite instability. The effect of different mutagens in a mismatch repair mutant of Drosophila (spel1)

Mismatch repair (MMR) process confers a type of genomic stability that maintains stable single repeated sequences, hence a failure of this process could deviate in cancer development. A characteristic phenotype of MMR-deficient cells is microsatellite instability (MSI) that could be modulated by mutagenic agents. The induction of MSI by the mutagens, bleomycin (BLM), hydrogen peroxide (H(2)O(2)), 2-acetylaminofluorene (2-AAF) and ethidium bromide (EB) was evaluated in vivo, by using a Drosophila melanogaster-null mutant of the msh2 mismatch repair gene (spel1). Whereas in the germ cells of the spel1 strain, we found microsatellite mutations in the five repeated sequences studied in untreated individuals, no alterations were found in the MMR-proficient strain. On the other hand, the data obtained from the treatment experiments show that BLM and 2-AAF induced a slight mutagenic effect in the MMR-deficient background but not in the normal one. These results indicate that the use of the Drosophila spel1 mutant (MMR-deficient) could be of relevant importance to identify environmental factors involved in carcinogenesis processes through genomic instability.     

Analysis of genomic damage in the mutagen-sensitive mus-201 mutant of Drosophila melanogaster by arbitrarily primed PCR (AP-PCR) fingerprinting.

DNA repair mechanisms are important to maintain the stability of the genome. In Drosophila melanogaster, the mus-201 gene is required in the excision repair process. To study the contribution of the mus-201 gene in the stability of the Drosophila genome, we have used the arbitrarily primed PCR fingerprinting method (AP-PCR). We have analysed the changes in the genomic DNA fingerprints from the progeny of wild-type males crossed with mus-201 repair-deficient or repair-proficient females. After induction of DNA damage with 2-acetylaminofluorene (2-AAF) in the wild-type parental males, quantitative and qualitative differences in the AP-PCR fingerprints were detected between the two crosses, and the estimate of the genomic damage detected by AP-PCR has clearly shown that the mus-201 repair deficiency is associated with an increase of genomic damage. The predominant type of alterations detected by AP-PCR under the mus-201 repair-deficient conditions agree with the results obtained in microsatellite PCR analysis, suggesting that the role of the mus-201 gene, necessary in excision repair, is not associated to the mismatch repair process. The work reported here demonstrates that the AP-PCR is a suitable technique to analyse genetic alterations in D. melanogaster and, consequently, can be used to compare the susceptibility to genomic damage of different DNA repair mutants.     

High mutation rate of TPE repeats: a microsatellite in the putative transposase of the hobo element in Drosophila melanogaster

The hobo transposable element contains a polymorphic microsatellite sequence located in its coding region, the TPE repeats. Previous surveys of natural populations of Drosophila melanogaster have detected at least seven different hobo transposons. These natural populations are geographically structured with regard to TPE polymorphism, and a scenario has been proposed for the invasion process. Natural populations have recently been completely invaded by hobo elements with three TPE repeats. New elements then appeared by mutation, triggering a new stage of invasion by other elements. Since TPE polymorphism appeared over a short period of time, we focused on estimating the mutation rate of these TPE repeats. We used transgenic lines harboring three TPE and/or five TPE hobo elements that had been evolving for at least 16 generations to search for a new TPE repeat polymorphism. We detected three mutants, with four, seven, and eight TPE repeats, respectively. The estimated mutation rate of the TPE repeats is therefore higher than that of neutral microsatellites in D. melanogaster (4.2 x 10-4 versus 6.5 x 10-6). The role of the transposition mechanism and the particular structure of the TPE repeats of the hobo element in this increase in the mutation rate are discussed.     

Polarity of mutations in tumor-associated microsatellite instability

Many factors have been implicated in influencing the rate of microsatellite mutations, including the length and base composition of the repeat motif, number of repeats, base composition of flanking sequences and, perhaps most importantly, degree of perfection of the repeats. The latter is of clinical relevance, since in both spino-cerebellar ataxia and fragile X syndrome, alleles with imperfect repeats appear to be much more stable than perfect ones. As yet, the relative importance of increased replication slippage and decreased mismatch repair efficiency in the preference of mutations to occur within perfect repeats has not been fully determined. D13S308E is an asymmetric trinucleotide repeat microsatellite with the sequence (CAT)₃CAC(CAT)CAC(CAT)₂CAC(CAT)CAC(CAT) ₁₅ , thus containing two parts: an 11-repeat imperfect portion (underlined above) and a 15-repeat perfect one (bold). We sequenced eight new mutant alleles of D13S308E from three human gastric tumors with instability in this and other microsatellites. In all mutations the size variation occurred exclusively in the perfect part of the microsatellite. These results constitute direct evidence that the molecular basis of microsatellite alterations seen in normal cells is similar to those that occur in human tumors with extensive microsatellite instability. The investigation of mechanisms involved in microsatellite mutations has been handicapped by the fact that they are rare events. The microsatellite instability observed in malignant tumors provides us with a useful general system to study these mechanisms. Received: 24 June 1997 / Accepted 7 October 1997     

Distribution of dinucleotide microsatellites in the Drosophila melanogaster genome.

Microsatellites, a special class of repetitive DNA, have become one of the most popular genetic markers. The progress of various genome projects has made it possible to study the genomic distribution of microsatellites and to evaluate the potential influence of several parameters on their genesis. We report the distribution of dinucleotide microsatellites in the genome of Drosophila melanogaster. When considering only microsatellites with five or more repeat units, the average length of dinucleotide repeats in D. melanogaster is 6.7 repeats. We tested a wide range of parameters which could potentially influence microsatellite density, and we did not detect a significant influence of recombination rate, number of exons, or total length of coding sequence. In concordance with the neutral expectation for the origin of microsatellites, a significant positive correlation between AT content and (AT/TA)n microsatellite density was detected. While this pattern may indicate that microsatellite genesis is a random process, we also found evidence for a nonrandom distribution of microsatellites. Average microsatellite density was higher on the X chromosome, but extreme heterogeneity was observed between different genomic regions. Such a clumping of microsatellites was also evident on a more local scale, as 38.9% of the contiguous sequences analyzed showed a deviation from a random distribution of microsatellites.     

DNA mismatch repair,microsatellite instability and cancer

Mismatch (MMR) repair system plays a significant role in restoration of stability in the genome. Mutations in mismatch repair genes hamper their activity thus bring about a defect in mismatch repair (MMR) mechanism thereby conferring instability in the microsatellite sequences of both the coding and non-coding regions of the genome. Mutated mismatch repair genes result in the expansion or contraction of microsatellite sequence and confer microsatellite unstable or replication error positive phenotype. Hypermethylation of promoter regions of some of the MMR genes also causes inactivation of these genes and thus contribute to MSI. Microsatellite instability is an indicator of MMR deficiency and is a prime cause of varied tumorogenesis.     

Detection of mosaicism in lymphocytes of parents of free trisomy 21 offspring

Genetic selection assays were developed to measure rates of deletion of one or more (CAG).(CTG) repeats, or an entire repeat tract, in Escherichia coli. In-frame insertions of >or=25 repeats in the chloramphenicol acetyltransferase (CAT) gene of pBR325 resulted in a chloramphenicol-sensitive (Cm(s)) phenotype. When (CAG)25 comprised the leading template strand, deletion of one or more repeats resulted in a chloramphenicol resistant (Cm(r)) phenotype at a rate of 4 x 10(-2) revertants per cell per generation. The mutation rates for plasmids containing (CAG)43 or (CAG)79 decreased significantly. When (CTG)n comprised the leading template strand the Cm(r) mutation rates were 100-1000 lower than for the opposite orientation. As an initial application of this assay, the effects of mutations influencing mismatch repair and recombination were examined. The methyl directed mismatch repair system increased repeat stability only when (CTG)n comprised the leading template strand. Replication errors made with the opposite repeat orientation were apparently not recognized. For the (CAG)n leading strand orientation, mutation rates were reduced as much as 3000-fold in a recA- strain. In a second assay, out-of-frame mutation inserts underwent complete deletion at rates ranging from about 5 x 10(-9) to 1 x 10(-7) per cell per generation. These assays allow careful quantitation of triplet repeat instability in E. coli and provide a way to examine the effects of mutations in replication, repair, and recombination on repeat instability.     

High density of long dinucleotide microsatellites in Drosophila subobscura

We isolated 96 dinucleotide repeats with five or more tandemly repeated units from a subgenomic Drosophila subobscura library. The mean repeat unit length of microsatellite clones in D. subobscura is 15, higher than that observed in other Drosophila species. Population variation was assayed in 32-40 chromosomes from Barcelona, Spain, using 18 randomly chosen microsatellite loci. Positive correlation between measures of variation and perfect repeat length measures (mean size, most common, and longest allele) is consistent with a higher mutation rate in loci with longer repeat units. Levels of microsatellite variation measured as variance in repeat number and heterozygosity in D. subobscura were similar to those of Drosophila pseudoobscura and higher than those of Drosophila melanogaster and Drosophila simulans. Our data suggest that higher levels of microsatellite variation, and possibly density, in D. subobscura compared with D. melanogaster are due to both a higher average effective population and a higher intrinsic slippage rate in the former species.     

Integrating genomics,bioinformatics, and classical genetics to study the effects of recombination on genome evolution

This study presents compelling evidence that recombination significantly increases the silent GC content of a genome in a selectively neutral manner, resulting in a highly significant positive correlation between recombination and "GC3s" in the yeast Saccharomyces cerevisiae. Neither selection nor mutation can explain this relationship. A highly significant GC-biased mismatch repair system is documented for the first time in any member of the Kingdom Fungi. Much of the variation in the GC3s within yeast appears to result from GC-biased gene conversion. Evidence suggests that GC-biased mismatch repair exists in numerous organisms spanning six kingdoms. This transkingdom GC mismatch repair bias may have evolved in response to a ubiquitous AT mutational bias. A significant positive correlation between recombination and GC content is found in many of these same organisms, suggesting that the processes influencing the evolution of the yeast genome may be a general phenomenon. Nonrecombining regions of the genome and nonrecombining genomes would not be subject to this type of molecular drive. It is suggested that the low GC content characteristic of many nonrecombining genomes may be the result of three processes (1) a prevailing AT mutational bias, (2) random fixation of the most common types of mutation, and (3) the absence of the GC-biased gene conversion which, in recombining organisms, permits the reversal of the most common types of mutation. A model is proposed to explain the observation that introns, intergenic regions, and pseudogenes typically have lower GC content than the silent sites of corresponding open reading frames. This model is based on the observation that the greater the heterology between two sequences, the less likely it is that recombination will occur between them. According to this "Constraint" hypothesis, the formation and propagation of heteroduplex DNA is expected to occur, on average, more frequently within conserved coding and regulatory regions of the genome. In organisms possessing GC-biased mismatch repair, this would enhance the GC content of these regions through biased gene conversion. These findings have a number of important implications for the way we view genome evolution and suggest a new model for the evolution of sex.