Microsatellite Allelic Homoplasy Due to Variable Flanking Sequences

Microsatellite DNA sequences have become the dominant source of nuclear genetic markers for most applications. It is important to investigate the basis of variation between alleles and to know if current assumptions about the mechanisms of microsatellite mutation (that is to say, variations involving simple changes in the number of repeat) are correct. We have characterized, by DNA sequencing, the human alleles of a new highly informative (CA)n repeat localized approximately 20 kb centromeric to the HLA-B gene. Although 12 alleles were identified based on conventional length criteria, sequencing of the alleles demonstrated that differences between alleles were found to be more complex than previously assumed: A high degree of microsatellite variability is due to variation in the region immediately flanking the repeat. These data indicate that the mutational process which generates polymorphism in this region has involved not only simple changes in the number of dinucleotide CA repeats but also perturbations in the nonrepeated 5′ and 3′ flanking sequences. Three families of alleles (not visible from the overall length of the alleles), with presumably separate evolutionary histories, exist and can yield to homoplasy of size. Effectively, we can observe alleles of the same size with different internal structures which are separated by a significant amount of variation. Although allelic homoplasy for noninterrupted microsatellite loci has been suggested between different species, it has not been unequivocally demonstrated within species. A strong association is noted between alleles defined at the sequence level and HLA-B alleles. The observation of several families of alleles at the population level provides information about the evolutionary history and mutation processes of microsatellites and may have implications for the use of these markers in phylogenetic, linkage disequilibrium studies, and gene mapping. Received: 14 May 1996 / Accepted: 9 September 1996     

Microsatellites in Zea - variability,patterns of mutations,and use for evolutionary studies

To evaluate the performance of microsatellites or simple sequence repeats (SSRs) for evolutionary studies in Zea, 46 microsatellite loci originally derived from maize were applied to diverse arrays of populations that represent all the diploid species of Zea and 101 maize inbreds. Although null phenotypes and amplification of more than two alleles per plant were observed at modest rates, no practical obstacle was encountered for applying maize microsatellites to other Zea species. Sequencing of microsatellite alleles revealed complex patterns of mutation including frequent indels in the regions flanking microsatellite repeats. In one case, all variation at a microsatellite locus came from indels in the flanking region rather than in the repeat motif. Maize microsatellites show great variability within populations and provide a reliable means to measure intraspecific variation. Phylogeographic relationships of Zea populations were successfully reconstructed with good resolution using a genetic distance based on the infinite allele model, indicating that microsatellite loci are useful in evolutionary studies in Zea. Microsatellite loci show a principal division between tropical and temperate inbred lines, and group inbreds within these two broad germplasm groups in a manner that is largely consistent with their known pedigrees. Received: 10 February 2001 / Accepted: 21 May 2001     

Phylogenetic Analysis of the Friedreich Ataxia GAA Trinucleotide Repeat

Friedreich ataxia is an autosomal recessive neurodegenerative disorder associated with a GAA repeat expansion in the first intron of the gene (FRDA) encoding a novel, highly conserved, 210 amino acid protein known as frataxin. Normal variation in repeat size was determined by analysis of more than 600 DNA samples from seven human populations. This analysis showed that the most frequent allele had nine GAA repeats, and no alleles with fewer than five GAA repeats were found. The European and Syrian populations had the highest percentage of alleles with 10 or more GAA repeats, while the Papua New Guinea population did not have any alleles carrying more than 10 GAA repeats. The distributions of repeat sizes in the European, Syrian, and African American populations were significantly different from those in the Asian and Papua New Guinea populations (p < 0.001). The GAA repeat size was also determined in five nonhuman primates. Samples from 10 chimpanzees, 3 orangutans, 1 gorilla, 1 rhesus macaque, 1 mangabey, and 1 tamarin were analyzed. Among those primates belonging to the Pongidae family, the chimpanzees were found to carry three or four GAA repeats, the orangutans had four or five GAA repeats, and the gorilla carried three GAA repeats. In primates belonging to the Cercopithecidae family, three GAA repeats were found in the mangabey and two in the rhesus macaque. However, an AluY subfamily member inserted in the poly(A) tract preceding the GAA repeat region in the rhesus macaque, making the amplified sequence approximately 300 bp longer. The GAA repeat was also found in the tamarin, suggesting that it arose at least 40 million years ago and remained relatively small throughout the majority of primate evolution, with a punctuated expansion in the human genome. Received: 18 August 2000 / Accepted: 10 November 2000     

Microsatellite Evolution: Testing the Ascertainment Bias Hypothesis

Previous studies suggest the median allele length of microsatellites is longest in the species from which the markers were derived, suggesting that an ascertainment bias was operating. We have examined whether the size distribution of microsatellite alleles between sheep and cattle is source dependent using a set of 472 microsatellites that can be amplified in both species. For those markers that were polymorphic in both species we report a significantly greater number of markers (P < 0.001) with longer median allele sizes in sheep, regardless of microsatellite origin. This finding suggests that any ascertainment bias operating during microsatellite selection is only a minor contributor to the variation observed. Received: 6 January 1997 / Accepted: 19 May 1997     

Sequence diversity and molecular evolution of the merozoite surface antigen 2 of plasmodium falciparum

Eleven new alleles of the Plasmodium falciparum merozoite surface antigen 2 (MSA2) from Papua New Guinea were analyzed by direct sequencing of polymerase chain reaction (PCR) products. We have used the sequence information to trace the molecular evolution of MSA2. The repeats of ten alleles belonging to the 3D7 allelic family differed considerably in size, nucleotide sequence, and repeat copy number. In the repeat region of these new alleles, codon usage was extremely biased with an exclusive use of NNT codons. Another new allele sequenced belonged to the FC27 family and confirmed the family-specific conserved structure of 96 and 36 bp repeats. In order to assess sequence microheterogeneity within samples defined as the same genotype by restriction fragment length polymorphism (RFLP), we have analyzed single-strand conformation polymorphism (SSCP) of different samples of the most frequent allele (D10 of the FC27 family) in the study population. No sequence heterogeneity could be detected within the repeat region. Based on analysis of the repeat regions in both allelic families, we discuss the hypothesis of a different evolutionary strategy being represented by each of the allelic families. Received: 8 February 1995 / Accepted: 24 March 1997     

Length polymorphism and allele structure of trinucleotide microsatellites in natural accessions of Arabidopsis thaliana

 The objective of this work was to assess the degree of trinucleotide microsatellite length polymorphism in the selfing species Arabidopsis thaliana. PCR amplifications of 12 microsatellite loci among 49 natural populations revealed between one to eight length variants (alleles) for each locus. The average number of alleles per locus was four and the average genetic diversity index was 0.43. Divergence between length variants was investigated at the nucleotide level. Several observations emerge from the sequence data: (1) for most loci, length polymorphism results only from variations in the number of trinucleotide repeats; (2) for a few others, some variability was noted in the flanking sequences; (3) for compound and interrupted loci containing two arrays of trinucleotide repeats, length variations preferentially affect the longest one. Five of the Arabidopsis thaliana accessions were clearly composed of two sublines. In 2 other accessions, some heterozygous individual plants, probably resulting from recent outcrosses, were found. A phylogenetic tree constructed on the basis of trinucleotide microsatellite allelic diversity shows that genetic relationships among the accessions are not correlated with their geographic origin. Received: 4 November 1997 / Accepted: 3 March 1998     

Divergent Human Y-Chromosome Microsatellite Evolution Rates

In this work, we analyze several characteristics influencing the low variability of the microsatellite DYS19 in the major founder Amerindian Y chromosome lineage containing the point mutation DYS199-T. Variation of DYS19 was compared with that of five other Y-linked tetranucleotide repeat loci (DYS389A, DYS389B, DYS390, DYS391, and DYS393) in the DYS199-T lineage. All the other microsatellites showed significantly higher levels of variability than DYS19 as measured by gene diversity and repeat number variance. Moreover, we had previously shown that DYS19 had high diversity in Brazilians and in several other populations worldwide. Thus, the slow DYS19 evolution in the DYS199-T lineage seems to be both locus and allele specific. To understand the slow DYS19 evolutionary rate, the microsatellite loci were compared according to their mapping on the Y chromosome and also on the basis of structural aspects such as the base composition of the repeat motif and flanking regions and the degree of perfection and size (repeat number) of the variable blocks. The only observed difference that might be related to the low DYS19 variability is its small average number of repeats, a value expected to be closer to the founder DYS19 allele in the DYS199-T lineage. These data were also compared to other derived Y lineages. The Tat-C lineage displayed a lower DYS19 variability correlated to a small average repeat number, while in the DYS234-G lineage, a high DYS19 variability was found associated to a larger average repeat number. This approach reveals that evolution of Y microsatellites in lineages defined by slowly evolving markers, such as point mutations, can be greatly influenced by the size (number of repeats of the variable block) of the founder allele in each microsatellite locus. Thus lineage-dating methods using microsatellite variation should be practiced with great care. Received: 7 November 1998 / Accepted: 9 April 1999     

Characterisation and analysis of microsatellite loci in a mangrove species, Avicennia marina (Forsk.) Vierh. (Avicenniaceae)

An enriched microsatellite library of the mangrove species Avicennia marina was constructed, in which 85.8% of the clones contained microsatellite sequences. Of the microsatellite repeat sequences isolated, 55.0% were di-nucleotides, 34.2% were tri-nucleotides, 50.0% were perfect, 24.2% were imperfect, and 15.0% were compound. Four different di-nucleotide repeats were isolated with repeat lengths ranging from 5 to 33; ten different tri-nucleotide repeats were isolated with repeat lengths ranging from 3 to 25. The most common di-nucleotide was the AC/TG repeat; the most common tri-nucleotide was the CCG/GGC repeat. Sixteen microsatellite sequences were selected for primer design, and 6 primers were selected to investigate the polymorphism detected among 15 individuals of A. marina from three natural populations in Australia. A total of 40 alleles were detected at 6 microsatellite loci. The number of alleles per microsatellite locus ranged from 5 to 13. On average, 7 alleles were detected per locus. All microsatellite loci showed high levels of gene diversity (heterozygosity), with values ranging from 0.53 to 0.88; the mean value of gene diversity was 0.70. Microsatellite loci were also tested for conservation across Avicennia species. There was a decline in amplification success with increasing divergence between Avicennia species. The results indicate that microsatellites are abundant in the Avicennia genome and can be valuable genetic markers for assessing the effects of deforestation and forest fragmentation in mangrove communities, which is an important issue for mangrove conservation and afforestation schemes. Received: 8 June 1999 / Accepted: 21 September 1999     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Evolution of Tandemly Repeated Sequences: What Happens at the End of an Array?

Tandemly repeated sequences are a major component of the eukaryotic genome. Although the general characteristics of tandem repeats have been well documented, the processes involved in their origin and maintenance remain unknown. In this study, a region on the paternal sex ratio (PSR) chromosome was analyzed to investigate the mechanisms of tandem repeat evolution. The region contains a junction between a tandem array of PSR2 repeats and a copy of the retrotransposon NATE, with other dispersed repeats (putative mobile elements) on the other side of the element. Little similarity was detected between the sequence of PSR2 and the region of NATE flanking the array, indicating that the PSR2 repeat did not originate from the underlying NATE sequence. However, a short region of sequence similarity (11/15 bp) and an inverted region of sequence identity (8 bp) are present on either side of the junction. These short sequences may have facilitated nonhomologous recombination between NATE and PSR2, resulting in the formation of the junction. Adjacent to the junction, the three most terminal repeats in the PSR2 array exhibited a higher sequence divergence relative to internal repeats, which is consistent with a theoretical prediction of the unequal exchange model for tandem repeat evolution. Other NATE insertion sites were characterized which show proximity to both tandem repeats and complex DNAs containing additional dispersed repeats. An ``accretion model' is proposed to account for this association by the accumulation of mobile elements at the ends of tandem arrays and into ``islands' within arrays. Mobile elements inserting into arrays will tend to migrate into islands and to array ends, due to the turnover in the number of intervening repeats. Received: 18 August 1997 / Accepted: 18 September 1998     

A Phylogenetic Perspective on Sequence Evolution in Microsatellite Loci

We examined the evolution of the repeat regions of three noncoding microsatellite loci in 58 species of the Polistinae, a subfamily of wasps that diverged over 140 million years ago. A phylogenetic approach allows two new kinds of approaches to studying microsatellite evolution: character mapping and comparative analysis. The basic repeat structure of the loci was highly conserved, but was often punctuated with imperfections that appear to be phylogenetically informative. Repeat numbers evolved more rapidly than other changes in the repeat region. Changes in number of repeats among species seem consistent with the stepwise mutation model, which is based on slippage during replication as the main source of mutations. Changes in repeat numbers can occur even when there are very few tandem repeats but longer repeats, especially perfect repeats led to greater rates of evolutionary change. Species phylogenetically closer to the one from which we identified the loci had longer stretches of uninterrupted repeats and more different motifs, but not longer total repeat regions. The number of perfect repeats increased more often than it decreased. However, there was no evidence that some species have consistently greater numbers of repeats across loci than other species have, once ascertainment bias is eliminated. We also found no evidence for a population size effect posited by one form of the directionality hypothesis. Overall, phylogenetic variation in repeat regions can be explained by adding neutral evolution to what is already known about the mutation process. The life cycle of microsatellites appears to reflect a balance between growth by slippage and degradation by an essentially irreversible accumulation of imperfections. Received: 13 April 1999 / Accepted: 8 September 1999     

Extensive Amplification and Transposition of a Novel Repetitive Element, Xstir, Together with Its Terminal Inverted Repeat in the Evolution of Xenopus

A DNA fragment containing short tandem repeat sequences (approximately 86-bp repeat) was isolated from a Xenopus laevis cDNA library. Southern blot and in situ hybridization analyses revealed that the repeat was highly dispersed in the genome and was present at approximately 1 million copies per haploid genome. We named this element Xstir (Xenopus short tandemly and invertedly repeating element) after its arrangement in the genome. The majority of the genomic Xstir sequences were digested to monomer and dimer sizes with several restriction enzymes. Their sequences were found to be highly homogeneous and organized into tandem arrays in the genome. Alignment analyses of several known sequences showed that some of the Xstir-like sequences were also organized into interspersed inverted repeats. The inverted repeats consisted of an inverted pair of two differently modified Xstirs separated by a short insert. In addition, these were framed by another novel inverted repeat (Xstir-TIR). The Xstir-TIR sequence was also found at the ends of tandem Xstir arrays. Furthermore, we found that Xstir-TIR was linked to a motif characterizing the T2 family which belonged to a vertebrate MITE (miniature inverted-repeat transposable element) family, suggesting the importance of Xstir-TIR for their amplification and transposition. The present study of 11 anuran and 2 urodele species revealed that Xstir or Xstir-like sequences were extensively amplified in the three Xenopus species. Genomic Xstir populations of X. borealis and X. laevis were mutually indistinguishable but significantly different from that of X. tropicalis. Received: 5 April 2000 / Accepted: 3 August 2000     

Mhc Allelic Diversity and Modern Human Origins

Thirty complete coding sequences of human major histocompatibility complex (Mhc) class II DRB alleles, spanning 237 codons, were analyzed for phylogenetic information using distance, parsimony, and likelihood approaches. Allelic genealogies derived from different parts of the coding sequence (exon 2, the 5′ and 3′ ends of exon 2, respectively, and exons 3–6) were compared. Contrary to prior assertions, a rigorous analysis of allelic genealogies in this gene family cannot be used to justify the claim that the lineage leading to modern humans contained on average at least 100,000 individuals. Phylogenetic inferences based upon the exon 2 region of the DRB loci are complicated by selection and recombination, so this part of the gene does not provide a complete and accurate view of allelic relationships. Attempts to reconstruct human history from genetic data must use realistic models which consider the complicating factors of nonequilibrium populations, recombination, and different patterns of selection. Received: 19 February 1997 / Accepted: 12 June 1997     

Conservation and evolution of microsatellite loci in primate taxa

Microsatellites are promising genetic markers for the study of demographic structure and phylogenetic history in populations. However, little information exists on the molecular nature of the repeats and their flanking sequences of a same microsatellite in a large range of species. In this study, we report polymorphism and consensus sequences of eight microsatellite loci using human primers in 20 primate species. The results show size polymorphism in almost all species and microsatellites. These loci are therefore useful markers for population genetic studies between populations of the same species. Insertion/deletion events are frequent in the flanking regions, the majority concerning several contiguous bases. This is in contrast with the more usual single base pair events in non-coding regions. The ranges of allele lengths in non-human primates often show no overlap with that of human, usually due to the deletion/insertion events in the flanking sequences, producing smaller allele lengths rather than smaller numbers of repeats. The use of length of PCR product will bias the inter-species interpretation reducing the number of observable alleles and treating as the same allele very divergent molecular sequences. Caution should be used when employing microsatellites in cross-species comparisons in which the species under study are separated by significant amounts of evolutionary time: in such cases allele comparison cannot be based on lengths alone.     

Mutational mechanisms,phylogeny, and evolution of a repetitive region within a clock gene ofDrosophila melanogaster

The D. melanogaster clock gene period (per) is an internally repetitive gene encoding a tandem array of Thr-Gly codons that are highly polymorphic in length in European natural populations. The two major length variants, (Thr-Gly)₂₀ and (Thr-Gly)₁₇, show a highly significant latitudinal cline. In this study we present the complete sequence of the Thr-Gly region of 91 individuals from 6 natural populations of D. melanogaster, 5 from Europe and 1 from North Africa. We further characterized these 91 individuals for polymorphic sites in two other regions, one upstream and one downstream of the Thr-Gly repeat. We used the haplotypic combinations of Thr-Gly allele with flanking markers in an attempt to identify the mechanisms involved in the evolution of the D. melanogaster Thr-Gly region and to infer the phylogenetic relationship existing among the Thr-Gly alleles. We observe evidence for both intra- and interallelic mutational mechanisms, including replication slippage, unequal crossing-over, and gene conversion. Received: 22 August 1995 / Accepted: 17 October 1995     

A Complex Organization of the Gene Encoding Cytochrome Oxidase Subunit 1 in the Mitochondrial Genome of the Dinoflagellate, Crypthecodinium cohnii: Homologous Recombination Generates Two Different cox1 Open Reading Frames

In the course of investigating mitochondrial genome organization in Crypthecodinium cohnii, a non-photosynthetic dinoflagellate, we identified four EcoRI fragments that hybridize to a probe specific for cox1, the gene that encodes subunit 1 of cytochrome oxidase. Cloning and sequence characterization of the four fragments (5.7, 5.1, 4.1, 3.5 kilobase pairs) revealed that cox1 exists in four distinct but related contexts in C. cohnii mtDNA, with a central repeat unit flanked by one of two possible upstream (flanking domain 1 or 2) and downstream (flanking domain 3 or 4) regions. The majority of the cox1 gene is located within the central repeat; however, the C-terminal portion of the open reading frame extends into flanking domains 3 and 4, thereby creating two distinct cox1 coding sequences. The 3′-terminal region of one of the cox1 reading frames can assume an elaborate secondary structure, which potentially could act to stabilize the mature mRNA against nucleolytic degradation. In addition, a high density of small inverted repeats (15–22 base pairs) has been identified at the 5′-end of cox1, further suggesting that hairpin structures could be important for gene regulation. The organization of cox1 in C. cohnii mtDNA appears to reflect homologous recombination events within the central repeat between different cox1 sequence contexts. Such recombining repeats are a characteristic feature of plant (angiosperm) mtDNA, but they have not previously been described in the mitochondrial genomes of protists. Received: 21 December 2000 / Accepted: 30 January 2001     

Allelic variation in the squirrel monkey x-linked color vision gene: biogeographical and behavioral correlates

Most Neotropical primate species possess a polymorphic X-linked and a monomorphic autosomal color vision gene. Consequently, populations are composed of both dichromatics and trichromatics. Most theories on the maintenance of this genetic system revolve around possible advantages for foraging ecology. To examine the issue from a different angle, we compared the numbers and relative frequencies of alleles at the X-linked locus among three species of Saimiri representing a wide range of geographical and behavioral variation in the genus. Exons 3, 4, and 5 of the X-linked opsin gene were sequenced for a large number of X chromosomes for all three species. Several synonymous mutations were detected in exons 4 and 5 for the originally reported alleles but only a single nonsynonymous change was detected. Two alleles were found that appeared to be the result of recombination events. The low occurrence of recombinant alleles and absence of mutations in the amino acids critical for spectral tuning indicates that stabilizing selection acts to maintain the combinations of critical sites specific to each allele. Allele frequencies were approximately the same for all Saimiri species, with a slight but significant difference between S. boliviensis and S. oerstedii. No apparent correlation exists between allele frequencies and behavioral or biogeographical differences between species, casting doubt on the speculation that the spectral sensitivities of the alleles have been maintained because they are specifically well-tuned to Saimiri visual ecology. Rather, the spectral tuning peaks might have been maintained because they are as widely spaced as possible within the limited range of middlewave to longwave spectra useful to all primates. This arrangement creates a balance between maximizing the distance between spectral tuning peaks (allowing the color opponency of the visual system to distinguish between peaks) and maximizing the number of alleles within a limited range (yielding the greatest possible frequency of heterozygotes).     

Ancient Origin of the Null Allele se 428 of the Human ABO-Secretor Locus (FUT2)

In human populations, a null allele having several nucleotide differences from the wild-type allele is segregating at the FUT2 locus (the ABO-Secretor locus) encoding α(1,2)fucosyltransferase. To estimate the age of the most recent common ancestor (MRCA) of these two alleles, we sequenced FUT2 homologues from chimpanzee, gorilla, orangutan, and green monkey. Since we did not detect acceleration or any heterogeneity in the substitution rate at this locus among these species, the age of the MRCA was estimated to be around 3 MYA, assuming the divergence time of human and chimpanzee to be 5 MYA. We developed a simple test to examine whether or not the old age of the MRCA of the FUT2 is consistent with that expected for two divergent neutral alleles sampled from a random mating population. An application of the test to the data at FUT2 indicated that the age of the MRCA is too old to be explained by the simple neutral assumptions, although our test depends on accurate estimation of the divergence time of human and chimpanzee in units of twice the human population size. Various possibilities including balancing selection are discussed to explain this old age of the MRCA. Received: 9 May 1999 / Accepted: 20 September 1999     

Characterization of the Sol3 family of nonautonomous transposable elements in tomato and potato

Sol3 transposons are mobile elements defined by long terminal inverted repeats which are found in tomato and potato. Members of the Sol3 family have been isolated from a variety of solanaceous species including Solanum tuberosum (potato), S. demissum, S. chacoense, Lycopersicon esculentum (tomato), and L. hirsutum. While highly conserved elements are found within different species, Sol3 terminal inverted repeats can also flank unrelated sequences. Southern blot analysis indicates that Sol3 elements are less prevalent in the potato (approximately 50 copies) than in the tomato (>100 copies) genome. No Sol3-hybridizing sequences were observed in tobacco. While a number of Sol3 elements ranging in size from 500 bp to 2 kbp were sequenced, no transposase coding domains could be identified within the internal regions of the elements. The data suggest that the Sol3 represent a heterogeneous family of nonautonomous transposable elements associated with an as-yet-unidentified autonomous transposon. Received: 18 September 1996 / Accepted: 11 March 1997     

Absence of Protein Polymorphism in the Ras Genes of Drosophila melanogaster

Sequence analysis of 27 alleles of each of the three Ras-related genes in Drosophila melanogaster indicates that they all have low levels of polymorphism but may experience slightly different evolutionary pressures. No amino acid replacement substitutions were indicated in any of the sequences, or in the sibling species D. simulans and D. mauritiana. The Dras1 gene, which is the major ras homologue in Drosophila, has less within-species variation in D. melanogaster relative to the amount of divergence from the sibling species than does Dras2, although the contrast was not significant by the HKA test. Dras2 appears to be maintaining two classes of haplotype in D. melanogaster, one of which is closer to the alleles observed in the sibling species, suggesting that this is not likely to be a pseudogene despite the absence of a mutant phenotype. Although differences in level of expression may affect the function of the genes, it is concluded that genetic variation in the Ras signal transduction pathways cannot be attributed to catalytic variation in the Ras proteins. Received: 5 November 1998 / Accepted: 26 March 1999