Developmental Variation in Amylase Allozyme Activity Associated with Chromosome Inversions in DROSOPHILA PERSIMILIS

The amylase locus in Drosophila persimilis is polymorphic for allozymes, two of which show associations with naturally occurring chromosome 3 inversions. Amy1.09 occurs at high frequencies only in Whitney (WT), while the other common arrangements-Standard (ST), Klamath (KL) and Mendocino (MD)-are predominantly Amy 1.00. We have examined numerous strains, representing various electromorphs and inversions, for variation in cis-specific activity expression in both third-instar larvae and adults. Comparisons of these two life stages also allows the survey of developmental variation in amylase activities. The amount of activity variation exceeds electrophoretic variation at this locus. Moreover, this variation is largely nonrandom and reveals more genic divergence among inversions. The 1.00 allozyme of MD is more active than 1.00 KL in larvae and adults and shows a different developmental pattern. The activity of the 1.00 allozyme of KL is greater than 1.00 allozyme of ST in larvae and adults, but these two arrangements have similar developmental patterns. WT 1 with a 1.00 allele is dramatically different from the 1.00 allozymes of other arrangements in its developmental pattern. The 1.09 allozymes has high activity in WT and KL, but these arrangements differ in their developmental pattern of expression, WT being more active in adults. F2 segregational analyses are consistent with the variation being due to either structural enzyme variants or closely linked cis-acting regulatory elements. We argue that the suppression of recombination between arrangements has allowed the divergence in amylase activity among inversions.     

Distinctions among Allelic Variants Associated with Chromosome 3 Inversions in DROSOPHILA PSEUDOOBSCURA and DROSOPHILA PERSIMILIS

Efforts were made to discriminate new genetic variants among electrophoretic alleles that are associated with chromosome 3 inversions of Drosophila pseudoobscura and D. persimilis. Apparent genetic similarities for electrophoretic alleles between these two species and among the common inversions they carry were reexamined by altering gel concentration and buffer pH. At the amylase locus, the 1.09 electrophoretic allele could be further separated into two allelic classes that differentiated the WT and KL arrangements. Similarly, the 0.84 electrophoretic allele was divided into two allelic classes, one characteristic of the Santa Cruz phylad arrangements, TL and SC, and the other found in strains of the Standard phylad arrangements and CH. Uncommon amylase alleles proved to be different alleles in the two species. No new allelic variants, however, could be found among strains with the amylase 1.00 allele, the commonest allele in the Standard phylad of both species. No major new allelic variation was detected for acid phosphatase-3 and larval protein-10 that revealed any further differentiation among species or inversions. Variation at all three loci in strains of the Bogota population remained genetically similar to variation in strains of mainland D. pseudoobscura.     

Examination of Allelic Variation at the Hexokinase Loci of DROSOPHILA PSEUDOOBSCURA and D. PERSIMILIS by Different Methods

Recently a number of electrophoretic techniques have been applied to reveal the presence of additional genetic variation among the electrophoretic mobility classes of the highly polymorphic xanthine dehydrogenase (XDH ) and esterase-5 (est-5) loci. We examined the hexokinase loci of Drosophila pseudoobscura and D. persimilis using a variety of techniques to determine whether further allelic variation could be revealed for these much less polymorphic loci and to analyze the nature of the known variation at the hexokinase-1 (hex-1) locus. The following studies were conducted: 135 strains of the two species from six localities were examined with buffer pH ranging from 5.5 to 10.0; 40 strains of D. pseudoobscura and 9 strains of D. persimilis from Mather were studied using starch gel concentrations ranging from 8.5 to 15.5% and were examined for differences in heat stability and reactivity to the thiol reagent pCMSA; strains were also tested for susceptibility to urea denaturation and differences in relative activities. Major findings of the work are: (1) No additional allelic variation could be detected at any of the hexokinase loci by applying these techniques. The finding of abundant hidden genetic variation in XDH and est-5 does not extend to all enzyme loci. (2) Evidence from studies using pCMSA indicates that the hex-1 alleles 0.6, 0.8, 1.0 and 1.2 of the two species form a series of unit charge steps. Since the 0.94 allele of D. persimilis has mobility intermediate between 0.8 and 1.0, it is argued that routine electrophoretic techniques are sensitive to at least some conservative amino acid substitutions. (3) Strong correlations were found at the hex-1 locus between low allelic frequency, reduced relative activity and reduced stability to heat and urea denaturation. Since the three sibling species, D. pseudoobscura, D. persimilis and D. miranda, all appear to share a common high frequency allele (1.0) at that locus, these findings are taken as evidence that the observed allelic frequencies are a result of directional selection and mutation, rather than any form of balancing selection.     

Genomic impacts of chromosomal inversions in parapatric Drosophila species

Chromosomal inversions impact genetic variation and facilitate speciation in part by reducing recombination in heterokaryotypes. We generated multiple whole-genome shotgun sequences of the parapatric species pair Drosophila pseudoobscura and Drosophila persimilis and their sympatric outgroup (Drosophila miranda) and compared the average pairwise differences for neutral sites within, just outside and far outside of the three large inversions. Divergence between D. pseudoobscura and D. persimilis is high inside the inversions and in the suppressed recombination regions extending 2.5 Mb outside of inversions, but significantly lower in collinear regions further from the inversions. We observe little evidence of decreased divergence predicted to exist in the centre of inversions, suggesting that gene flow through double crossovers or gene conversion is limited within the inversion, or selection is acting within the inversion to maintain divergence in the face of gene flow. In combination with past studies, we provide evidence that inversions in this system maintain areas of high divergence in the face of hybridization, and have done so for a substantial period of time. The left arm of the X chromosome and chromosome 2 inversions appear to have arisen in the lineage leading to D. persimilis approximately 2 Ma, near the time of the split of D. persimilis-D. pseudoobscura-D. miranda, but likely fixed within D. persimilis much more recently, as diversity within D. persimilis is substantially reduced inside and near these two inversions. We also hypothesize that the inversions in D. persimilis may provide an empirical example of the 'mixed geographical mode' theory of inversion origin and fixation, whereby allopatry and secondary contact both play a role.     

The Speciation History of Drosophila Pseudoobscura and Close Relatives: Inferences from DNA Sequence Variation at the Period Locus

Thirty-five period locus sequences from Drosophila pseudoobscura and its siblings species, D. p. bogotana, D. persimilis, and D. miranda, were studied. A large amount of variation was found within D. pseudoobscura and D. persimilis, consistent with histories of large effective population sizes. D. p. bogotana, however, has a severe reduction in diversity. Combined analysis of per with two other loci, in both D. p. bogotana and D. pseudoobscura, strongly suggest this reduction is due to recent directional selection at or near per within D. p. bogotana. Since D. p. bogotana is highly variable and shares variation with D. pseudoobscura at other loci, the low level of variation at per within D. p. bogotana can not be explained by a small effective population size or by speciation via founder event. Both D. pseudoobscura and D. persimilis have considerable intraspecific gene flow. A large portion of one D. persimilis sequence appears to have arisen via introgression from D. pseudoobscura. The time of this event appears to be well after the initial separation of these two species. The estimated times since speciation are one mya for D. pseudoobscura and D. persimilis and 2 mya since the formation of D. miranda.     

Genetics of hybrid male sterility among strains and species in the Drosophila pseudoobscura species group

Taxa in the early stages of speciation may bear intraspecific allelic variation at loci conferring barrier traits in hybrids such as hybrid sterility. Additionally, hybridization may spread alleles that confer barrier traits to other taxa. Historically, few studies examine within- and between-species variation at loci conferring reproductive isolation. Here, we test for allelic variation within Drosophila persimilis and within the Bogota subspecies of D. pseudoobscura at regions previously shown to contribute to hybrid male sterility. We also test whether D. persimilis and the USA subspecies of D. pseudoobscura share an allele conferring hybrid sterility in a D. pseudoobscura bogotana genetic background. All loci conferred similar hybrid sterility effects across all strains studied, although we detected some statistically significant quantitative effect variation among D. persimilis alleles of some hybrid incompatibility QTLs. We also detected allelism between D. persimilis and D. pseudoobscura USA at a second chromosome hybrid sterility QTL. We hypothesize that either the QTL is ancestral in D. persimilis and D. pseudoobscura USA and lost in D. pseudoobscura bogotana, or gene flow transferred the QTL from D. persimilis to D. pseudoobscura USA. We discuss our findings in the context of population features that may contribute to variation in hybrid incompatibilities.     

An Experimental Investigation of the Unit Charge Model of Protein Polymorphism and Its Relation to the Esterase-5 Locus of DROSOPHILA PSEUDOOBSCURA, DROSOPHILA PERSIMILIS, and DROSOPHILA MIRANDA

The relationship between charge changes and electrophoretic mobility changes is investigated experimentally. The charge of several proteins is altered by reaction with small molecules of known structure and the change in electrophoretic mobility is measured. The method of Ferguson plots is used to separate charge and shape components of mobility differences. The average effect of an amino acid charge change on the mobility of the esterase-5( 1.00) allele of Drosophila pseudoobscura is estimated to be 0.046. This estimate is then used to apply the step model of Ohta and Kimura (1973) to electrophoretic mobility data for the esterase-5 locus of D. pseudoobscura and D. miranda. The variation in electrophoretic mobility at this locus was found to be in agreement with the predictions of the step model.     

Molecular evolution of the Sex-Ratio inversion complex in Drosophila pseudoobscura: analysis of the Esterase-5 gene region

The Sex-Ratio chromosome in Drosophila pseudoobscura is subject to meiotic drive. It is associated with a series of three nonoverlapping paracentric inversions on the right arm of the X chromosome. The esterase-5 gene region has been localized to section 23 within the subbasal inversion of the Sex-Ratio inversion complex, making esterase- 5 a convenient locus for molecular evolutionary analyses of the Sex- Ratio inversion complex and the associated drive system. A 504-bp fragment of noncoding, intergenic DNA from the esterase-5 gene region was amplified and sequenced from 14 Sex-Ratio and 14 Standard X chromosomes of D. pseudoobscura, and from 9 X chromosomes of its two sibling species, Drosophila persimilis and Drosophila miranda. There is extensive sequence differentiation between the Sex-Ratio and Standard chromosomal types. The common Standard chromosome is highly polymorphic, while, as expected from either the neutral mutation theory or the selective sweep hypothesis, the rarer Sex-Ratio chromosome has much less within-chromosome nucleotide polymorphism. We estimate that the Standard and Sex-Ratio chromosomes in D. pseudoobscura diverged between 700,000 and 1.3 Mya, or at least 2 million generations ago. The clustering of D. pseudoobscura Sex-Ratio chromosomes in a neighbor- joining phylogeny indicates a fairly old, monophyletic origin in this species. It appears from these data that Sex-Ratio genes were present prior to the divergence of D. pseudoobscura and D. persimilis and that both the Standard and Sex-Ratio chromosomes of D. persimilis were derived from the Standard chromosome of D. pseudoobscura after the inversion events that isolated the D. pseudoobscura Sex-Ratio chromosome.      

Unusual pattern of single nucleotide polymorphism at the exuperantia2 locus of Drosophila pseudoobscura

We have investigated the pattern of DNA sequence variation at the exuperantia2 locus in Drosophila pseudoobscura. This adds to the increasing dataset of genetic variation in D. pseudoobscura, a useful model species for evolutionary genetic studies. The level of silent site nucleotide diversity and the divergence from an outgroup Drosophila miranda are comparable with those for other X-linked loci. One peculiar pattern at the exu2 locus of D. pseudoobscura is a complete linkage disequilibrium between two SNPs, one of which is a replacement site. As a result, there are two distinct haplotype groups in our dataset. Based upon the comparisons with the outgroup sequences from D. miranda and Drosophila persimilis, we show that the newly derived haplotype group has lower diversity than the ancestral haplotype group. The pattern of protein evolution at exu2 shows some deviation from the neutral model. Together, these and other characteristics of the exu2 locus suggest the action of selection on the pattern of SNP variation, consistent with a partial selective sweep associated with the newly derived haplotype.     

Amylase gene expression in intraspecific and interspecific somatic transformants of Drosophila.

The Amylase locus in Drosophila melanogaster normally contains two copies of the structural gene for alpha-amylase, a centromere-proximal copy, Amy-p, and a distal copy, Amy-d. Products of the two genes may display discrete electrophoretic mobilities, but many strains known to carry the Amy duplication are characterized by a single amylase electromorph, e.g., Oregon-R, which produces the mobility variant AMY-1. A transient expression assay was used in somatic transformation experiments to test the functional status of the Amy genes from an Oregon-R strain. Plasmid constructs containing either the proximal or distal copy were tested in amylase-null hosts. Both genes produced a functional AMY-1 isozyme. Constructs were tested against an AMY-3 reference activity produced by a coinjected plasmid that contains the Amy-d3 allele from a Canton-S strain. With reference to the internal control, the Amy-p and Amy-d genes from Oregon-R expressed different relative activity levels for AMY-1 in transient assays. The transient expression assay was successfully used to test the functional status of Amy-homologous sequences from strains of other species of Drosophila characterized by a single amylase elctromorph, namely, Drosophila pseudoobscura ST and Drosophila miranda S 204. The amylase-null strain of D. melanogaster provided the hosts for these interspecific somatic transformation experiments.     

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.     

Frequency Distribution of Esterase-5 Alleles in Two Populations of DROSOPHILA PSEUDOOBSCURA

Statistical tests comparing allele frequencies in natural populations with those predicted by various theories of genic variation depend critically on the accurate enumeration of alleles. This study used a series of five sequential electrophoretic conditions to characterize the allele frequency distributions of esterase-5 in two large population samples of Drosophila pseudoobscura from California. In Standard chromosome lines 12 electromorphs were discriminated using a single electrophoretic condition. When four additional criteria were used, the number of electromorphs increased to 41, 33 in one population and 22 in the other. Both populations had the same two alleles in high frequency, with other alleles present in frequencies of 6% or less. Although each population had a number of unique alleles, a χ² contingency test demonstrated no significant genetic divergence between them. A statistical comparison of allele frequencies in both populations with that predicted by neutral models suggests that the individual and combined distributions deviate from neutrality in the direction of purifying selection.—Sex-Ratio chromosomes differed markedly from Standard chromosomes in both allelic content and diversity. In 32 Sex-Ratio chromosomes from one population only three alleles were found, all of which were detected under the initial "standard" electrophoretic conditions. Moreover, none of these alleles was found in the Standard chromosome lines.     

Dosage compensation and dietary glucose repression of larval amylase activity inDrosophila miranda

The functional locus for alpha-amylase (Amy) in Drosophila miranda is in the evolutionarily new X2 chromosome. X2 evolved from an autosome in response to an ancestral autosome-Y translocation that gave rise to the "neo-Y" chromosome of this species. Y-linked Amy, if still present in the ancestrally translocated element, is unexpressed. Dosage compensation for amylase activity was examined in larvae of the S 204 strain. Since dietary glucose is known to repress Amy expression in Drosophila melanogaster, dosage compensation of amylase activity in male larvae of D. miranda was tested by rearing larvae of both sexes on yeast diets with or without a glucose supplement. The WT 10 strain of Drosophila persimilis, a sibling species in which Amy is autosomally linked, was used as a reference for tests of amylase activity differences between the sexes. On the diet with glucose, Amy expression was repressed in both WT 10 and S 204 larvae and male larvae of S 204 displayed dosage compensation for amylase activity. On the nonrepressing diet consisting of yeast alone, S 204 continued to display dosage compensation.     

Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajima's D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.     

A Comparative Study of the Esterase-5 Locus in DROSOPHILA PSEUDOOBSCURA,D. PERSIMILIS and D. MIRANDA

Electrophoretic phenotypes of the esterase-5 locus were examined in the sibling species D. pseudoobscura, D. persimilis and D. miranda. D. persimilis alleles were found to have uniformly higher charge on monomers than corresponding alleles of either D. pseudoobscura or D. miranda. Consequently, D. persimilis shares no alleles in common with either D. pseudoobscura or D. miranda, while the latter two species share a number of alleles. It was discovered that by increasing the concentration of acrylamide gel and increasing the length of migration, more allelic differences could be distinguished. Also more alleles were discovered by examining monomer mobility in addition to dimer mobility. In D. persimilis and D. miranda it was found that the previously known high frequency allelic classes broke down into several allelic classes. A test of goodness-of-fit to the infinite alleles model was done and a rough agreement with the model was found.     

Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.     

Mitochondrial DNA variation in Drosophila pseudoobscura and related species in Pacific northwest populations

We have analysed mitochondrial DNA (mtDNA) from Pacific Northwest populations of Drosophila pseudoobscura, D. persimilis, and D. miranda using six restriction enzymes. We find that HpaII restriction sites are hypervariable compared to the other enzymes used. This hypervariability allows construction of a maximum parsimony map linking each mtDNA genotype. Small insertions, possibly tandem duplications, appear to have arisen concomitantly with, or subsequent to, speciation events, perhaps within the A + T rich region. Convergence of mtDNA genotypes is also evident. Unlike findings for other populations of these species, we find little evidence of mitochondrial introgression between D. pseudoobscura and D. persimilis, despite their ability to produce fertile hybrid females.     

Lack of Genic Similarity between Two Sibling Species of Drosophila as Revealed by Varied Techniques

Acrylamide gel electrophoresis was performed on the enzyme xanthine dehydrogenase in sixty isochromosomal lines of Drosophila persimilis from three geographic populations. Sequential electrophoretic analysis using varied gel concentrations and buffers revealed twenty-three alleles in this species where only five had been described previously. These new electrophoretic techniques also detected a profound increase in divergence of gene frequencies at this locus between D. persimilis and its sibling species D. pseudoobscura. The implications of these results for questions of speciation and the maintenance of genetic variability are discussed.     

Variation in Biochemical Properties of Allozymes of Xanthine Dehydrogenase in DROSOPHILA PSEUDOOBSCURA

Twenty-six D. pseudoobscura strains isogenic for xanthine dehydrogenase alleles from Mesa Verde, Colorado, were tested for differences in the biochemical properties of different allelic forms of xanthine dehydrogenase. No significant differences in binding affinity (K_m) or substrate specificity of the enzyme were found. Significant variation among strains, in activity (V_max) and among electromorphs, as well as among strains, in thermolability was found. For the few strains tested, the activity and thermolability differences were shown to co-segregate with the electrophoretic mobility of the variant allele.