Lack of nucleotide polymorphism in the Y-linked sperm flagellar dynein gene Dhc-Yh3 of Drosophila melanogaster and D. simulans

We studied levels of intra- and interspecific nucleotide variation associated with a Y-linked gene in five members of the Drosophila melanogaster subgroup. Using published sequence for 348 bp of the Dhc-Yh3 gene, and degenerate PCR primers designed from comparisons of the sea urchin and Chlamydomonas flagellar dynein genes, we recovered a 1738-bp region in D. melanogaster. Analyses of sequence variation in a worldwide collection of 11 lines of D. melanogaster and 10 lines of D. simulans found only a single silent polymorphism in the latter species. The synonymous site divergence per site for Dhc-Yh3 is comparable to values for X and autosomal genes. Assuming a Wright-Fisher population model, the lack of variation is statistically less than expected using appropriately reduced estimates of theta from the X and autosomes. Because the Y chromosome encodes only six known genes, genetic hitchhiking associated with background selection is unlikely to explain this low variation. Conversely, adaptive hitchhiking, as associated with sex-ratio chromosomes, or a large variance in male fertility may reduce the polymorphism on the Y chromosome. Codon bias is very low, as seen for other genes in regions of low recombination.     

Inviability of hybrids between D. melanogaster and D. simulans results from the absence of simulans X not the presence of simulans Y chromosome

Interspecific crosses between D. melanogaster and D. simulans or its sibling species result in unisexual inviability of the hybrids. Mostly, crosses of D. melanogaster females X D. simulans males produce hybrid females. On the other hand, only hybrid males are viable in the reciprocal crosses. A classical question is the cause of the unisexual hybrid inviability on the chromosomal level. Is it due to the absence of a D. simulans X chromosome or is it due to the presence of a D. simulans Y chromosome? A lack of adequate chromosomal rearrangements available in D. simulans has made it difficult to answer this question. However, it has been assumed that the lethality results from the absence of the D. simulans X rather than the presence of the D. simulans Y. Recently I synthesized the first D. simulans compound-XY chromosome that consists of almost the entire X and Y chromosomes. Males carrying the compound-XY and no free Y chromosome are fertile. By utilizing the compound-XY chromosome, the viability of hybrids with various constitutions of cytoplasm and sex chromosomes has been examined. The results consistently demonstrate that the absence of a D. simulans X chromosome in hybrid genome, and not the presence of the Y chromosome, is a determinant of the hybrid inviability.     

Inviability of hybrids between D. melanogaster and D. simulans results from the absence of simulans X not the presence of simulans Y chromosome.

Interspecific crosses between D. melanogaster and D. simulans or its sibling species result in unisexual inviability of the hybrids. Mostly, crosses of D. melanogaster females x D. simulans males produce hybrid females. On the other hand, only hybrid males are viable in the reciprocal crosses. A classical question is the cause of the unisexual hybrid inviability on the chromosomal level. Is it due to the absence of a D. simulans X chromosome or is it due to the presence of a D. simulans Y chromosome? A lack of adequate chromosomal rearrangements available in D. simulans has made it difficult to answer this question. However, it has been assumed that the lethality results from the absence of the D. simulans X rather than the presence of the D. simulans Y. Recently I synthesized the first D. simulans compound-XY chromosome that consists of almost the entire X and Y chromosomes. Males carrying the compound-XY and no free Y chromosome are fertile. By utilizing the compound-XY chromosome, the viability of hybrids with various constitutions of cytoplasm and sex chromosomes has been examined. The results consistently demonstrate that the absence of a D. simulans X chromosome in hybrid genome, and not the presence of the Y chromosome, is a determinant of the hybrid inviability.     

Polymorphism and divergence at the prune locus in Drosophila melanogaster and D. simulans

The prune locus of Drosophila melanogaster lies at the tip of the X chromosome, in a region of reduced recombination in which nearby loci show reduced variation relative to evolutionary divergence from D. simulans. DNA sequencing of prune alleles from D. melanogaster and D. simulans reveals extremely low variation in D. melanogaster but greater variation in D. simulans. Divergence between the two species is not reduced. This pattern may be explained by either positive selection leading to hitchhiking of neutral variation or background selection against deleterious mutations. The pattern of silent versus replacement polymorphism and divergence at prune is consistent with either a model of weakly deleterious selection against amino acid substitutions or balancing selection.      

Rates of DNA sequence evolution are not sex-biased in Drosophila melanogaster and D. simulans

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

Elevated polymorphism and divergence in the class C scavenger receptors of Drosophila melanogaster and D. simulans

Scavenger receptor proteins are involved in the cellular internalization of a broad variety of foreign material, including pathogenic bacteria during phagocytosis. I find here that nonsynonymous divergence in three class C scavenger receptors (Sr-C's) between Drosophila melanogaster and D. simulans and between each of these species and D. yakuba is approximately four times the typical genome average. These genes also exhibit unusually high levels of segregating nonsynonymous polymorphism in D. melanogaster and D. simulans populations. A fourth Sr-C is comparatively conserved. McDonald-Kreitman tests reveal a significant excess of replacement fixations between D. melanogaster and D. simulans in the Sr-C's, but tests of polymorphic site frequency spectra do not support models of directional selection. It is possible that the molecular functions of SR-C proteins are sufficiently robust to allow exceptionally high amino acid substitution rates without compromising organismal fitness. Alternatively, SR-Cs may evolve under diversifying selection, perhaps as a result of pressure from pathogens. Interestingly, Sr-CIII and Sr-CIV are polymorphic for premature stop codons. Sr-CIV is also polymorphic for an in-frame 101-codon deletion and for the absence of one intron.     

hobo transposable elements in Drosophila melanogaster and D. simulans.

Genomic patterns of occurrence of the transposable element hobo are polymorphic in the sibling species Drosophila melanogaster and D. simulans. Most tested strains of both species have apparently complete (3.0 kb) and smaller hobo elements (H lines), but in both species some strains completely lack such canonical hobo elements (E lines). The occurrence of H and E lines in D. simulans as well as in D. melanogaster implies that an hypothesis of recent introduction in the latter species is inadequate to explain the phylogenetic occurrence of hobo. Particular internally deleted elements, the approximately 1.5 kb Th1 and Th2 elements, are abundant in many lines of D. melanogaster, and an analogous 1.1 kb internally deleted element, h del sim, is abundant in most lines of D. simulans. Besides the canonical hobo sequences, both species (and their sibling species D. sechellia and D. mauritiana) have many hobo-hybridizing sequences per genome that do not appear to be closely related to the canonical hobo sequence.     

Nucleotide variation at the runt locus in Drosophila melanogaster and Drosophila simulans.

Intra- and interspecific nucleotide variation for the major developmental gene runt in Drosophila was studied in D. melanogaster and D. simulans. The 1.5-kb protein-coding region and the 0.4-kb intron of the runt gene were sequenced for 11 alleles in each species. The D. melanogaster alleles originated from east Africa. Estimated parameters of intraspecific variation in D. melanogaster (exons: theta = 0.018, pi = 0.018; intron: theta = 0.014, pi = 0.014) and D. simulans (exons: theta = 0.007, pi = 0.005; intron: theta = 0.008, pi = 0.005) were below average for other X-linked genes, while divergence between species (exons: D = 0.094; intron: D = 0.069) fell within the normal range for both silent and replacement changes. This estimate for runt, along with published values for three other genes in regions of normal recombination, show east African D. melanogaster to be roughly twice as polymorphic as D. simulans. The majority of nucleotide variation, silent and replacement, in both species was found to be selectively neutral using various statistical tests (HKA, McDonald-Kreitman, Tajima, and Fu and Li tests). Monte Carlo simulations of the coalescent process significantly rejected a Wright-Fisher model with respect to an amino acid polymorphism and the distribution of polymorphic sites among the D. simulans lines. This indicated an old lineage and may reflect ancestral population substructuring in D. simulans.     

Multiple allelism of the net gene in Drosophila melanogaster and Drosophila simulans

The net gene mutations are known to cause abnormal pattern of veining in all wing regions except for the first posterior cells. In natural populations of Drosophila melanogaster, the net alleles were identified, which differ in phenotypic expression from standard mutations. The mutants net-extra-analis from a population Belokurikha-2000 have only a single additional vein in the third posterior cell. A line from Chernobyl-1986 population have another nontypical allele netCh86 and shows a lower degree of abnormalities than that usually observed. About 10% of these flies have an additional vein fragment in the first posterior cell. In both males and females of D. simulans population Tashkent-2001, which exhibit netST91 mutation, a net of additional veins is formed as a specific additional fragment in the first posterior cell. The pattern of veining conferred by alleles net-extra-analis and netCh86 is altered to a lesser extent; these alleles are dominant with respect to alleles net2-45 and netST91, which cause more abnormalities. The heterozygotes for alleles netST9 and netCh86 and for Df(2) net62 deletion have an additional fragment in the first posterior cell and show similarly strong deviations from normal wing vein pattern. The natural net alleles correspond, presumably, to different molecular gene defects involved into uncertain local interactions with numerous modifying factors and other genes that specify the wing vein pattern.     

DNA sequence polymorphism and divergence at the erect wing and suppressor of sable loci of Drosophila melanogaster and D. simulans

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (approximately 50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high FST values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.     

Comparison of alcohol dehydrogenase expression in Drosophila melanogaster and D. simulans

Alcohol dehydrogenase (ADH) gene expression was analyzed in Drosophila melanogaster and its sibling species D. simulans. The levels of ADH activity, ADH-cross-reacting material (CRM), and ADH-mRNA were analyzed for several strains of each species, which derive from diverse geographic locations around the world. There is considerable quantitative variation in ADH activity, CRM level, and RNA level among strains within species at all developmental stages. However, the only consistent differences between the two species are in pupal RNA level and in late-adult activity and CRM level. Late-adult melanogaster flies that are homozygous for the Slow allozyme have approximately twice the level of ADH activity and CRM as do simulans flies. The regression of activity on CRM over strains is highly significant and essentially the same for each species, which means that most, if not all, of the activity difference between the species is due to a difference in concentration of the ADH protein. In contrast, there is no significant regression of CRM level on mRNA level in adults of either species; nor is there a significant difference in RNA level between species. Therefore, the difference in ADH protein concentration is not due to RNA template availability. Thus, the interspecific difference in ADH level in adults must be due either to a difference in the rate of translation of the two RNAs or to a difference in protein stability.     

Comparative Life Histories and Ecophysiology of Drosophila melanogaster and D. simulans

Numerous laboratory investigations have compared Drosophila melanogaster and D. simulans for various life history traits and fitness related ecophysiological parameters. From presently available information, it is however difficult to get a general comparative pattern describing the divergence of their ecological niches and understanding their demographic success. Two environmental factors seem however to have played a major role: temperature and alcoholic resources. From an ecophysiological approach, D. simulans may be described as generally more sensitive to stresses; other results point to this species as more cold adapted than its sibling; in some cases, however, D. simulans may appear as better adapted to a warm environment. When investigated, ecophysiological traits show a lesser geographic variability in D. simulans than in D. melanogaster. Presently available information does not explain the ecological prevalence of D. simulans in many places with a mild temperate or subtropical climate. This is presumably due to the fact that most comparisons have been done at a single, standard temperature of 25 degrees C. Comparative studies should be undertaken, spanning the thermal ranges of the two species, and the phenotypic plasticity of ecophysiological traits should now be considered.     

Variation in response to CO2 in Drosophila melanogaster and D. simulans

A type of slow-recovery from high temperature CO2 exposure is a wild-type characteristic of Drosophila simulans. Species hybrids of wild-type D. melanogaster and D. simulans have intermediate recovery times. The response to nitrogen exposure at both high and low temperatures is the same for the two species suggesting that the observed CO2 response is more than just anoxia. Since there is no female bias in hybrids, and since injections of D. simulans extract into D. melanogaster produced no increase in CO2 sensitivity in the recipients, we conclude that the prolonged recovery time results from genomic differences between D. simulans and D. melanogaster, and is not due to a sigma-like infectious agent.     

Low genic variation in male-reproductive-tract proteins of Drosophila melanogaster and D. simulans

We report results, using two-dimensional gel electrophoresis (2DE), of natural population surveys of allelic variation in approximately 300 male-reproductive-tract polypeptides in both Drosophila melanogaster and its sibling species, D. simulans. Despite our efforts to maximize operational sensitivity of our 2DE gels to polymorphism, variation estimates in both species were low (proportion of polymorphic loci [P] = 9%, and average heterozygosity [H] = 1%-3%), compared with those by one-dimensional gel electrophoresis (1DE) (P = 29%-55%; H = 8%-19%) in the same populations. However, H of polymorphic loci was very similar for 2DE and 1DE proteins; and for 17 of a total of 54 polymorphic proteins, 2DE detected three or four distinct alleles. The results suggest that the differing levels of variability widely seen with 1DE and 2DE are real and reflect differing intensities of functional constraint between different classes of structural loci. However, the alternative possibility remains that 2DE has a greater between-locus unevenness of variant detection sensitivity than does 1DE.     

Genetics of Drosophila simulans male mating discrimination in crosses with D. melanogaster

The genetic bases of sexual isolation between Drosophila melanogaster and D. simulans have been mainly studied in females, and there is little information about the role of the males in interspecific mating discrimination. Using D. simulans synthetic lines with compound chromosomes from a population of the Seychelles Islands (high frequency of interspecific mating) and a multimarker strain (low frequency), we show that D. simulans males play an important role in discriminating D. melanogaster females. The genetics of male discrimination fits well with the inheritance mode of a single locus, dominant for sexual isolation, located in chromosome II near the net mutation (2L-0.0). The heterospecific mating success of the male was not related to his sexual vigor. The specific load of male cuticular hydrocarbons was counted as a possible source of discrimination used by the D. melanogaster female.     

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.     

Rapid evolution at the Drosophila telomere: transposable element dynamics at an intrinsically unstable locus

Drosophila telomeres have been maintained by three families of active transposable elements (TEs), HeT-A, TAHRE, and TART, collectively referred to as HTTs, for tens of millions of years, which contrasts with an unusually high degree of HTT interspecific variation. While the impacts of conflict and domestication are often invoked to explain HTT variation, the telomeres are unstable structures such that neutral mutational processes and evolutionary tradeoffs may also drive HTT evolution. We leveraged population genomic data to analyze nearly 10,000 HTT insertions in 85  Drosophila melanogaster genomes and compared their variation to other more typical TE families. We observe that occasional large-scale copy number expansions of both HTTs and other TE families occur, highlighting that the HTTs are, like their feral cousins, typically repressed but primed to take over given the opportunity. However, large expansions of HTTs are not caused by the runaway activity of any particular HTT subfamilies or even associated with telomere-specific TE activity, as might be expected if HTTs are in strong genetic conflict with their hosts. Rather than conflict, we instead suggest that distinctive aspects of HTT copy number variation and sequence diversity largely reflect telomere instability, with HTT insertions being lost at much higher rates than other TEs elsewhere in the genome. We extend previous observations that telomere deletions occur at a high rate, and surprisingly discover that more than one-third do not appear to have been healed with an HTT insertion. We also report that some HTT families may be preferentially activated by the erosion of whole telomeres, implying the existence of HTT-specific host control mechanisms. We further suggest that the persistent telomere localization of HTTs may reflect a highly successful evolutionary strategy that trades away a stable insertion site in order to have reduced impact on the host genome. We propose that HTT evolution is driven by multiple processes, with niche specialization and telomere instability being previously underappreciated and likely predominant.     

Variable post‐zygotic isolation in Drosophila melanogaster/D. simulans hybrids

The study of hybrid inviability reveals cryptic divergence between the genetic interactions that maintain stable phenotypes in the pure species . We characterized the effects of natural variation on the penetrance of hybrid inviability phenotypes in crosses between Drosophila melanogaster and two species of the D. simulans subcomplex, D. simulans and D. sechellia. Using a panel of wild‐caught lines, we studied the levels of genetic variance present in D. simulans and D. sechellia affecting prezygotic and post‐zygotic isolation in hybridizations with D. melanogaster females. We observed extensive variability in the viability of hybrid individuals, dependent on the genotype of the parents, suggesting that intraspecific natural variation manifests directly in hybrid phenotypes. Furthermore, we found that genetic background significantly affects the penetrance of a well‐studied determinant of hybrid inviability: the interaction between Hmr^mel–Lhr^sim. Our results suggest that hybrid inviability – and reproductive isolation generally – can be modified by polymorphisms at multiple loci segregating within the parental species. Just as the penetrance of most mutant phenotypes can be modified by the genetic background within the pure species, the penetrance of hybrid inviability phenotypes is highly influenced by the parental genotypes.