Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster.

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.     

Molecular and Phenotypic Variation of the Zw Locus Region in Drosophila Melanogaster

Restriction map polymorphism in a 13-kb region of the Zw locus in Drosophila melanogaster was investigated for 64 X chromosome lines with seven 6-cutter and ten 4-cutter restriction enzymes. A total of 203 restriction sites were scored, of which 20 were found to be polymorphic. The estimated nucleotide variation for this region for overall data (pi = 0.003 and 0.001, and theta = 0.003 and 0.002, for 4-cutter and 6-cutter studies, respectively) was smaller than that reported for most regions studied in D. melanogaster. It was found that the Slow allozyme has a larger nucleotide variation and haplotype diversity than the Fast allozyme. Results suggest the relatively recent divergence of the Fast allozyme from the Slow allozyme. Glucose 6-phosphate dehydrogenase (G6PD) activity was measured as a phenotype of the Zw locus. A significant difference in G6PD activity between allozymes was detected. The between-line effect was highly significant within the Slow allozyme, but was not significant within the Fast allozyme. Although a direct causative link could not be established, these results suggest an association between the amounts of quantitative and molecular genetic variation at the Zw locus region.     

Molecular Variation of Adh and P6 Genes in an African Population of Drosophila Melanogaster and Its Relation to Chromosomal Inversions

Four-cutter molecular polymorphism of Adh and P6, and chromosome inversion polymorphism of chromosome II were investigated in 95 isogenic lines of an Ivory Coast population of Drosophila melanogaster, a species assumed to have recently spread throughout the world from a West African origin. The P6 gene showed little linkage disequilibrium with the In(2L)t inversion, although it is located within this inversion. This suggests that the inversion and the P6 locus have extensively exchanged genetic information through either double crossover or gene conversion. Allozymic variation in ADH was in linkage disequilibrium with In(2L)t and In(2R)NS inversions. Evidence suggests either that inversion linkage with the Fast allele is selectively maintained, or that this allele only recently appeared. Molecular polymorphism at the Adh locus in the Ivory Coast is not higher than in North American populations. New haplotypes specific to the African population were found, some of them connect the ``Wa(s)-like' haplotypes found at high frequencies in the United States to the other slow haplotypes. Their relation with In(2L)t supports the hypothesis that Wa(s) recently recombined away from an In(2L)t chromosome which may be the cause of its divergence from the other haplotypes.     

The Coalescent Process in Models with Selection and Recombination

The statistical properties of the process describing the genealogical history of a random sample of genes at a selectively neutral locus which is linked to a locus at which natural selection operates are investigated. It is found that the equations describing this process are simple modifications of the equations describing the process assuming that the two loci are completely linked. Thus, the statistical properties of the genealogical process for a random sample at a neutral locus linked to a locus with selection follow from the results obtained for the selected locus. Sequence data from the alcohol dehydrogenase (Adh) region of Drosophila melanogaster are examined and compared to predictions based on the theory. It is found that the spatial distribution of nucleotide differences between Fast and Slow alleles of Adh is very similar to the spatial distribution predicted if balancing selection operates to maintain the allozyme variation at the Adh locus. The spatial distribution of nucleotide differences between different Slow alleles of Adh do not match the predictions of this simple model very well.     

DNA polymorphism at the cytosolic phosphoglucose isomerase (PgiC) locus of the wild plant Arabidopsis thaliana

DNA variation in a 4.7-kb region of the cytosolic phosphoglucose isomerase (PgiC) locus was investigated for 21 ecotypes of Arabidopsis thaliana. The estimated nucleotide diversity was 0.0038, which was one-third of those in previously investigated loci. Since most of the nucleotide variations (93%) were singleton and doubleton, Tajima's test statistic was significantly negative. About 50% of nucleotide polymorphisms in exons were replacement, which caused significance in McDonald and Kreitman's test when compared with Arabis gemmifera and Cardaminopsis petraea. These results indicated that DNA polymorphism at the PgiC locus was not under neutrality. There were two divergent sequence types in the PgiC region, which were associated with allozyme variation. The Fast allozyme was shown to have originated from the Slow allozyme, since two outgroup species had the Slow form. A phylogenetic tree of ecotypes with the Fast allozyme had the shape of a star phylogeny. Mismatch distribution of the Fast allozyme ecotypes resembled that expected under an expanding population model. These results suggest positive selection for the Fast allozyme of the PGIC in A. thaliana.     

Comparative analysis of S haplotypes with very similar SLG alleles in Brassica rapa and Brassica oleracea

Self-incompatibility in Brassica is controlled by a single multi-allelic locus (the S locus) which harbors at least two highly polymorphic genes, SLG and SRK. SRK is a putative transmembrane receptor kinase and its amino acid sequence of the extracellular domain of SRK (the S domain) exhibits high homology to that of SLG. The amino acid sequences of the SLGs of S8 and S46 haplotypes of B. rapa are very similar and those of S23 and S29 haplotypes of B. oleracea were also found to be almost identical. In both cases, SLG and the S domain of SRK of the same haplotype were less similar. This seems to contradict the idea that SLG and SRK of the same haplotype have the same self-recognition specificity. In the transmembrane-kinase domain, the SRK alleles of the S8 and S46 haplotypes had almost identical nucleotide sequences in spite of their lower homology in the S domain. Such a cluster of nucleotide substitutions is probably due to recombination or related events, although recombination in the S locus is thought to be suppressed. Based on our observations, the recognition mechanism and the evolution of self-incompatibility in Brassica are discussed.     

The genetic basis of adaptive pigmentation variation in Drosophila melanogaster

In a broad survey of Drosophila melanogaster population samples, levels of abdominal pigmentation were found to be highly variable and geographically differentiated. A strong positive correlation was found between dark pigmentation and high altitude, suggesting adaptation to specific environments. DNA sequence polymorphism at the candidate gene ebony revealed a clear association with the pigmentation of homozygous third chromosome lines. The darkest lines sequenced had nearly identical haplotypes spanning 14.5 kb upstream of the protein-coding exons of ebony. Thus, natural selection may have elevated the frequency of an allele that confers dark abdominal pigmentation by influencing the regulation of ebony.     

Molecular population genetics of ref(2)P, a locus which confers viral resistance in Drosophila

The ref(2)P locus (2-54.2) is polymorphic for two allelic forms in natural populations of Drosophila melanogaster, ref(2)Po and ref(2)Pp. The latter allele confers resistance to the rhabdovirus sigma infecting wild populations. Previous work, based on a small sample of prescreened restrictive (resistant) and permissive (susceptible) alleles, identified a large number of amino acid replacement changes (7) relative to synonymous changes (1). Such protein variability could be the result of variation-enhancing selection. To further test the selection hypothesis, we have examined the DNA sequences of ten randomly chosen lines of D. melanogaster and one line of D. simulans. Nine of the ten lines are permissive; D. simulans does not harbor the virus. The melanogaster alleles contain 4 synonymous changes, 19 noncoding changes, and 13 amino acid replacement changes, indicating a relatively high level of polymorphism. Three sequenced restrictive alleles have nearly identical sequences, indicating that they are relatively young. Compared to the permissive alleles, they share only a complex deletion at codon 34, CAG-AAT to GGA, which our analysis indicates to be the site conferring the restrictive phenotype. Patterns of polymorphism and divergence differ from neutral predictions by several criteria for the amino terminal region, which contains the complex deletion (codons 1-91), but not the remainder of the protein (codons 92-599). We find a higher rate of evolution on the D. melanogaster lineage than on the D. simulans lineage. The relatively large amount of both replacement and silent polymorphism in the permissive alleles and the lack of divergence between permissive and restrictive alleles suggests that the sigma virus and ref(2)P may be engaged in an evolutionary race in which new restrictive alleles are continually arising but are relatively short-lived.      

Molecular variation at the In(2L)t proximal breakpoint site in natural populations of Drosophila melanogaster and D. simulans

A previous study of nucleotide polymorphism in a Costa Rican population of Drosophila melanogaster found evidence for a nonneutral deficiency in the number of haplotypes near the proximal breakpoint of In(2L)t, a common inversion polymorphism in this species. Another striking feature of the data was a window of unusually high nucleotide diversity spanning the breakpoint site. To distinguish between selective and neutral demographic explanations for the observed patterns in the data, we sample alleles from three additional populations of D. melanogaster and one population of D. simulans. We find that the strength of associations among sites found at the breakpoint varies between populations of D. melanogaster. In D. simulans, analysis of the homologous region reveals unusually elevated levels of nucleotide polymorphism spanning the breakpoint site. As with American populations of D. melanogaster, our D. simulans sample shows a marked reduction in the number of haplotypes but not in nucleotide diversity. Haplotype tests reveal a significant deficiency in the number of haplotypes relative to the neutral expectation in the D. simulans sample and some populations of D. melanogaster. At the breakpoint site, the level of divergence between haplotype classes is comparable to interspecific divergence. The observation of interspecific polymorphisms that differentiate major haplotype classes in both species suggests that haplotype classes at this locus are considerably old. When considered in the context of other studies on patterns of variation within and between populations of D. melanogaster and D. simulans, our data appear more consistent with the operation of selection than with simple demographic explanations.     

Genetic polymorphism at two linked loci,Sod and Est-6, in Drosophila melanogaster

We have examined the patterns of polymorphism at two linked loci, Sod and Est-6, separated by nearly 1000 kb on the left arm of chromosome 3 of Drosophila melanogaster. The evidence suggests that natural selection has been involved in shaping the polymorphisms. At the Sod locus, a fairly strong (s>0.01) selective sweep, started ≥2600 years ago, increased the frequency of a rare haplotype, F(A), to about 50% frequency in populations of Europe, Asia, and the Americas. More recently, an F(A) allele mutated to an S allele, which has increased to frequencies 5–15% in populations of Europe, Asia and North America. All S alleles are identical (or very nearly) in sequence and differ by one nucleotide substitution (which accounts for the F→S electrophoretic difference) from F(A) alleles. At the Est-6 locus, the evidence indicates both directional and balancing selection impacting separately the promoter and the coding regions of the gene, with linkage disequilibrium occurring within each region. Some linkage disequilibrium also exists between the two genes.     

Measuring Selection Coefficients Affecting the Alcohol Dehydrogenase Polymorphism in DROSOPHILA MELANOGASTER

This paper describes a perturbation experiment on the frequency of the F and S Alcohol dehydrogenase (Adh) alleles of D. melanogaster. Fifty-four iso-female lines set up from three wild populations and with initial F frequencies of either 0.25, 0.50 or 0.75 were maintained on standard laboratory food medium at 22 degrees. At generations 4, 12 and 20 the lines were again scored for Adh gene frequencies. Maximum likelihood procedures were used to estimate selection coefficients for the Adh genotypes. An analysis of deviance was used to compare the coefficients against expectations under the hypotheses of neutrality and of constant values for the three base populations, and for the three initial gene frequency classes. Highly-significant departures from neutrality were observed; over all 54 lines, the set of relative fitnesses for S/S:F/S:F/F was estimated as 1.00:1.08:1.08. In addition, there were significant differences between lines in the outcome of selection which were not attributable to differences between base populations or initial F frequencies. These residual between-line differences, as well as some between-generation, within-line differences are discussed in terms of linkage disequilibria with background genes and electrophoretically cryptic variation at the Adh locus.     

Diversity and locus specificity of chicken MHC B class I sequences

The major histocompatibility complex B (MHC B) region in a standard haplotype of Leghorn chickens contains two closely linked class I loci, B-FI and B-FIV. Few sequences of B-FI alleles are available, and therefore alleles of the two loci have not been compared with regard to sequence diversity or locus specificity. Here, we report eight new B-F alpha 1/alpha 2-coding sequences from broiler chicken MHC B haplotypes, and a unique recombinant between the two B-F loci. The new sequences were combined with existing B-F sequences from Leghorn and broiler haplotypes for analysis. On the basis of phylogenetic analysis and conserved sequence motifs, B-F sequences separated into two groups (Groups A and B), corresponding to B-FIV and B-FI locus, respectively. Every broiler haplotype had one B-F sequence in Group A and the second B-F sequence, if it existed, clustered in Group B. Group B (presumptive B-FI locus) sequences identified in broiler haplotypes resembled the human MHC class I HLA-C locus in their distinctive pattern of allelic polymorphism. Compared with B-FIV, B-FI alleles were less polymorphic and possessed a conserved locus-specific motif in the alpha1 helix, but nevertheless demonstrated evidence of diversifying selection. One B-FI alpha 1/alpha 2-coding nucleotide sequence was completely conserved in four different broiler haplotypes, but each allele differed in the exon encoding the alpha 3 domain.     

Intrahaplotype polymorphism at the Brassica S locus

The S locus receptor kinase and the S locus glycoproteins are encoded by genes located at the S locus, which controls the self-incompatibility response in Brassica. In class II self-incompatibility haplotypes, S locus glycoproteins can be encoded by two different genes, SLGA and SLGB. In this study, we analyzed the sequences of these genes in several independently isolated plants, all of which carry the same S haplotype (S(2)). Two groups of S(2) haplotypes could be distinguished depending on whether SRK was associated with SLGA or SLGB. Surprisingly, SRK alleles from the two groups could be distinguished at the sequence level, suggesting that recombination rarely occurs between haplotypes of the two groups. An analysis of the distribution of polymorphisms along the S domain of SRK showed that hypervariable domains I and II tend to be conserved within haplotypes but to be highly variable between haplotypes. This is consistent with these domains playing a role in the determination of haplotype specificity.     

Mitochondrial DNA in the Drosophila Melanogaster Complex

Mitochondrial DNA in the complex Drosophila melanogaster was among the first studied in metazoans. The variability of the molecule was extensively studied using restriction enzymes, gene sequence and recently sequence of the whole coding region. Within the complex, seven major haplotypes have been described, one (me) in D. melanogaster, three in D. simulans (siI, siII, siIII), two in D. mauritiana (maI, maII), and one in D. sechellia (se). The molecular distance between the haplotypes is comprised between 1 and 5%, except for siII and maI, which are virtually identical. The nucleotide diversity within each of these haplotypes is very low, varying from 0 to 0.0005. Most of the cytoplasms are infected by the bacterium Wolbachia and different bacterial strains infect cytoplasms harboring different mtDNA types. mtDNA polymorphism is discussed in relation with Wolbachia, nuclear polymorphism and speciation events.     

Sixteen newH-2 haplotypes derived from wild mice

Wild mice captured in Texas, Scotland, Federal Republic of Germany, Denmark, Spain, Greece, Israel, Egypt, and Chile were mated to inbred strains and through successive backcross matings and H-2 typing lines homozygous for wild-derived H-2, haplotypes were established. The lines, which are neither congenic nor inbred, were then typed with antibodies defining known H-2 alleles at class I and class II loci. In addition, antisera were produced by the immunization of inbred strains with tissues of the new lines. Sixteen of the lines were characterized in this manner. The characterization resulted in the identification of 16 new H-2 haplotypes, 11 new K alleles, 10 new D alleles, and 21 new class I antigenic determinants, most of them of the private type. Most of the haplotypes represent natural recombinants sharing segments of the H-2 complex with previously identified haplotypes. A number of haplotypes are recombinants between the K and the A loci, which in genetic studies have proved difficult to separate. The lines, however, also provide evidence for preservation of blocks of genes in the H-2 complex, particularly in the class II region. Some of class I alleles previously found in wild mice from Michigan have now been found again in these mice. Several class II alleles of these lines appear to be the same as those found in inbred strains. Identical or nearly identical class I and class II alleles thus commonly occur in different populations. These findings strengthen the argument that in populations, H-2 alleles are relatively stable.