Patterns of variation in wing morphology in the cactophilic Drosophila buzzatii and its sibling D. koepferae

Drosophila buzzatii and D. koepferae are two sibling species that breed on the necrotic tissues of several cactus species and show a certain degree of niche overlap. Also, they show differences in several life history traits, such as body size and developmental time, which probably evolved as a consequence of adaptation to different host plants. In this work we investigate the ecological and genetic factors affecting wing morphology variation both within and between species. Three wing traits were scored, distal and proximal wing length and width in isofemale lines reared in two of the most important host cacti: Opuntia sulphurea and Trichocereus terschekii. Our results revealed that differences between species and sexes in wing size and shape were significant, whereas the cactus factor was only significant for wing size. Intraspecific analyses showed that differences among isofemale lines were highly significant for both size and shape in both species, suggesting that an important fraction of variation in wing morphology has a genetic basis. Moreover, the line by cactus interaction, which can be interpreted as a genotype by environment interaction, also accounted for a significant proportion of variation. In summary, our study shows that wing size is phenotypically plastic and that populations of D. buzzatii and D. koepferae harbour substantial amounts of genetic variation for wing size and shape. Interspecific differences in wing size and shape are interpreted in terms of spatial predictability of the different host plants in nature.     

Oviposition preference and life history traits in cactophilic Drosophila koepferae and D. buzzatii in association with their natural hosts

Drosophila koepferae and D. buzzatii are two closely related cactophilic species inhabiting the arid lands of southern South America. Previous studies have shown that D. buzzatii breeds primarily on the necrotic cladodes of several Opuntia cacti and D. koepferae on the rotting stems of columnar cacti of the genera Trichocereus and Cereus. In this paper, we analyze the patterns of host plant utilization in a locality where both Drosophila species are sympatric. Field studies showed an absence of differential attraction of adult flies to the rots of two major host cacti: O. sulphurea and T. terschekii. However, the proportion of D. buzzatii flies emerged from the rotting cladodes of O. sulphurea was significantly higher than in T. terschekii. In laboratory experiments, egg to adult viability in single species cultures varied when both Drosophila species were reared in media prepared with O. sulphurea or T. terschekii. In addition, between-species comparisons of flies emerged from single species cultures showed that D. buzzatii adults were smaller and developed faster than D. koepferae. Furthermore, analysis of flies emerged in mixed species cultures showed differences in oviposition preference and oviposition behavior. We discuss the observed between-species differences and suggest that these traits are the result of adaptation to specific patterns of spatial and temporal predictability of their respective preferred host plants: columnar are less dense and less ephemeral resources, whereas the opuntias are more abundant, and fast rotting cacti. This revised version was published online in July 2006 with corrections to the Cover Date.     

Oviposition acceptance and fecundity schedule in the cactophilic sibling species Drosophila buzzatii and D. koepferae on their natural hosts

We tested for the occurrence of oviposition acceptance for different media prepared with cactus tissues of three alternative cactus hosts: Opuntia sulphurea, O. quimilo and Trichocereus terschekii for 4 consecutive days in lines of two Drosophila buzzatii populations and one population of D. koepferae. Our results showed that the former laid significantly more eggs on both Opuntia cacti than on T. terschekii, whereas D. koepferae preferred T. terschekii. In addition, fecundity schedules differed between species: D. buzzatii laid similar numbers of eggs along the four-day sampling period, whereas D. koepferae showed an oviposition peak on the second day of egg collection on T. terschekii. We suggest that the between-species disparities observed in oviposition acceptance and fecundity schedule may be related to the temporal and spatial predictability of Opuntia versus T. terschekii (cardón) as part of the different adaptive strategies that have evolved after the split of D. koepferae and D. buzzatii from their recent common ancestor. Therefore, the willingness to accept hosts would be an important factor in the habitat selection and in the maintenance of species diversity.     

Comparative molecular population genetics of the Xdh locus in the cactophilic sibling species Drosophila buzzatii and D. koepferae

The Xdh (rosy) gene is one of the best studied in the Drosophila genus from an evolutionary viewpoint. Here we analyze nucleotide variation in a 1875-bp fragment of the second exon of Xdh in Argentinian populations of the cactophilic D. buzzatii and its sibling D. koepferae. The major electrophoretic alleles of D. buzzatii not only lack diagnostic amino acids in the region studied but also differ on average from each other by four to 13 amino acid changes. Our data also suggest that D. buzzatii populations belonging to different phytogeographic regions are not genetically differentiated, whereas D. koepferae exhibits a significant pattern of population structure. The Xdh region studied is twice as polymorphic in D. buzzatii as in D. koepferae. Differences in historical population size or in recombinational environment between species could account for the differences in the level of nucleotide variation. In both species, the Xdh region exhibits a great number of singletons, which significantly departs from the frequency spectrum expected under neutrality for nonsynonymous sites and also for synonymous sites in D. buzzatii. These excesses of singletons could be the signature of a recent population expansion in D. buzzatii, whereas they may be simply explained as the result of negative selection in D. koepferae.     

Wing morphology and fluctuating asymmetry depend on the host plant in cactophilic Drosophila

As in most insect groups, host plant shifts in cactophilic Drosophila represent environmental challenges as flies must adjust their developmental programme to the presence of different chemical compounds and/or to a microflora that may differ in the diversity and abundance of yeasts and bacteria. In this context, wing morphology provides an excellent opportunity to investigate the factors that may induce changes during development. In this work, we investigated phenotypic plasticity and developmental instability of wing morphology in flies on the cactophilic Drosophila buzzatii and Drosophila koepferae raised on alternative breeding substrates. We detected significant differences in wing size between and within species, and between flies reared on different cactus hosts. However, differences in wing shape between flies emerged from different cactus hosts were not significant either in D. buzzatii or in D. koepferae. Our results also showed that morphological responses involved the entire organ, as variation in size and shape correlated between different portions of the wing. Finally, we studied the effect of the rearing cactus host on developmental instability as measured by the degree of fluctuating asymmetry (FA). Levels of FA in wing size were significantly greater in flies of both species reared in non-preferred when compared with those reared in preferred host cacti. Our results are discussed in the framework of an integrative view aimed at investigating the relevance of host plant shifts in the evolution of the guild of cactophilic Drosophila species that diversified in South America.     

Host-related life history traits in interspecific hybrids of cactophilic Drosophila

In the genus Drosophila (Diptera: Drosophilidae), interspecific hybridization is a rare phenomenon. However, recent evidence suggests a certain degree of introgression between the cactophilic siblings Drosophila buzzatii Patterson & Wheeler and Drosophila koepferae Fontdevila & Wasserman. In this article, we analyzed larval viability and developmental time of hybrids between males of D. buzzatii and females of D. koepferae, raised in media prepared with fermenting tissues of natural host plants that these species utilize in nature as breeding sites. In all cases, developmental time and larval viability in hybrids was not significantly different from parental lines and, depending on the cross, hybrids developed faster than both parental species or than the slowest species. When data of wing length were included in a discriminant function analysis, we observed that both species can be clearly differentiated, while hybrids fell in two categories, one intermediate between parental species and the other consisting of extreme phenotypes. Thus, our results point out that hybrid fitness, as measured by developmental time and viability, is not lower than in the parental species.     

Intraspecific variation ofDrosophila buzzatii larval breeding success onOpuntia: A yeast-plant-insect relationship

The cactophilic species,Drosophila buzzatii, normally breeds in decaying pockets ofOpuntia cladodes, in which there is a complex interaction with the microbial flora, especially yeast species. Isofemale lines were used to estimate genetic variation among larvae reared on their natural feeding substrate. Four naturally occurring cactophilic yeast species isolated from the same Tunisian oasis as theDrosophila population were used. Two fitness components were studied for each line, viability and developmental time. Genetic variations amongD. buzzatii lines were observed for both traits. A significant yeast species x isofemale line interaction for viability was also evidenced, suggesting the occurrence of specialized genotypes for the utilization of breeding substrates. This genetic heterogeneity in the natural population may favor a better adaptation to the patchily distribution of yeasts.     

GENETICS OF HOST-CACTUS RESPONSE AND LIFE-HISTORY EVOLUTION AMONG ANCESTRAL AND DERIVED POPULATIONS OF CACTOPHILIC DROSOPHILA MOJAVENSIS

The extent of host-specific genetic variation for two life-history traits, egg to adult developmental time and viability, and one morphological trait closely tied to fitness, adult thorax size, was exposed by employing a nested half-sib/full-sib breeding design with Baja and mainland populations of Drosophila mojavensis recently extracted from nature. This study was motivated by the presence of substantial variation in life histories among populations of D. mojavensis that use the fermenting tissues of particular species of columnar cacti for feeding and breeding in the Sonoran Desert. Full-sib progeny from all sire-dam crosses were split into cultures of agria cactus, Stenocereus gummosus, and organ pipe cactus, S. thurberi, to examine patterns of genotype-by-environment interaction for these fitness components. Baja flies expressed shorter egg-to-adult developmental times, higher viabilities, and smaller body sizes than mainland flies consistent with previous studies. Significant sire and dam components of variance were exposed for developmental time and thorax size. Genotype-by-environment interactions were significant at the level of dams for developmental time and nearly significant for viability (P = 0.09). Narrow- and broad-sense heritabilities were influenced by host cactus, sex, and population. No strong pattern of genetic correlation emerged among fitness components suggesting that host-range expansion has not been accompanied by formation of coadapted life histories, yet the ability to estimate genetic correlations and their standard errors was compromised by the unbalanced nature of the data set. Genetic correlations in performance across cacti were slightly positive, evidence for ecological generalism among populations explaining the observed pattern of multiple host cactus use within the species range of D. mojavensis.     

Direct and correlated responses to artificial selection on developmental time and wing length in Drosophila buzzatii

Abstract.— Developmental time and body size are two positively correlated traits closely related to fitness in many organisms including Drosophila . Previous work suggested that these two traits are involved in a trade-off that may result from a negative genetic correlation between their effects on pre-adult and adult fitness. Here, we examine the evolution of developmental time and body size (indexed by wing length) under artificial selection applied to one or both traits in replicated D. buzzatii populations. Directional changes in both developmental time and wing length indicate the presence of substantial additive genetic variance for both traits. The strongest response to selection for fast development was found in lines selected simultaneously to reduce both developmental time and wing length, probably as an expected consequence of a synergistic effect of indirect selection. When selection was applied in the direction opposite to the putative genetic correlation, that is, large wing length but fast development, no responses were observed for developmental time. Lines selected to reduce both wing length and developmental time diverged slightly faster from the control than lines selected to increase wing length and reduce developmental time. However, wing length did not diverge from the control in lines selected only for fast development. These results suggest a complex genetic basis of the correlation between developmental time and wing length, but are generally consistent with the hypothesis that both traits are related in a trade-off. However, we found that this trade-off may disappear under uncrowded conditions, with fast-developing lines exhibiting a higher pre-adult viability than other lines when tested at high larval density.     

Genetic Variation in Foraging Traits and Life-History Traits of the Predatory Mite Neoseiulus womersleyi (Acari: Phytoseiidae) among Isofemale Lines

The predatory mite Neoseiulus womersleyi shows a significant correlation between its olfactory response and dispersal tendency in different geographical populations. This study investigated the genetic background of the relationship using isofemale lines. Y-tube olfactometer tests confirmed that there was a genetic component in predator response to herbivore-induced plant volatiles. Wind tunnel tests in the absence of the herbivore-induced plant volatiles revealed that the dispersal tendencies of N. womersleyi exhibited genetic variation among isofemale lines, and other experiments revealed the existence of significant differences in prey consumption rate, fecundity, and developmental time. However, there was no genetic correlation between behavioral traits (olfactory response, innate dispersal) and the other traits, suggesting that the positive correlation between the behavioral traits was not caused by genetic factors.     

Host shifts and the beginning of signal divergence

Divergence between populations adapting to different environments may be facilitated when the populations differ in their sexual traits. We tested whether colonizing a novel environment may, through phenotypic plasticity, change sexual traits in a way that could alter the dynamics of sexual selection. This hypothesis has two components: changes in mean phenotypes across environments, and changes in the genetic background of the phenotypes that are produced -- or genotype x environment interaction (G x E). We simulated colonization of a novel environment and tested its effect on the mating signals of a member of the Enchenopa binotata species complex of treehoppers (Hemiptera: Membracidae), a clade that has diverged in a process involving host plant shifts and signal diversification. We found substantial genetic variation and G x E in most signal traits measured, with little or no change in mean signal phenotypes. We suggest that the expression of extant genetic variation across old and novel environments can initiate signal divergence.     

Developmental time and size-related traits in Drosophila buzzatii along an altitudinal gradient from Argentina

Clinal analysis for fitness-related traits provides a well-known approach to investigate adaptive evolution. Several fitness-related traits (developmental time, thorax length, wing length and wing loading) were measured at two laboratory generations (G7 and G33) of D. buzzatii from an altitudinal gradient from northwestern Argentina, where significant thermal differences persist. Developmental time (DT) was positively correlated with altitude of origin of population. Further, DT was negatively correlated with maximal mean temperature at the site of origin of population, and this thermal variable decreases with altitude. Wing loading tended to be larger in highland than in lowland populations, suggesting that flight performance is subject to stronger selection pressure in highland populations. Developmental time showed a significant increase with laboratory generation number. There was no significant correlation between developmental time and body size across populations along the altitudinal cline of DT. This result illustrates that developmental time and body size do not always evolve in the same direction, even though both traits are often positively and genetically correlated in a well-known tradeoff in Drosophila.     

LIMITS TO THE SOUTHERN BORDER OF DROSOPHILA SERRATA: COLD RESISTANCE,HERITABLE VARIATION,AND TRADE-OFFS

There are a number of evolutionary hypotheses about why species distributions are limited, but very little empirical information to test them. We present data examining whether the southern distribution of Drosophila serrata is limited by cold responses. Species comparisons were undertaken for cold resistance, development time, and viability at 15°C and 25°C for D. serrata and other species with a more southerly distribution (D. melanogaster, D. simulans, and D. immigrans). Relative to the other species, D. serrata had a long development time at both temperatures and a low level of cold resistance. Using isofemale lines collected in different seasons, central and marginal populations were compared for cold resistance, as well as development time and viability at 14°C. The border population had a relatively higher resistance to cold shock in postwinter collections, but there was no population differentiation for prewinter collections or for the other traits. The presence of variation among isofemale lines within the border populations suggests that genetic variation as measured in the laboratory is unlikely to limit range expansion. Population cages were used in the field to determine if D. serrata persisted over winter at borders. Although all cages yielded adult offspring at northern sites, only a few produced offspring at or just south of the border. In contrast, all cages with D. simulans produced adult offspring, suggesting that climatic factors limited D. serrata numbers. Offspring from surviving adults showed a phenotypic trade-off between fecundity and cold resistance. Comparisons of the cold resistance of field males and females with their laboratory-reared offspring provided evidence for heritable variation in field-reared flies. Overall, the results suggest that cold stress is important in limiting the southern distribution of D. serrata, but it seems unlikely that a lack of genetic variation restricts range expansion.     

The evolutionary history of Drosophila buzzatii. XXIII. High content of nonsatellite repetitive DNA in D. buzzatii and in its sibling D. koepferae.

The frequency and types of repetitive nonsatellite DNA of two sibling species of the repleta group of Drosophila, D. buzzatii, and D. koepferae have been determined. For each species, the analysis is based on a sample of more than 100 clones (400 kb) obtained from genomic DNA. A theoretical model has been developed to correct for the presence of a mixture of repetitive and unique DNA in these clones. After correction, a high content of repetitive DNA has been demonstrated for both species (D. buzzatii, 19-26%; D. koepferae, 27-32%). The repetitive sequences have been classified according to their hybridization pattern when used as probes against genomic DNA and by their in situ hybridization signals on polytene chromosomes. Data suggest that the main nonsatellite component of these species is simpler and more repetitive than that of D. melanogaster, pointing to a wide variability in content and class size distribution of repetitive DNA among Drosophila species.     

Environmental heterogeneity and the maintenance of genetic variation for reproductive diapause in Drosophila melanogaster

Drosophila melanogaster has colonized temperate habitats on multiple continents over a historical time period, and many traits vary predictably with latitude. Despite considerable attention paid to clinal variation in Drosophila, the mechanisms generating such patterns in nature remain largely unidentified. In D. melanogaster, the expression of reproductive diapause can be induced by exposure to low temperatures and shortened photoperiods. Both diapause expression itself and the underlying genetic variance for diapause expression have widespread impacts on organismal fitness, and diapause incidence exhibits a 60% cline in frequency in the eastern United States. The major aim of this study was to evaluate whether the relative fitness of diapause and nondiapause genotypes varies predictably with environment. In experimental population cages in the laboratory, the frequency of genotypes that express diapause increased over time when flies were exposed to environmental stress, whereas the frequency of nondiapause genotypes increased when flies were cultured under benign control conditions. Other fitness traits correlated with the genetic variance for diapause expression (longevity, mortality rates, stress resistance, lipid content, preadult viability, fecundity profiles, and development time) also diverged between experimental treatments. Similarly, sampling of isofemale lines from natural populations revealed that the frequency of diapause incidence cycled over time in seasonal habitats: diapause expression was at high frequency following the winter season and subsequently declined throughout the summer months. In contrast, diapause expression was low and temporally homogeneous in isofemale line collections from human-associated urban habitats. These data suggest that genetic variation underlying the diapause-nondiapause dichotomy may be actively maintained by selection pressures that vary spatially and temporally in natural populations.     

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.     

The role of the use of different host plants in the maintenance of the inversion polymorphism in the cactophilic Drosophila buzzatii

Abstract. Inversion polymorphisms often have been associated with fitness variation. Cactophilic Drosophila buzzatii has been used widely for the study of the maintenance of chromosomal variation. The purpose of this paper is to address the relative importance of variable selection regimes associated with the use of three different host cacti and antagonistic pleiotropy in the maintenance of chromosomal variation. Using homokaryotypic stocks derived from several lines homozygous for four second-chromosome arrangements, we show that inversions significantly affect first-instar larva to adult viability (VT), developmental time (DT) and adult thorax length (TL). We also show that the effects of inversions on DT and VT are dependent on the cactus rearing media. The effects of polymorphic gene arrangements on life-history traits suggest the existence of trade-offs between early and late fitness components. The dosage of arrangement 2st , the ancestral gene order, was negatively correlated with DT and TL, whereas flies carrying the derived arrangements 2j and 2jq⁷ had longer DTs and larger TLs. Arrangements 2st and 2jq⁷ increased viability, at least in one of the cactus media tested. Our results suggest that environmental heterogeneity, as represented by the use of different cactus hosts and the trade-off between DT and TL, may be involved in the maintenance of the polymorphism. In addition, our data suggest that the chromosomal phylogeny may be decoupled from the evolution of the genes affecting life-history traits linked to the inversion system.     

Transpecific polymorphisms in an inversion linked esterase locus in Drosophila buzzatii

Nucleotide variation was studied in a 1.1 kb section of the coding region of an Esterase gene (Est-A) that maps in the center of the segments rearranged by polymorphic inversions in the cactophilic Drosophila buzzatii. We examine 30 homozygous second-chromosome lines differing in gene arrangement and three D. koepferae isofemale lines as outgroups. Our data show that Est-A is a highly polymorphic gene at both synonymous and replacement sites. Significant departures from homogeneity in the distribution of the ratio of silent polymorphism to divergence predicted by the neutral theory reveals a local excess of silent polymorphism. This is consistent with the presence of two apparent narrow peaks of elevated silent polymorphism surrounding nonconservative amino acid substitutions. These polymorphisms as well as others at synonymous and nonsynonymous sites are shared with D. koepferae. We suggest that the presence of shared nucleotide polymorphisms is probably due to interspecific gene flow and/or balancing selection acting on replacement variants and/or to a decreased probability of loss of ancestral polymorphisms caused by linkage to an adaptive inversion polymorphism. Recurrent mutation and persistence of neutral ancestral polymorphisms cannot, however, be ruled out. The analysis of the distribution of nucleotide variation among the three chromosomal arrangements sampled reveals that derived arrangements (J and JZ(3)) are less polymorphic than the ancestral ST, and that the widely distributed ST and J arrangements are genetically differentiated. However, a significant number of polymorphisms are shared between arrangements, suggesting frequent exchange either from gene conversion or from double crossovers in heterokaryotypes. Finally, our present results in combination with data of sequence variation at the breakpoints of inversion J suggest that this old gene arrangement has risen in frequency in relatively recent times.     

A system of complementary genes in hybrids between Drosophila koepferae and D. buzzatii: a Markov chain model allows to make inferences about their number and relationships

In backcrosses between D. koepferae and D. buzzatii, the disruption of a system of species-specific complementary factors brings about hybrid male inviability. This system consists of a lethal factor, hmi-1, linked to the X chromosome of D. koepferae, and several conspecific autosomal suppressors. However, hmi-1 hybrid males can also be rescued by factors present in some strains of D. buzzatii. The present work aims to estimate the number of hmi-1 suppressors in one of these strains by means of Markov chains. The obtained results allow discarding models with one or more chromosomes having independent suppressor effect. On the other hand, models having n chromosomes that interact in groups of r, being 1 < r < or = n, to produce rescue effect, provide good approximations to the observed results. The best fit to the data is obtained with four or five chromosomes with suppressor effect, interacting epistatically in groups of three to rescue the viability of hmi-1 males.     

Genotype?��?environment interaction is weaker in genitalia than in mating signals and body traits in Enchenopa treehoppers (Hemiptera: Membracidae)

Theory predicts that selection acting across environments should erode genetic variation in reaction norms; i.e., selection should weaken genotype × environment interaction (G × E). In spite of this expectation, G × E is often detected in fitness-related traits. It thus appears that G × E is at least sometimes sustained under selection, a possibility that highlights the need for theory that can account for variation in the presence and strength of G × E. We tested the hypothesis that trait differences in developmental architecture contribute to variation in the expression of G × E. Specifically, we assessed the influence of canalization (robustness to genetic or environmental perturbations) and condition-dependence (association between trait expression and prior resource acquisition or vital cellular processes). We compared G × E across three trait types expected to differ in canalization and condition-dependence: mating signals, body size-related traits, and genitalia. Because genitalia are expected to show the least condition-dependence and the most canalization, they should express weaker G × E than the other trait types. Our study species was a member of the Enchenopa binotata species complex of treehoppers. We found significant G × E in most traits; G × E was strongest in signals and body traits, and weakest in genitalia. These results support the hypothesis that trait differences in developmental architecture (canalization and condition-dependence) contribute to variation in the expression of G × E. We discuss implications for the dynamics of sexual selection on different trait types.