Do parasitoid preferences for different host species match virulence?

Abstract.  Leptopilina boulardi is a parasitoid wasp specialist of Drosophila larvae of the melanogaster subgroup. In Mediterranean areas, natural populations are highly virulent against their main host Drosophila melanogaster . In Congo, populations are less virulent against D. melanogaster but are able to develop successfully inside the tropical African species Drosophila yakuba . Host preferences are compared between two laboratory isofemale lines of L. boulardi , obtained from populations of Congo and Tunisia, respectively, and differing in virulence levels against D. melanogaster and D. yakuba . Host selection is studied by offering female parasitoids a choice between larvae of the two host species. In agreement with optimal foraging models, the line highly virulent against D. melanogaster shows a clear preference for this host species. The other line, less virulent against D. melanogaster but more virulent against D. yakuba , prefers to oviposit on D. yakuba . Such preferences can be observed after a period of host-patch exploitation only, suggesting that experience plays an important role in the host-selection process. These results evidence the existence of intraspecific variability in preference between two host species in L. boulardi , a major requisite in theoretical models of parasite specialization by the host. They also sustain the hypothesis that intraspecific variation in parasitoid preferences between host species might mirror intraspecific variation in virulence.     

High density of long dinucleotide microsatellites in Drosophila subobscura

We isolated 96 dinucleotide repeats with five or more tandemly repeated units from a subgenomic Drosophila subobscura library. The mean repeat unit length of microsatellite clones in D. subobscura is 15, higher than that observed in other Drosophila species. Population variation was assayed in 32-40 chromosomes from Barcelona, Spain, using 18 randomly chosen microsatellite loci. Positive correlation between measures of variation and perfect repeat length measures (mean size, most common, and longest allele) is consistent with a higher mutation rate in loci with longer repeat units. Levels of microsatellite variation measured as variance in repeat number and heterozygosity in D. subobscura were similar to those of Drosophila pseudoobscura and higher than those of Drosophila melanogaster and Drosophila simulans. Our data suggest that higher levels of microsatellite variation, and possibly density, in D. subobscura compared with D. melanogaster are due to both a higher average effective population and a higher intrinsic slippage rate in the former species.     

Rapid divergence of microsatellite abundance among species of Drosophila

Among major taxonomic groups, microsatellites exhibit considerable variation in composition and allele length, but they also show considerable conservation within many major groups. This variation may be explained by slow microsatellite evolution so that all species within a group have similar patterns of variation, or by taxon-specific mutational or selective constraints. Unfortunately, comparing microsatellites across species and studies can be problematic because of biases that may exist among different isolation and analysis protocols. We present microsatellite data from five Drosophila species in the Drosophila subgenus: D. arizonae, D. mojavensis, and D. pachea (three cactophilic species), and D. neotestacea and D. recens (two mycophagous species), all isolated at the same time using identical protocols. For each species, we compared the relative abundance of motifs, the distribution of repeat size, and the average number of repeats. Dimers were the most abundant microsatellites for each species. However, we found considerable variation in the relative abundance of motif size classes among species, even between sister taxa. Frequency differences among motifs within size classes for the three cactophilic species, but not the two mycophagous species, are consistent with other studied Drosophila. Frequency distributions of repeat number, as well as mean size, show significant differences among motif size classes but not across species. Sizes of microsatellites in these five species are consistent with D. virilis, another species in the subgenus Drosophila, but they have consistently higher means than in D. melanogaster, in the subgenus Sophophora. These results confirm that many aspects of microsatellite variation evolve quickly but also are subject to taxon-specific constraints. In addition, the nature of microsatellite evolution is dependent on temporal and taxonomic scales, and some variation is conserved across broad taxonomic levels despite relatively high rates of mutation for these loci.     

Tissue-Specific Expression Phenotypes of Hawaiian Drosophila Adh Genes in Drosophila Melanogaster Transformants

Interspecific differences in the tissue-specific patterns of expression displayed by the alcohol dehydrogenase (Adh) genes within the Hawaiian picture-winged Drosophila represent a rich source of evolutionary variation in gene regulation. Study of the cis-acting elements responsible for regulatory differences between Adh genes from various species is greatly facilitated by analyzing the behavior of the different Adh genes in a homogeneous background. Accordingly, the Adh gene from Drosophila grimshawi was introduced into the germ line of Drosophila melanogaster by means of P element-mediated transformation, and transformants carrying this gene were compared to transformants carrying the Adh genes from Drosophila affinidisjuncta and Drosophila hawaiiensis. The results indicate that the D. affinidisjuncta and D. grimshawi genes have relatively higher levels of expression and broader tissue distribution of expression than the D. hawaiiensis gene in larvae. All three genes are expressed at similar overall levels in adults, with differences in tissue distribution of enzyme activity corresponding to the pattern in the donor species. However, certain systematic differences between Adh gene expression in transformants and in the Hawaiian Drosophila are noted along with tissue-specific position effects in some cases. The implications of these findings for the understanding of evolved regulatory variation are discussed.     

X chromosome DNA variation in Drosophila virilis.

Comparisons of polymorphism patterns between distantly related species are essential in order to determine their generality. However, most work on the genus Drosophila has been done only with species of the subgenus Sophophora. In the present work, we have sequenced one intron and surrounding coding sequences of 6 X-linked genes (chorion protein s36, elav, fused, runt, suppressor of sable and zeste) from 21 strains of wild-type Drosophila virilis (subgenus Drosophila). From these data, we have estimated the average level of DNA polymorphism, inferred the effective population size and population structure of this species, and compared the results with those obtained for other Drosophila species. There is no reduction in variation at two loci close to the centromeric heterochromatin, in contrast to Drosophila melanogaster.     

Effects of extreme temperatures on phenotypic variation and developmental stability in Drosophila melanogaster and Drosophila buzzatii

The effect of stressful and nonstressful rearing temperatures on phenotypic variation of four quantitative characters (thorax length, wing length, number of sternopleural chaetae, number of arista branches) and on developmental stability (fluctuating asymmetry) of the three latter characters was estimated in two Drosophila species: Drosophila melanogaster and Drosophila buzzatii . In both species, a general trend for increasing of phenotypic variation and fluctuating asymmetry at stress temperatures was observed; in fluctuating asymmetry, this effect was more pronounced. An increase of phenotypic variation under stress was shown for all characters examined except sternopleural chaeta number in D. buzzatii . Comparison of species responses suggests that the increase of variation in D. melanogaster was somewhat higher than in D. buzzatii .     

Wolbachia-mediated persistence of mtDNA from a potentially extinct species

Drosophila quinaria is polymorphic for infection with Wolbachia, a maternally transmitted endosymbiont. Wolbachia-infected individuals carry mtDNA that is only distantly related to the mtDNA of uninfected individuals, and the clade encompassing all mtDNA haplotypes within D. quinaria also includes the mtDNA of several other species of Drosophila. Nuclear gene variation reveals no difference between the Wolbachia-infected and uninfected individuals of D. quinaria, indicating that they all belong to the same interbreeding biological species. We suggest that the Wolbachia and the mtDNA with which it is associated were derived via interspecific hybridization and introgression. The sequences in the Wolbachia and the associated mtDNA are ≥6% divergent from those of any known Drosophila species. Thus, in spite of nearly complete species sampling, the sequences from which these mitochondria were derived remain unknown, raising the possibility that the donor species is extinct. The association between Wolbachia infection and mtDNA type within D. quinaria suggests that Wolbachia may be required for the continued persistence of the mtDNA from an otherwise extinct Drosophila species. We hypothesize that pathogen-protective effects conferred by Wolbachia operate in a negative frequency-dependent manner, thus bringing about a stable polymorphism for Wolbachia infection.     

Evolution of the response to heat shock in genus Drosophila

Thermotolerance was studied in a wide spectrum of Drosophila species and strains originating from different climatic zones and considerably differing from one another in the ambient temperature of their habitats. The species that lived in hot climate have a higher thermotolerance. Most species of the virilis group exhibited positive correlation between the HSP70 accumulation after heat exposure and thermotolerance; however, this correlation was absent in some species and strains. For example, the D. melanogaster Oregon R strain, which had the highest sensitivity to heat shock (HS) among all strains and species studied, displayed the maximum level of HSP70 proteins after HS. The patterns of induction of various heat shock protein (HSP) families after heat exposure in a wide spectrum of Drosophila species were compared. The results obtained suggest that the HSP40 and low-molecular-weight HSPs (lmwHSPs) play a significant role in thermotolerance and adaptation to hot climate. Polymorphism in hsp70 gene clusters of Drosophila and variation in the numbers of gene copies and hsp70 isoforms in group virilis were found. The evolutionary role of the variation in the number of hsp70 gene copies observed in the strains and species of genus Drosophila is discussed.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Evolution of arginine kinase within the genus Drosophila

The rate of evolution and extant levels of polymorphism for arginine kinase were examined in species of the genus Drosophila. Surveys of 11 species for electrophoretic variation revealed an average heterozygosity of 0.003. Using antisera prepared against arginine kinase purified from a representative of each of three major subdivisions of the genus, immunological distances were measured by microcomplement fixation (MC'F) assay for 20 species of Drosophila. These data are consistent with the broad outline of the phylogenetic relationships within the genus suggested by other kinds of data. The unit evolutionary period for Drosophila arginine kinase was estimated to be 59 million years (MY). This contrasts with an estimate of 30 MY as the unit evolutionary period for vertebrate creatine kinase. It is suggested that this difference in evolutionary rate arises because a single gene encodes Drosophila arginine kinase, whereas at least three different genes encode the various forms of vertebrate creatine kinase.     

Clinal variation in Drosophila serrata for stress resistance and body size

Clines for size and stress resistance traits have been described for several Drosophila species and replicable clines across different species may indicate climatic selection. Here we consider clines in stress resistance traits in an Australian endemic species, D. serrata, by comparing levels of variation within and among isofemale lines initiated with flies collected from the eastern coast of Australia. We also consider clinical variation in chill coma recovery, a trait that has recently been shown to exhibit high levels of variation among Drosophila species. Patterns were compared with those in the cosmopolitan species D. melanogaster from the same area. Both desiccation and starvation resistance showed no clinical pattern despite heritable variation among isofemale lines. In contrast chill coma resistance exhibited a linear cline in the anticipated direction, resistance increasing with latitude. Body size was measured as wing length and body weight. Both traits showed geographic variation and strong non-linear clines with a sharp reduction in size in the tropics. These results are discussed in the context of climatic selection and evolutionary processes limiting species borders.     

A meristic trait in Drosophila melanogaster x simulans hybrids

A quantitative trait, abdominal bristle number, is described in a number of stocks of the two sibling species Drosophila melanogaster and Drosophila simulans and in interspecific hybrids obtained from these stocks. Coefficients of variation for the trait in the parental species and the hybrids are compared, and in the case of the latter are found to fall in the same range as those of the pure species. The results are discussed in the light of previous findings on other bristle systems.     

Comparative genomic analysis of allatostatin-encoding (Ast) genes in Drosophila species and prediction of regulatory elements by phylogenetic footprinting

The role of the YXFGLa family of allatostatin (AST) peptides in dipterans is not well-established. The recent completion of sequencing of genomes for multiple Drosophila species provides an opportunity to study the evolutionary variation of the allatostatins and to examine regulatory elements that control gene expression. We performed comparative analyses of Ast genes from seven Drosophila species (Drosophila melanogaster, Drosophila simulans, Drosophila ananassae, Drosophila yakuba, Drosophila pseudoobscura, Drosophila mojavensis, and Drosophila grimshawi) and used phylogenetic footprinting methods to identify conserved noncoding motifs, which are candidates for regulatory regions. The peptides encoded by the Ast precursor are nearly identical across species with the exception of AST-1, in which the leading residue may be either methionine or valine. Phylogenetic footprinting predicts as few as 3, to as many as 17 potential regulatory sites depending on the parameters used during analysis. These include a Hunchback motif approximately 1.2 kb upstream of the open reading frame (ORF), overlapping motifs for two Broad-complex isoforms in the first intron, and a CF2-II motif located in the 3'-UTR. Understanding the regulatory elements involved in Ast expression may provide insight into the function of this neuropeptide family.     

Unraveling Selection in the Mitochondrial Genome of Drosophila

We examine mitochondrial DNA variation at the cytochrome b locus within and between three species of Drosophila to determine whether patterns of variation conform to the predictions of neutral molecular evolution. The entire 1137-bp cytochrome b locus was sequenced in 16 lines of Drosophila melanogaster, 18 lines of Drosophila simulans and 13 lines of Drosophila yakuba. Patterns of variation depart from neutrality by several test criteria. Analysis of the evolutionary clock hypothesis shows unequal rates of change along D. simulans lineages. A comparison within and between species of the ratio of amino acid replacement change to synonymous change reveals a relative excess of amino acid replacement polymorphism compared to the neutral prediction, suggestive of slightly deleterious or diversifying selection. There is evidence for excess homozygosity in our world wide sample of D. melanogaster and D. simulans alleles, as well as a reduction in the number of segregating sites in D. simulans, indicative of selective sweeps. Furthermore, a test of neutrality for codon usage shows the direction of mutations at third positions differs among different topological regions of the gene tree. The analyses indicate that molecular variation and evolution of mtDNA are governed by many of the same selective forces that have been shown to govern nuclear genome evolution and suggest caution be taken in the use of mtDNA as a ``neutral' molecular marker.     

Variation of the mitochondrial genome in the evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. The genome contains no introns involved in recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely low and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as the inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophila made it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

Overdispersion of the molecular clock: temporal variation of gene-specific substitution rates in Drosophila

Simple models of molecular evolution assume that sequences evolve by a Poisson process in which nucleotide or amino acid substitutions occur as rare independent events. In these models, the expected ratio of the variance to the mean of substitution counts equals 1, and substitution processes with a ratio greater than 1 are called overdispersed. Comparing the genomes of 10 closely related species of Drosophila, we extend earlier evidence for overdispersion in amino acid replacements as well as in four-fold synonymous substitutions. The observed deviation from the Poisson expectation can be described as a linear function of the rate at which substitutions occur on a phylogeny, which implies that deviations from the Poisson expectation arise from gene-specific temporal variation in substitution rates. Amino acid sequences show greater temporal variation in substitution rates than do four-fold synonymous sequences. Our findings provide a general phenomenological framework for understanding overdispersion in the molecular clock. Also, the presence of substantial variation in gene-specific substitution rates has broad implications for work in phylogeny reconstruction and evolutionary rate estimation.     

An interspecific QTL study of Drosophila wing size and shape variation to investigate the genetic basis of morphological differences

The Drosophila wing has been used as a model in studies of morphogenesis and evolution; the use of such models can contribute to our understanding of mechanisms that promote morphological divergence among populations and species. We mapped quantitative trait loci (QTL) affecting wing size and shape traits using highly inbred introgression lines between D. simulans and D. sechellia, two sibling species of the melanogaster subgroup. Eighteen QTL peaks that are associated with 12 wing traits were identified, including two principal components. The wings of D. simulans and D. sechellia significantly diverged in size; two of the QTL peaks could account for part of this interspecific divergence. Both of these putative QTLs were mapped at the same cytological regions as other QTLs for intraspecific wing size variation identified in D. melanogaster studies. In these regions, one or more loci could account for intra- and interspecific variation in the size of Drosophila wings. Three other QTL peaks were related to a pattern of interspecific variation in wing size and shape traits that is summarized by one principal component. In addition, we observed that female wings are significantly larger and longer than male wings and the second, fourth and fifth longitudinal veins are closer together at the distal wing area. This pattern was summarized by another principal component, for which one QTL was mapped.     

High hemocyte load is associated with increased resistance against parasitoids in Drosophila suzukii, a relative of D. melanogaster

Among the most common parasites of Drosophila in nature are parasitoid wasps, which lay their eggs in fly larvae and pupae. D. melanogaster larvae can mount a cellular immune response against wasp eggs, but female wasps inject venom along with their eggs to block this immune response. Genetic variation in flies for immune resistance against wasps and genetic variation in wasps for virulence against flies largely determines the outcome of any fly-wasp interaction. Interestingly, up to 90% of the variation in fly resistance against wasp parasitism has been linked to a very simple mechanism: flies with increased constitutive blood cell (hemocyte) production are more resistant. However, this relationship has not been tested for Drosophila hosts outside of the melanogaster subgroup, nor has it been tested across a diversity of parasitoid wasp species and strains. We compared hemocyte levels in two fly species from different subgroups, D. melanogaster and D. suzukii, and found that D. suzukii constitutively produces up to five times more hemocytes than D. melanogaster. Using a panel of 24 parasitoid wasp strains representing fifteen species, four families, and multiple virulence strategies, we found that D. suzukii was significantly more resistant to wasp parasitism than D. melanogaster. Thus, our data suggest that the relationship between hemocyte production and wasp resistance is general. However, at least one sympatric wasp species was a highly successful infector of D. suzukii, suggesting specialists can overcome the general resistance afforded to hosts by excessive hemocyte production. Given that D. suzukii is an emerging agricultural pest, identification of the few parasitoid wasps that successfully infect D. suzukii may have value for biocontrol.