Historical effective size and the level of genetic diversity in Drosophila melanogaster and Drosophila pseudoobscura

We report the results of a sequential gel electrophoretic study of protein variation in Drosophila melanogaster and its comparison with D. pseudoobscura. The number of alleles and mean heterozygosity were lower in D. melanogaster than in D. pseudoobscura. On the other hand, geographical populations of Drosophila melanogaster have been shown to be much more differentiated than those of D. pseudoobscura. The results suggest that in D. melanogaster low-frequency alleles have been lost during the colonization process and that major alleles have become differentiated among populations. Population bottlenecks, due to various causes, appear to have played a significant role in the shaping of genetic variation in natural populations of many species. It is proposed that a comparison of genetic variation at homologous gene loci between related species can bring out effects of historical bottlenecks and provide an alternative approach for analyzing causes of genetic variation in natural populations.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Gene Flow and the Geographical Distribution of a Molecular Polymorphism in DROSOPHILA PSEUDOOBSCURA

This paper discusses the relation between the geographical distribution of an enzyme polymorphism and population structure in Drosophila pseudoobscura. California populations of this species living in very different montane and lowland habitats separated by several kilometers are similar to each other in the frequency of an esterase allele. Previous estimates suggest that gene flow is too limited to account for this homogeneity of genetic structure, so that it must reflect some balancing force of natural selection. We show, however, that dispersal over unfavorable habitats is much greater than earlier supposed. Isolated populations of D. pseudoobscura separated by 15 km from other populations are subject to large amounts of immigration. This is shown by changes in the seasonal abundance of this species and in the annual pattern of lethal alleles in such populations. The genetic structure of an experimentally perturbed isolated population in an oasis returned to normal within a single year, suggesting that such populations are ephemeral and that the oasis is subject to annual recolonization by distant migrants. Direct assessment of marked flies shows that they can move at least 10 km in 24 hours over a desert. Such extensive gene flow may help explain the distribution of the esterase allele, and is relevant to the high level of molecular polymorphism and its general lack of geographic differentiation throughout the range of D. pseudoobscura.     

Gene Flow in Natural Populations of DROSOPHILA MELANOGASTER with Special Reference to Lethal Allelism Rates and Protein Variation

A simultaneous survey of 14 protein loci, together with frequencies and within- and between-population allelism rates of lethal chromosomes, was carried out in five (four Japanese and one Korean) natural populations and one cage population of Drosophila melanogaster. It was found that lethal allelism rates decrease rapidly as geographic distance between two populations increases, while variation at protein loci shows a remarkable similarity over all populations examined. These findings suggest that there are very high levels of gene flow in these natural populations and that selection at protein loci which can maintain substantial geographic variation, if present, is overshadowed by gene flow. There is no indication that invasion of D. melanogaster to the Far East occurred so recently that the frequencies of lethal chromosomes are still in nonequilibrium.     

Genic Variation in Abundant Soluble Proteins of DROSOPHILA MELANOGASTER and DROSOPHILA PSEUDOOBSCURA

Genic variation was surveyed for 20 proteins of Drosophila melanogaster and 18 proteins of D. pseudoobscura. Analysis was by extraction and one-dimensional polyacrylamide gel electrophoresis under nondenaturing conditions, followed by staining with Coomassie Brilliant Blue to detect soluble proteins present in relatively large amounts ("abundant soluble proteins"). D. melanogaster was polymorphic for 65% of its protein loci and an individual was heterozygous for 10% of its loci. The respective figures for D. pseudoobscura were 61% and 11%. These estimates of genic variation fall between previously published estimates obtained for these species by one-dimensional electrophoresis of soluble enzymes and those obtained by two-dimensional electrophoresis of solubilized abundant proteins. However, variation for both species could be strongly partitioned between loci, on the basis of tissue and stage expression of the proteins. The results are discussed with respect to their bearing on the possibility that abundant proteins constitute a distinct class of proteins less polymorphic than soluble enzymes.     

A Comprehensive Study of Genic Variation in Natural Populations of Drosophila melanogaster. III. Variations in Genetic Structure and Their Causes between Drosophila melanogaster and Its Sibling Species Drosophila simulans

The natural populations of Drosophila melanogaster and Drosophila simulans were compared for their genetic structure. A total of 114 gene-protein loci were studied in four mainland (from Europe and Africa) and an island (Seychelle) populations of D. simulans and the results were compared with those obtained on the same set of homologous loci in fifteen worldwide populations of D. melanogaster. The main results are as follows: (1) D. melanogaster shows a significantly higher proportion of loci polymorphic than D. simulans (52% vs. 39%, P<0.05), (2) both species have similar mean heterozygosity and mean number of alleles per locus, (3) the two species share some highly polymorphic loci but they do not share loci that show high geographic differentiation, and (4) D. simulans shows significantly less geographic differentiation than D. melanogaster. The differences in genetic differentiation between the two species are limited to loci located on the X and second chromosomes only; loci on the third chromosome show similar level of geographic differentiation in both species. These two species have previously been shown to differ in their pattern of variation for chromosomal polymorphisms, quantitative and physiological characters, two-dimensional electrophoretic (2DE) proteins, middle repetitive DNA and mitochondrial DNA. Variation in niche-widths and/or genetic "strategies" of adaptation appear to be the main causes of differences in the genetic structure of these two species.     

Allozyme variation in Greek wild populations of Drosophila melanogaster and D. simulans along a North-South gradient

Allozyme frequency data from five Greek wild sympatric populations of Drosophila melanogaster and D. simulans along a North-South gradient were analyzed for genotype-environment relationships. The regression coefficient of genetic distance on geographic distance indicates that there is a significant relationship between these parameters for D. melanogaster only. Highly significant differences in specific alleles at certain loci were found between the various local populations studied. The changes in Gpdh ^F of D. melanogaster and Est-6 ^F of D. simulans exhibited clinal patterns in allele frequencies. In addition, analysis of D. melanogaster Gpdh ^Fand Adh ^F allele frequencies shows that the Greek data do not have regression coefficients (regressing allele frequency on degrees North of latitude) of the same sign as East-and West-Coast United-States populations. These contradictory data are discussed in relation to what is known about the maintenance of the Adh and Gpdh polymorphisms.     

An evolutionary model for the duplication and divergence of esterase genes in Drosophila

Summary The esterase 5 (Est-5 = gene, EST 5 = protein) enzyme in Drosophila pseudoobscura is encoded by one of three paralogous genes, Est-5A, Est5B, and Est-5C, that are tightly clustered on the right arm of the X chromosome. The homologous Est-6 locus in Drosophila melanogaster has only one paralogous neighbor, Est-P. Comparisons of coding and flanking DNA sequences among the three D. pseudoobscura and two D. melanogaster genes suggest that two paralogous genes were present before the divergence of D. pseudoobscura from D. melanogaster and that, later, a second duplication occurred in D. pseudoobscura. Nucleotide sequences of the coding regions of the three D. pseudoobscura genes showed 78–85% similarity in pairwise comparisons, whereas the relatedness between Est-6 and Est-P was only 67%. The higher degree of conservation in D. pseudoobscura likely results from the comparatively recent divergence of Est-5B and Est-5C and from possible gene conversion events between Est-5A and Est-5B. Analyses of silent and replacement site differences in the two exons of the paralogous and orthologous genes in each species indicate that common selective forces are acting on all five loci. Further evidence for common purifying selective constraints comes from the conservation of hydropathy profiles and proposed catalytic residues. However, different levels of amino acid substitution between the paralogous genes in D. melanogaster relative to those in D. pseudoobscura suggest that interspecific differences in selection also exist.Offprint requests to: R.C. Richmond     

A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement‐specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.     

Contrasting patterns of geographic variation in the cosmopolitan sibling species Drosophila melanogaster and Drosophila simulans

An electrophoretic study was carried out to compare the geographic pattern of genetic variation in Drosophila simulans with that of its sibling species, Drosophila melanogaster. An identical set of 32 gene-protein loci was studied in four geographically distant populations of D. simulans and two populations of D. melanogaster, all originating from Europe and Africa. The comparison yielded the following results: (1) tropical populations of D. simulans were, in terms of the number of unique alleles, average heterozygosity per locus, and percentage of loci polymorphic, more variable than conspecific-temperate populations; (2) some loci in both species showed interpopulation differences in allele frequencies that suggest latitudinal clines; and (3) temperate-tropical genetic differentiation between populations was much less in D. simulans than in D. melanogaster. Similar differences between these two species have previously been shown for chromosomal, quantitative, physiological, and middle-repetitive DNA variation. Estimates of N _m (number of migrants per generation) from the spatial distribution of rare alleles suggest that both species have similar levels of interpopulation gene flow. These observations lead us to propose two competing hypotheses: the low level of geographic differentiation in D. simulans is due to its evolutionarily recent worldwide colonization and, alternatively, D. simulans has a narrower niche than D. melanogaster. Geographic variation data on different genetic elements (e.g., mitochondrial DNA, two-dimensional proteins, etc.) are required before these hypotheses can be adequately tested.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Odorant receptor (Or) genes: polymorphism and divergence in the D. melanogaster and D. pseudoobscura lineages

Background

In insects, like in most invertebrates, olfaction is the principal sensory modality, which provides animals with essential information for survival and reproduction. Odorant receptors are involved in this response, mediating interactions between an individual and its environment, as well as between individuals of the same or different species. The adaptive importance of odorant receptors renders them good candidates for having their variation shaped by natural selection.

Methodology/Principal Findings

We analyzed nucleotide variation in a subset of eight Or genes located on the 3L chromosomal arm of Drosophila melanogaster in a derived population of this species and also in a population of Drosophila pseudoobscura. Some heterogeneity in the silent polymorphism to divergence ratio was detected in the D. melanogaster/D. simulans comparison, with a single gene (Or67b) contributing ∼37% to the test statistic. However, no other signals of a very recent selective event were detected at this gene. In contrast, at the speciation timescale, the MK test uncovered the footprint of positive selection driving the evolution of two of the encoded proteins in both D. melanogaster —OR65c and OR67a —and D. pseudoobscura —OR65b1 and OR67c.

Conclusions

The powerful polymorphism/divergence approach provided evidence for adaptive evolution at a rather high proportion of the Or genes studied after relatively recent speciation events. It did not provide, however, clear evidence for very recent selective events in either D. melanogaster or D. pseudoobscura.     

P1 clones from Drosophila melanogaster as markers to study the chromosomal evolution of Muller's A element in two species of the obscura group of Drosophila

Thirty P1 clones from the X chromosome (Muller's A element) of Drosophila melanogaster were cross-hybridized in situ to Drosophila subobscura and Drosophila pseudoobscura polytene chromosomes. An additional recombinant phage Dsuby was also used as a marker. Twenty-three (77%) of the P1 clones gave positive hybridization on D. pseudoobscura chromosomes bat only 16 (53%) did so with those of D. subobscura. Eight P1 clones gave more than one hybridization signal on D. pseudoobscura and/or D. subobscura chromosomes. All P1 clones and Dsuby hybridized on Muller's A element (X chromosome) of D. subobscura. In contrast, only 18 P1 clones and Dsuby hybridized on Muller's D element (XR chromosomal arm) of D. pseudoobscura; 4 additional P1 clones hybridized on Muller's D element (XR chromosomal arm) of this species and the remaining P1 clone gave on hybridization signal on each arm of the X chromosome. This latter clone may contain one breakpoint of a pericentric inversion that may account for the interchange of genetic material between Muller's A and D elements in D. pseudoobscura. In contrast to the rare interchange of genetic material between chromosomal elements, profound differences in the order and spacing of markers were detected between D. melanogaster, D. pseudoobscura and D. subobscura. In fact, the number of chromosomal segments delimited by identical markers and conserved between pairwise comparisons is small. Therefore, extensive reorganization within Muller's A element has been produced during the divergence of the three species. Rough estimates of the number of cytologically detectable inversions contributing to differentiation of Muller's A element were obtained. The most reliable of these estimates is that obtained from the D. pseudoobscura and D. melanogaster comparison since a greater number of markers have been mapped in both species. Tentatively, one inversion breakpoint about every 200 kb has been produced and fixed during the divergence of D. pseudoobscura and D. melanogaster.     

Dynamics of Correlated Genetic Systems. V. Rates of Decay of Linkage Disequilibria in Experimental Populations of DROSOPHILA MELANOGASTER

The dynamic behavior of four-locus gametic frequency distributions was studied in five replicate cage populations of Drosophila melanogaster for up to 50 generations. The joint frequency distributions were resolved into gene frequencies and various disequilibrium measures. In addition, F statistics for marginal single-locus genotypic frequency distributions were followed through time. The gene frequency, disequilibrium and F statistics were obtained for four chromosome 3 enzyme marker loci [isocitrate dehydrogenase (3–27.1), esterase-6 (3–36.8), phosphoglucomutase (3–43.4) and esterase-C (3–49.0)]. The initial structure of the experimental populations featured random mating proportions, and two complementary gametic types with respect to the marker loci, thus assuring complete pairwise linkage disequilibrium among the markers.——The experimental results indicate: (1) the between-replicate variance in gene frequency varied substantially among loci, with isocitrate dehydrogenase showing the greatest between-replicate variance, and esterase-C the least. (2) The F statistics initially were strongly negative but decayed to the neighborhood of zero for all marker loci except esterase-C. The rate at which the F statistics approached zero varied among the marker loci, indicating substantial differences in the distribution of selective effects along the chromosome. The centromeric region, marked by esterase-C, shows the strongest selective effects. (3) The rate of decay of linkage disequilibrium was much faster than expected for pairs of neutral loci, averaging 1.82 times the neutral rate over all replicates and pairs of loci. This acceleration, which was observed for all six pairwise combinations of loci, was interpreted as resulting from the interaction between selection and recombination. Our experimental results are consistent with many investigations of linkage disequilibrium in natural populations of Drosophila melanogaster that show little or no disequilibrium among enzyme loci. (4) A fortuitous contamination of two cages revealed an apparent regulatory interaction between the migrant and nonmigrant chromosomes at the esterase-C locus. The migrant chromosomes were very rapidly absorbed into the recipient populations, despite this interaction. This result suggests that the dynamics of migration in populations may be phenomenologically richer than anticipated by simple theory.     

Sequence and evolution of the Drosophila pseudoobscura glycerol-3-phosphate dehydrogenase locus

The Gpdh genomic region has been cloned and sequenced in Drosophila pseudoobscura. A total of 6.8 kb of sequence was obtained, encompassing all eight exons of the gene. The exons have been aligned with the sequence from D. melanogaster, and the rates of synonymous and nonsynonymous substitution have been compared to those of other genes sequenced in these two species. Gpdh has the lowest rate of nonsynonymous substitution yet seen in genes sequenced in both D. pseudoobscura and D. melanogaster. No insertion/deletion events were observed, and the overall architecture of the gene (i.e., intron sites, etc.) is conserved. An interesting amino acid reversal was noted between the D. melanogaster Fast allele and the D. pseudoobscura gene.     

Conservatism of sites of tRNA loci among the linkage groups of severalDrosophila species

Summary The sites of seven tRNA genes (Arg-2, Lys-2, Ser-2b, Ser-7, Thr-3, Thr-4, Val-3b) were studied by in situ hybridization.¹²⁵I-labeled tRNA probes fromDrosophila melanogaster were hybridized to spreads of polytene chromosomes prepared from fourDrosophila species representing different evolutionary lineages (D. melanogaster, Drosophila hydei, Drosophila pseudoobscura, andDrosophila virilis). Most tRNA loci occurred on homologous chromosomal elements of all four species. In some cases the number of hybridization sites within an element varied and sites on nonhomologous elements were found. It was observed that both tRNA ₂ ^Arg and tRNA ₂ ^Lys hybridized to the same site on homologous elements in several species. These data suggest a limited amount of exchange among different linkage groups during the evolution ofDrosophila species.     

Estimation of Genetic Variability in Natural Populations of DROSOPHILA SIMULANS by Two-Dimensional and Starch Gel Electrophoresis

Genic variation in natural populations of Drosophila simulans was surveyed using allozymic and two-dimensional electrophoretic techniques. Consistent with some previous reports, allozymic heterozygosity appeared lower than in the sibling species D. melanogaster (0.07 vs. 0.16). No variation was detected by two-dimensional electrophoresis of 19 lines scored for 70 abundant proteins. This is consistent with reported reductions in estimates of genic heterozygosity by two-dimensional electrophoresis in D. melanogaster, Mus musculus, and man. Although the amount of intraspecific variation detected in abundant proteins was lower than that detected for allozymes in D. simulans and D. melanogaster, the genetic distances between the sibling species calculated from the two data sets are not significantly different (0.35 and 0.20). The allozyme and two-dimensional electrophoresis data confirmed the impression from other measures of genetic variation (mitochondrial DNA restriction maps and inversion polymorphisms) that D. simulans is substantially less variable than D. melanogaster.     

Length polymorphism in the threonine-glycine-encoding repeat region of theperiod gene inDrosophila

Summary Single-fly polymerase chain reaction amplification and direct DNA sequencing revealed high levels of length polymorphism in the threonine-glycine encoding repeat region of theperiod (per) gene in natural populations ofDrosophila melanogaster. DNA comparison of two alleles of identical lengths gave a high number of synonymous substitutions suggesting an ancient time of separation. However detailed examination of the sequences of different Thr-Gly length variants indicated that this divergence could be understood in terms of four deletion/insertion events. InDrosophila pseudoobscura a length polymorphism is observed in a five-amino acid degenerate repeat, which corresponds tomelanogaster's Thr-Gly domain. In spite of the differences betweenD. melanogaster andD. pseudoobscura in the amino acid sequence of the repeats, the predicted secondary structures suggest evolutionary and mechanistic constraints on theper protein of these two species.     

Direct Evidence of Genic Differentiation between Sex Ratio and Standard Gene Arrangements of X Chromosome in DROSOPHILA PSEUDOOBSCURA

Evidence of strong association of alleles at the Esterase-5 and adult Acid phosphatase-6 loci with the Sex Ratio and the Standard X chromosome arrangements in D. pseudoobscura is presented. There is no differentiation at these two loci, either in the Sex Ratio or the Standard arrangement from two widely separated geographic populations. At the Esterase-5 locus, alleles characteristic of one arrangement are found in low frequencies in the other arrangement. It is concluded, from these and other observations that the differentiation between the Sex Ratio and the Standard arrangement, at these two loci, is a result of natural selection favoring different coadapted gene blocks in these arrangements. The average heterozygosity at these two loci in different karyotypes is in the following order: ST/SR > ST/ST >> SR/SR which agrees very well with the fitness estimates of these karyotypes.     

Polymorphism and differentiation of multilocus DNA markers in natural populations of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Using multilocus (RAPD) markers, variation and divergence of genomic DNA was examined in two Drosophila melanogaster populations from Russia and three populations from Ukraine. The populations were found to exhibit high polymorphism at RAPD markers. Estimation of genetic distances between the populations showed low differentiation of geographically distant populations of D. melanogaster. Significant gene flow between the D. melanogaster populations was found, which depended on the geographical distance between them.     

Transgenictools for members of the genus Drosophila with sequenced genomes

The sequencing of the genomes of 12 Drosophila species has created an opportunity for much in the way of comparative molecular analyses amongst these species. To aid that endeavor, we have made several transformation vectors based on the piggyBac transposon with 3xP3-EGFP and -ECFP transgenic markers that should be useful for mutagenesis and establishing the GAL4/UAS system in these species. We have tested the ability of mini-white to be used as a marker for insertional mutagenesis, and have observed mini-white derived pigmentation of the testes sheath in a subset of lines from D. pseudoobscura and D. virilis. We have incorporated a source of piggyBac transposase into nine Drosophila species, and have demonstrated the functionality of these transposase lines for mobilization of marked inserts in vivo. Additionally, we tested the ability of a D. melanogaster nanos enhancer element to drive expression of GAL4 in D. melanogaster, D. simulans, D. erecta, D. yakuba, D. pseudoobscura, and D. virilis. The efficacy of the nos-Gal4 transgene was determined by measuring the response of UAS-EGFPtub in all six species. Our results show that D. melanogaster nos-Gal4 drives expression in other species, to varying degrees, in similar spatiotemporal domains in the ovaries, testes, and embryos as seen in D. melanogaster. However, expression levels are variable, demonstrating the possible need to use species-specific promoters in some cases. In summary, we hope to provide a set of guidelines and basic tools, based upon this work, for both insertional mutagenesis and GAL4/UAS system-based experiments in multiple species of Drosophila.