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.     

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.     

A microsatellite-based multilocus phylogeny of the Drosophila melanogaster species complex

Uncovering the genealogy of closely related species remains a major challenge for phylogenetic reconstruction. It is unlikely that the phylogeny of a single gene will represent the phylogeny of a species as a whole [1], but DNA sequence data across a large number of loci can be combined in order to obtain a consensus tree [2]. Long sequences are needed, however, to minimize the effect of (infrequent) base substitutions, and sufficient individuals must be sequenced per species to account for intraspecific polymorphisms, an overwhelming task using current DNA sequencing technology. By contrast, microsatellites are easy to type [3], allowing the analysis of many loci in multiple individuals. Despite their successful use in mapping [4,5], behavioural ecology [6] and population genetics [7], their usefulness for the phylogenetic reconstruction of closely related taxa has never been demonstrated, even though microsatellites are often conserved across species [8-10]. One drawback to microsatellite use is their high mutation rate (10(-4)-10(-2)), combined with an incomplete understanding of their mutation patterns. Many microsatellites are available for Drosophila melanogaster, and they are distributed throughout the genome [11]. Most can be amplified in the D. melanogaster species complex [12,13] and have low mutation rates [14, 15]. We show that microsatellite-specific distance measurements [16] correlate with other multilocus distances, such as those obtained from DNA-DNA hybridization data. Thus microsatellites may provide an ideal tool for building multilocus phylogenies. Our phylogenetic reconstruction of the D. melanogaster complex provides strong evidence that D. sechellia arose first, followed by a split between D. simulans and D. mauritiana.     

Polymorphism and divergence at a Drosophila pseudogene locus

The larval cuticle protein (Lcp) cluster in Drosophila melanogaster contains four functional genes and a closely related pseudogene. A 630-bp fragment including the larval cuticle pseudogene locus (Lcp psi) was nucleotide sequenced in 10 strains of D. melanogaster and a 458-bp Lcp psi fragment from D. simulans was also sequenced. We used these data to test the hypotheses that the rates of synonymous and nonsynonymous substitution are equal, that the absolute levels of variation are higher than in functional genes, and that intraspecific polymorphism is correlated with interspecific divergence. As predicted, synonymous and nonsynonymous substitution rates were equivalent, and overall nucleotide divergence between D. melanogaster and D. simulans (Jukes-Cantor distance = 0.149 +/- 0.150) was extremely high. However, within-species DNA sequence comparisons at Lcp psi revealed lower levels of polymorphism (theta = 0.001 +/- 0.001) than at many functional loci in D. melanogaster. Using the HUDSON, KREITMAN, and AGUADE (HKA) test, we show that the level of polymorphism in Lcp psi within D. melanogaster is lower than expected given the amount of divergence between D. melanogaster and D. simulans when the pseudogene data are compared to the Adh 5' flanking region. Because the Lcp psi lies in a region of relatively infrequent recombination, we suggest that the low level of within-species polymorphism is the result of background selection.     

Regularities in the development of children''s causality beliefs about school performance across six sociocultural contexts

The yellow (y) gene maps near the telomere of the X chromosome in Drosophila melanogaster but not in D. subobscura. Thus the strong reduction in the recombination rate associated with telomeric regions is not expected in D. subobscura. To study the divergence of a gene whose recombination rate differs between two species, the y gene of D. subobscura was sequenced. Sequence comparison between D. melanogaster and D. subobscura revealed several elements conserved in noncoding regions that may correspond to putative cis-acting regulatory sequences. Divergence in the y gene coding region between D. subobscura and D. melanogaster was compared with that found in other genes sequenced in both species. Both, yellow and scute exhibit an unusually high number of synonymous substitutions per site (ps). Also for these genes, the extent of codon bias differs between both species, being much higher in D. subobscura than in D. melanogaster. This pattern of divergence is consistent with the hitchhiking and background selection models that predict an increase in the fixation rate of slightly deleterious mutations and a decrease in the rate of fixation of slightly advantageous mutations in regions with low recombination rates such as in the y-sc gene region of D. melanogaster.     

Evolutionary inferences from DNA variation at the 6-phosphogluconate dehydrogenase locus in natural populations of drosophila: selection and geographic differentiation

Several allozyme-coding genes in Drosophila melanogaster show patterns suggesting that polymorphisms at these loci are targets of balancing selection. An important question is whether these genes have similar distributions of underlying DNA sequence variation which would indicate similar evolutionary processes occurring in this class of loci. One such locus, 6-phosphogluconate dehydrogenase (Pgd), has previously been shown to exhibit clinal variation for Fast/Slow electromorph variation in the United States and Australia, unusually large electromorph frequency differences between the United States and Africa, and other patterns indicative of selection. We measured four-cutter DNA restriction site and allozyme variation at Pgd among 142 D. melanogaster X chromosomes collected from several geographic regions including North Carolina, California, and Zimbabwe (Africa). We also sequenced a representative sample of 13 D. melanogaster Pgd genes collected in North Carolina and a single copy of Pgd from the sibling species, Drosophila simulans. While some population genetic models predict excess DNA polymorphism in genes which are targets of balancing selection, the D. melanogaster samples from the United States had significantly reduced levels of DNA polymorphism and extraordinarily high levels of linkage disequilibrium, providing evidence of hitchhiking effects of advantageous mutants at Pgd or at linked sites. Therefore, while selection has probably influenced the distribution of DNA variation at Pgd, the precise nature of these selective events remains obscure. Since the Pgd region appears to have low rates of crossing over, the reduced level of variation at this locus supports the idea that recombination rates are important determinants of levels of DNA polymorphism in natural populations. Furthermore, while patterns of allozyme variation are very similar at Pgd and Adh, the DNA data show that the evolutionary histories of these genes are dramatically different. We observed extensive differences in the amount and distribution of variation in D. melanogaster Pgd samples from the United States and Zimbabwe which cannot be explained by differential selection on the Fast/Slow polymorphism in these two geographic regions. Thus, genetic drift among partially isolated populations has also been an important factor in determining the distribution of variation at Pgd in D. melanogaster. Finally, we assayed four-cutter variation at Pgd in a sample of 19 D. simulans X chromosomes and observed reduced levels of DNA variability and high levels of linkage disequilibrium. These patterns are consistent with predictions of some hitchhiking models.     

The mutation rate and the distribution of mutational effects of viability and fitness in Drosophila melanogaster

The empirical distributions of the average viability and fitness of mutation accumulation lines of Drosophila melanogaster were analyzed using minimum distance estimation. Data come from two different experimental designs where mutations were allowed to accumulate: 1) in copies of chromosome II protected from natural selection and recombination (viability: Mukai et al., 1972; Ohnishi, 1977; fitness: Houle et al., 1992), 2) in inbred lines derived from the same isogenic stock (viability: Fernández & López-Fanjul, 1996; fitness: this paper). Information from all data sets converged, indicating that the mutational rates were small, about 1% for viability and 3% for fitness. For both traits, the rate of mutational decline appears to be smaller than suggested by previous studies (about one-fifth of the latter), the average mutational effect was neither severe nor very slight, ranging from -0.1 to -0.3, and the distribution of mutant effects was, at most, slightly leptokurtic. Therefore, the mutational load in natural populations is one to two orders of magnitude smaller than previously thought (as based upon analyses conditional to estimates of the mutational decline of viability or fitness that appear to be biased upward). Over 95% of the mutational variance of each trait was contributed by non-slightly deleterious mutations (absolute homozygous effect larger than 0.03 or 0.1, depending on the data set considered) occurring at a rate not higher than 0.025 per haploid genome and generation. Our data suggest that most deleterious mutations affecting fitness act mainly through a single component-trait.     

Reduced natural selection associated with low recombination in Drosophila melanogaster

Synonymous codons are not used equally in many organisms, and the extent of codon bias varies among loci. Earlier studies have suggested that more highly expressed loci in Drosophila melanogaster are more biased, consistent with findings from several prokaryotes and unicellular eukaryotes that codon bias is partly due to natural selection for translational efficiency. We link this model of varying selection intensity to the population-genetics prediction that the effectiveness of natural selection is decreased under reduced recombination. In analyses of 385 D. melanogaster loci, we find that codon bias is reduced in regions of low recombination (i.e., near centromeres and telomeres and on the fourth chromosome). The effect does not appear to be a linear function of recombination rate; rather, it seems limited to regions with the very lowest levels of recombination. The large majority of the genome apparently experiences recombination at a sufficiently high rate for effective natural selection against suboptimal codons. These findings support models of the Hill-Robertson effect and genetic hitchhiking and are largely consistent with multiple reports of low levels of DNA sequence variation in regions of low recombination.     

Mutation rate of a microsatellite sequence in normal human fibroblasts

Dinucleotide repeats, because of their repetitive nature, are prone to frameshift mutations, most likely via a DNA-polymerase slippage mechanism. Mutation rates in microsatellite DNA sequences are high in mismatch repair-defective cells. In normal cells, only estimates of maximal rates of mutation in microsatellites have been possible previously, because of the low sensitivity of screening assays for mutations in endogenous sequences. We have measured the spontaneous mutation rate of a dinucleotide repeat in diploid human foreskin fibroblasts. In our system, the mutation target is a (CA)17 repeat contained within a stably integrated plasmid. The repeat disrupts the reading frame of a neomycin (neo) resistance gene within the plasmid. Cells containing frameshift mutations in the CA repeat that correct the reading frame of the neo gene are selected using the neo analogue G418. This system of measuring microsatellite mutation rates is highly sensitive, because there is a specific target within which mutations can be selected. Fluctuation analysis of cells containing the target DNA yielded mutation rates of <3.1 x 10(-8) to 44.8 x 10(-8) mutations/cell/generation. This is the first report of a direct measurement of a spontaneous mutation rate of a microsatellite sequence in normal human cells.     

Inferences on genome-wide deleterious mutation rates in inbred populations of Drosophila and mice

A theoretical analysis was carried out on the mutation load observed in long-maintained inbred lines from two experiments with Drosophila and mice. The rate of decline in fitness and its sampling distribution were predicted for both experiments using Monte Carlo simulation with a range of mutational parameters and models. The predicted rates of change in fitness were compared to the empirical observed rates, which were close to zero. The classical hypothesis of many deleterious mutations (about one event per genome per generation) of small effect (1-2%) resulting in a mutation pressure for fitness of about 1% per generation is incompatible with the data. Recent estimates suggesting an overall mutation pressure for fitness traits of about 0.1% are, however, compatible with the observed load.     

Diffuse luminance flicker increases retinal vessel diameter in humans

To determine whether male- or female-biased mutation rates have affected the molecular evolution of Drosophila melanogaster and D. simulans, we calculated the male-to-female ratio of germline cell divisions ([symbol: see text]) from germline generation data and the male-to-female ratio of mutation rate ([symbol: see text]) by comparing chromosomal levels of nucleotide divergence. We found that the ratio of germline cell divisions changes from indicating a weak female bias to indicating a weak male bias as the age of reproduction increases. The range of [symbol: see text] values that we observed, however, does not lead us to expect much, if any, difference in mutation rate between the sexes. Silent and intron nucleotide divergence were compared between nine loci on the X chromosome and nine loci on the second and third chromosomes. The average levels of nucleotide divergence were not significantly different across the chromosomes, although both silent and intron sites show a trend toward slightly more divergence on the X. These results indicate a lack of sex- or chromosome-biased molecular evolution in D. melanogaster and D. simulans.     

Alterations of the enteric nervous system in neonatal necrotizing enterocolitis revealed by whole-mount immunohistochemistry

We report the results of a comprehensive search of Drosophila melanogaster DNA sequences in GenBank for di-, tri-, and tetranucleotide repeats of more than four repeat units, and a DNA library screen for dinucleotide repeats. Dinucleotide repeats are more abundant (66%) than tri- (30%) or tetranucleotide (4%) repeats. We estimate that 1917 dinucleotide repeats with 10 or more repeat units are present in the euchromatic D. melanogaster genome and, on average, they occur once every 60 kb. Relative to many other animals, dinucleotide repeats in D. melanogaster are short. Tri- and tetranucleotide repeats have even fewer repeat units on average than dinucleotide repeats. Our WorldWide Web site (http://www.bio.cornell.edu/genetics/aquadro/+ ++aquadro.html) posts the complete list of 1298 microsatellites (> or = five repeat units) identified from the GenBank search. We also summarize assay conditions for 70 D. melanogaster microsatellites characterized in previous studies and an additional 56 newly characterized markers.     

High genomic deleterious mutation rates in hominids

It has been suggested that humans may suffer a high genomic deleterious mutation rate. Here we test this hypothesis by applying a variant of a molecular approach to estimate the deleterious mutation rate in hominids from the level of selective constraint in DNA sequences. Under conservative assumptions, we estimate that an average of 4.2 amino-acid-altering mutations per diploid per generation have occurred in the human lineage since humans separated from chimpanzees. Of these mutations, we estimate that at least 38% have been eliminated by natural selection, indicating that there have been more than 1.6 new deleterious mutations per diploid genome per generation. Thus, the deleterious mutation rate specific to protein-coding sequences alone is close to the upper limit tolerable by a species such as humans that has a low reproductive rate, indicating that the effects of deleterious mutations may have combined synergistically. Furthermore, the level of selective constraint in hominid protein-coding sequences is atypically low. A large number of slightly deleterious mutations may therefore have become fixed in hominid lineages.     

Apoptotic repair of genotoxic tissue damage and the role of p53 gene

The spontaneous mutation rate per replication per genome is nearly invariant in microbes; however, the rate of spontaneous genomic mutations in higher eukaryotes is much higher. Furthermore, the mutation rates per locus per generation among Drosophila, mice and humans are similar, despite the large differences in generation time. A simple explanation for these findings is that mice and humans have a specific antimutagenic mechanism that is lacking in Drosophila. I propose that apoptotic repair-deletion of genotoxic damage-bearing cells-operates in mammalian germ cells and that it works more accurately in humans than in mice because of a slower rate of cell turnover and a longer generation time. It has been a long-standing puzzle that germline mutation frequencies decrease markedly as the dose-rate of radiation is lowered in mice but not in Drosophila. This can be readily explained by p53-dependent apoptotic repair, because the p53 gene is absent from the genome of Drosophila. Fetuses of p53+/+ mice have proficient apoptotic repair capacity for X-ray-induced teratogenic damage, but p53-null fetuses completely lack this capacity. Further, I propose that the primary role of the p53 gene is to guard germ cells and embryos from genotoxic damage. This implies that the tumour suppressor function of the p53 gene results from p53-dependent apoptotic deletion of cells with genotoxic damage. The reasoning behind this proposal is given by reviewing reports that Drosophila larvae are insensitive to tumour formation after irradiation. Finally, I discuss the genetic effects of radiation in humans.     

Deleterious mutation accumulation in organelle genomes

It is well established on theoretical grounds that the accumulation of mildly deleterious mutations in nonrecombining genomes is a major extinction risk in obligately asexual populations. Sexual populations can also incur mutational deterioration in genomic regions that experience little or no recombination, i.e., autosomal regions near centromeres, Y chromosomes, and organelle genomes. Our results suggest, for a wide array of genes (transfer RNAs, ribosomal RNAs, and proteins) in a diverse collection of species (animals, plants, and fungi), an almost universal increase in the fixation probabilities of mildly deleterious mutations arising in mitochondrial and chloroplast genomes relative to those arising in the recombining nuclear genome. This enhanced width of the selective sieve in organelle genomes does not appear to be a consequence of relaxed selection, but can be explained by the decline in the efficiency of selection that results from the reduction of effective population size induced by uniparental inheritance. Because of the very low mutation rates of organelle genomes (on the order of 10(-4) per genome per year), the reduction in fitness resulting from mutation accumulation in such genomes is a very long-term process, not likely to imperil many species on time scales of less than a million years, but perhaps playing some role in phylogenetic lineage sorting on time scales of 10 to 100 million years.     

A genetic system to detect mitotic recombination between repeated chromosomal sequences in Drosophila Schneider line 2 cells

In order to study mitotic homologous recombination in somatic Drosophila melanogaster cells in vitro and to learn more on the question how recombination is influenced by mutagens, a genetic system was developed where spontaneous and drug-induced recombination could be monitored. Two recombination reporter substrates were stably introduced in multiple copies into the genome of established D. melanogaster Schneider line 2 cells: one plasmid (pSB310) contained the 5' and 3' deleted neomycin phosphoribosyltransferase alleles neoL and neoR as direct repeats; the other (pSB485) contained similar deletions (lacZL and lacZR) of the beta-galactosidase gene (lacZ). Restoration of a functional neo gene upon mitotic recombination between homologous sequences allowed direct selection for the event, whereas recombination in single cells harbouring the integrated lacZ-based reporter plasmid was detected by histochemical staining or flow cytometric analysis (FACS). The neo-based construct in the clonal transgenic cell line 44CD4 showed a spontaneous recombination frequency of 2.9 x 10(-4), whereas the 485AD1 cell line harbouring the lacZ-based construct exhibited a frequency of 2.8 x 10(-4). The alkylating agents EMS and MMS and the clastogen mitomycin C were able to induce recombination in the 485AD1 cell line in a dose-dependent manner. The results obtained from these studies suggest that the transgenic cell lines are potentially useful tools for identifying agents which stimulate direct repeat recombination in somatic Drosophila cells.     

Evolution of the mitochondrial ATPase 6 gene in Drosophila: unusually high level of polymorphism in D. melanogaster

We have determined 1990 bp mitochondrial DNA sequence which extends from 3' end of the cytochrome oxidase subunit I (COI) gene to 5' end of the COIII gene from two sibling species of Drosophila, D. simulans and D. mauritiana. Analyses of the sequences and part of the NADH dehydrogenase subunit 2 gene and the COI gene together with those from D. melanogaster and D. yakuba revealed that amino-acid substitution rate of the ATPase 6 gene seems to be higher in some strains of D. melanogaster than in the other species. High level of amino-acid polymorphism in this gene was observed in D. melanogaster. Synonymous substitution rate is relatively constant in all the genes examined, suggesting that mutation rate is not higher in the ATPase 6 gene of D. melanogaster. The amino-acid substitutions found specifically in D. melanogaster are at the sites which are not conserved among mammals, yeast and E. coli. These sites of the ATPase 6 gene might lose the selective constraint in D. melanogaster, and the amino-acid substitutions can be explained by neutral mutations and random genetic drift.     

Polygenic mutation in Drosophila melanogaster: genotype x environment interaction for spontaneous mutations affecting bristle number

A highly inbred line of Drosophila melanogaster was subdivided into replicate sublines that were subsequently maintained independently with 10 pairs of parents per generation. The parents were randomly sampled for 19 'unselected' sublines and artificially selected for high or low abdominal or sternopleural bristle number for 12 'selected' sublines (with 3 replicate selection lines/trait/direction of selection). Divergence in mean bristle number among the unselected sublines, and response of the selected sublines to selection, are attributable to the accumulation of new mutations affecting bristle number. The input of mutational variance per generation, VM, can be estimated from the magnitude of response or divergence, assuming neutrality of mutations affecting the bristle traits. We reared unselected lines at generations 222 and 224, and selected lines at generations 182-184 of mutation accumulation at each of three temperatures (18 degrees C, 25 degrees C, 28 degrees C), and estimated the mutational variance common to all environments and the mutational variance from genotype x environment interaction. For sternopleural bristle number, the mutational interaction variance was 26% of the mutational variance common to all temperatures, and the interaction variance was due to temperature x line interaction. For abdominal bristle number, the mutational interaction variance was 142% of the mutational variance common to all temperatures, and the interaction variance was due to interactions of temperature x line, sex x line, and temperature x sex x line. It is possible that segregating variation for bristle number is maintained partly by genotype x environment interaction, but information on the fitness profiles of mutations affecting bristle number in each environment will be necessary to evaluate this hypothesis quantitatively.     

Naturally occurring variation in copia expression is due to both element (cis) and host (trans) regulatory variation

Significant differences in levels of copia [Drosophila long terminal repeat (LTR) retrotransposon] expression exist among six species representing the Drosophila melanogaster species complex (D. melanogaster, Drosophila mauritiana, Drosophila simulans, Drosophila sechellia, Drosophila yakuba, and Drosophila erecta) and a more distantly related species (Drosophila willistoni). These differences in expression are correlated with major size variation mapping to putative regulatory regions of the copia 5' LTR and adjacent untranslated leader region (ULR). Sequence analysis indicates that these size variants were derived from a series of regional duplication events. The ability of the copia LTR-ULR size variants to drive expression of a bacterial chloramphenicol acetyltransferase reporter gene was tested in each of the seven species. The results indicate that both element-encoded (cis) and host-genome-encoded (trans) genetic differences are responsible for the variability in copia expression within and between Drosophila species. This finding indicates that models purporting to explain the dynamics and distribution of retrotransposons in natural populations must consider the potential impact of both element-encoded and host-genome-encoded regulatory variation to be valid. We propose that interelement selection among retrotransposons may provide a molecular drive mechanism for the evolution of eukaryotic enhancers which can be subsequently distributed throughout the genome by retrotransposition.     

Contrasting patterns of nucleotide sequence variation at the glucose dehydrogenase (Gld) locus in different populations of Drosophila melanogaster

We have analyzed nucleotide sequence variation at the Glucose dehydrogenase (Gld) locus from four populations of Drosophila melanogaster from four continents. All four population samples show a significant reduction in silent variation compared to the neutral expectation. The levels of silent variation across all four populations are consistent with the predictions of the background selection model; however, Zimbabwe has a remarkably low level of variation. In the face of dramatically reduced silent polymorphism, an amino acid variant, leading to the common allozyme polymorphism at Gld, remains in low to intermediate frequency in all non-African samples. In the Chinese population sample, the ratio of replacement to silent variation is significantly elevated compared to the neutral expectation. The difference in patterns of variation across these population samples suggests that selection on Gld (or the Gld region) has been different in the Chinese population than in the other three.