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.     

Testing chromosomal phylogenies and inversion breakpoint reuse in Drosophila

A combination of cytogenetic and bioinformatic procedures was used to test the chromosomal phylogeny relating Drosophila buzzatii with D. repleta. Chromosomes X and 2, harboring most of the inversions fixed between these two species, were analyzed. First, chromosomal segments conserved during the divergence of the two species were identified by comparative in situ hybridization to the D. repleta chromosomes of 180 BAC clones from a BAC-based physical map of the D. buzzatii genome. These conserved segments were precisely delimited with the aid of clones containing inversion breakpoints. Then GRIMM software was used to estimate the minimum number of rearrangements necessary to transform one genome into the other and identify all possible rearrangement scenarios. Finally, the most plausible inversion trajectory was tested by hybridizing 12 breakpoint-bearing BAC clones to the chromosomes of seven other species in the repleta group. The results show that chromosomes X and 2 of D. buzzatii and D. repleta differ by 12 paracentric inversions. Nine of them are fixed in chromosome 2 and entail two breakpoint reuses. Our results also show that the cytological relationship between D. repleta and D. mercatorum is closer than that between D. repleta and D. peninsularis, and we propose that the phylogenetic relationships in this lineage of the repleta group be reconsidered. We also estimated the rate of rearrangement between D. repleta and D. buzzatii and conclude that rates within the genus Drosophila vary substantially between lineages, even within a single species group.     

Local adaptation of stress related traits in Drosophila buzzatii and Drosophila simulans in spite of high gene flow

We addressed the question if local adaptation to a thermal gradient is possible in spite of a high gene flow among closely spaced populations of two species of Drosophila from the island of La Gomera (Canary Islands). Variation in multiple traits related to stress resistance in different life stages was measured in both species in flies collected from five localities at different altitudes and thereby with different climatic conditions. Based on microsatellite loci, the populations were not genetically differentiated. However, 18 of the 24 independent traits measured showed significant differentiation among populations of Drosophila buzzatii, but only nine of 25 for Drosophila simulans. This difference in the number of traits might reflect higher habitat specificity and thus higher potential for local adaptation of D. buzzatii than D. simulans. We found clinal variation, as some traits showed significant linear regressions on altitude, but more on altitude cubed.     

On the pBuM189 satellite DNA variability among South American populations of Drosophila buzzatii

The pBuM189 satellite DNA was analysed in Drosophila buzzatii populations that cover most of the species distribution in South America. This satDNA consists of A + T-rich monomers of 189 bp and previous studies showed a fast rate of evolutionary change of this component of D. buzzatii genome. A total of 63 pBuM189 repetition units from 14 D. buzzatii populations (9 from Brazil and 5 from Argentina) were studied. The average nucleotide variability among the 63 repeats is 4.2 %. At least one repeat (Juan/4) seems to be part of another pBuM189 satDNA subfamily. The nucleotide alignment of all 63 repeats revealed no specific nucleotide substitutions, or indels, that could discriminate each population or groups of geographically close populations. Such lack of satDNA interpopulational differentiation is congruent with previous mtDNA data that indicate a high gene flow and very little population differentiation throughout most of the D. buzzatii distribution in South America. Gene flow might have been possible during glaciation events in the Pleistocene, such as the one occurred between 13,000 and 18,000 years ago, when D. buzzatii probably had a more continuous distribution than what is observed today.     

Remating and sperm displacement in a natural population of Drosophila buzzatii inferred from mother-offspring analysis of microsatellite loci

Prospects for estimation of parameters of models of sperm competition from field data have improved recently with the development of methods that employ multilocus genotype data from brood-structured samples. Sperm competition in Drosophila buzzatii is of special interest because it is possible to directly observe the breeding behaviour of this species in its natural habitat of rotting cactus. Previous laboratory experiments showed that this species exhibits an unusual pattern of frequent remating and sperm partitioning. This paper reports the first attempt to estimate the frequency of female remating and sperm competition in natural populations of D. buzzatii. For the Australian population studied, the mean remating frequency was lower (alpha = 2.12-2.20) than previously estimated in laboratory experiments with the same population, whereas mean sperm displacement (beta = 0.69-0.71) fell within the limits of previous laboratory results. The evolution of the D. buzzatii mating system is discussed.     

Chromosomal inversions effect body size and shape in different breeding resources in Drosophila buzzatii

The cactophilic Drosophila buzzatii provides an excellent model for the study of reaction norms across discrete environments because it breeds on rotting tissues (rots) of very different cactus species. Here we test the possible effects of second chromosome inversions on body size and shape (wing loading) across suitable natural breeding substrates. Using homokaryotypic stocks derived from several lines homozygous for four naturally occurring chromosomal inversions, we show that arrangements significantly affect size-related traits and wing loading. In addition, karyotypes show differing effects, across natural breeding resources, for wing loading. The 2st and 2jz(3) arrangements decrease and the 2j arrangement increases wing loading. For thorax length and wing loading, karyotypic correlations across host plants are slightly lower in females than in males. These results support the hypothesis that these traits have a genetic basis associated with the inversion polymorphism.     

Gametic associations between inversion and allozyme polymorphisms in Drosophila buzzatii.

Gametic disequilibria between second chromosome polymorphic arrangements and seven linked allozyme loci were estimated in seven populations of Drosophila buzzatii from Argentina. Significant and consistent associations across populations were detected for Est-1, Est-2, Aldox, and XDH: Phenograms based on Nei's genetic distance showed that chromosomes carrying the 2ST arrangement were more similar to each other, irrespective of the population from which they were extracted, than to chromosomes carrying the derived 2J and 2JZ3.Restriction of recombination in heterokaryotypes seems to be the best explanation for the significant linkage disequilibria between inversions and the loci located inside the rearranged segments, for example, Est-1 and Aldox, or close to the break points, for example, Est-2. However, epistatic interactions between Xdh, which is outside the inversions and not near the break points, and loci tightly linked to the inversions, is the most likely explanation for the association between Xdh and chromosomal arrangements. Some of the associations detected in endemic Argentinean populations are coincident with data obtained in colonizing populations of the Old World and Australia. Thus historical processes that took place in the original area of the species' distribution can account for these linkage disequilibria in colonized populations of D. buzzatii.     

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.     

Inversion frequencies in Drosophila serrata along an eastern Australian transect.

Clinal patterns over broad geographic regions provide a way of identifying characteristics of species under selection and are increasingly being used in quantitative trait locus mapping of adaptive genetic variation in Drosophila. However, interpretations of clinal patterns can be complicated by inversions that also vary clinally and reduce recombination in some parts of the genome. Drosophila serrata (Malloch) is an Australian endemic species being used to investigate the genetic basis of geographic variation in climatic adaptation and mate recognition. Here we describe inversions in D. serrata populations from the east coast of Australia, covering tropical and temperate regions. Seven autosomal paracentric inversions and 1 apparently complex X chromosome arrangement were identified from these populations. All inverted arrangements were relatively more common in tropical populations; 2 common inversions showed clinal patterns over part of the range of D. serrata. Inversion polymorphism was relatively higher in tropical populations and almost absent in populations near the cooler southern border, in agreement with findings on other Drosophila species. While these patterns will complicate mapping of adaptive variation in D. serrata, they suggest that this species will be useful in investigatingthe dynamics of inversion-trait associations in natural populations.     

The Evolutionary History of DROSOPHILA BUZZATII. III. Cytogenetic Relationships between Two Sibling Species of the Buzzatii Cluster

Drosophila buzzatii has been found sympatric in Argentina with a closely-related sibling species, D. serido. The biogeographical, reproductive and chromosomal data allow us to combine these species into an evolutionary unit, the buzzatii cluster. Salivary gland chromosomes also have been used to determine their phylogenetic relationships with other closely related species, showing that the buzzatii cluster species share two inversions—2d² and 2s⁶—with the species of the martensis cluster. Both clusters arose from South American populations of the ancestor of the mulleri complex, and we propose to include D. buzzatii and D. serido in the mulleri complex of the repleta group.     

The foldback-like transposon Galileo is involved in the generation of two different natural chromosomal inversions of Drosophila buzzatii

Chromosomal inversions are the most common type of genome rearrangement in the genus Drosophila. Although the potential of transposable elements (TEs) for generating inversions has been repeatedly demonstrated in the laboratory, little is known on their role in the generation of natural inversions, which are those effectively contributing to the adaptation and/or evolution of species. We have cloned and sequenced the two breakpoints of the polymorphic inversion 2q7 of D. buzzatii. The sequence analysis of the breakpoint regions revealed the presence in the inverted chromosomes of large insertions, formed by complex assemblies of transposons, that are absent from the chromosomes without the inversion. Among the transposons inserted, the Foldback-like element Galileo, that was previously found responsible of the generation of the widespread inversion 2j of D. buzzatii, is present at both 2q7 breakpoints and is the most likely inducer of the inversion. A detailed study of the nucleotide and structural variation in the breakpoint regions of six chromosomal lines with the 2q7 inversion detected no nucleotide differences between them, which suggests a monophyletic and recent origin. In contrast, a remarkable degree of structural variation was observed in the same six chromosomal lines. It thus appears that the two breakpoints of the inverted chromosomes have become genetically unstable hotspots, as was previously found for the 2j inversion breakpoints. The possibility that this instability is caused by structural properties of Foldback elements is discussed.     

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 .     

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.     

Characterization of a species-specific satellite DNA of Drosophila buzzatii

The pBuM189 satellite DNA family was found to be species specific for Drosophila buzzatii . It consists of slightly AT-rich tandemly arranged repeats with a high copy number in the genome and shows a very high level of intraspecific sequence similarity. pBuM189 repeats cannot be detected in the genomes of closely related species such as Drosophila serido , Drosophila borborema and Drosophila koepferae . The data support the marginal systematic position of D. buzzatii within the buzzatii cluster of the Drosophila repleta group.     

Genetic mapping of the Adh locus in the repleta group of Drosophila by in situ hybridization

A biotinylated probe of the Adh (alcohol dehydrogenase) gene of Drosophila melanogaster was used for in situ hybridization on polytene chromosomes of D. mojavensis and D. buzzatii, two species of the repleta group of the genus Drosophila. Hybridization showed that the Adh gene maps at the G1a band of the third chromosome. This is in accordance with a previous result obtained through the use of interspecific hybrid asynapsis as a cytological marker and establishes the limits of the precision of this method.     

Genetic Architecture of Autosome-Mediated Hybrid Male Sterility in Drosophila

Several estimators have been developed for assesing the number of sterility factors in a chromosome based on the sizes of fertile and sterile introgressed fragments. Assuming that two factors are required for producing sterility, simulations show that one of these, twice the inverse of the relative size of the largest fertile fragment, provides good average approximations when as few as five fertile fragments are analyzed. The estimators have been used for deducing the number of factors from previous data on several pairs of species. A particular result contrasts with the authors' interpretations: instead of the high number of sterility factors suggested, only a few per autosome are estimated in both reciprocal crosses involving Drosophila buzzatii and D. koepferae. It has been possible to map these factors, between three and six per chromosome, in the autosomes 3 and 4 of these species. Out of 203 introgressions of different fragments or combinations of fragments, the outcome of at least 192 is explained by the mapped zones. These results suggest that autosome-mediated sterility in the male hybrids of these species is mediated by a few epistatic factors, similarly to X-mediated sterility in the hybrids of other Drosophila species.     

The evolutionary history of Drosophila buzzatii. XXXV. Inversion polymorphism and nucleotide variability in different regions of the second chromosome

Inversions are portions of a chromosome where the gene order is reversed relative to a standard reference orientation. Because of reduced levels of recombination in heterokaryotypes, inversions have a potentially important effect on patterns of nucleotide variability in those genomic regions close to, or included in, the inverted fragments. Here we report sequence variation at three anonymous regions (STSs) located at different positions in relation to second-chromosome inversion breakpoints in 29 isochromosomal lines derived from an Argentinean population of Drosophila buzzatii. In agreement with previous findings in Drosophila, gene flux (crossing over and/or gene conversion) between arrangements seems to appreciably increase as we approach the middle sections of inversion 2j, and patterns of nucleotide variability within, as well as genetic differentiation between chromosome arrangements, are comparable to those observed at the molecular marker outside the inverted fragments. On the other hand, nucleotide diversity near the proximal breakpoint of inversion 2j is reduced when contrasted with that found at the other regions, particularly in the case of derived inverted chromosomes. Using the data from the breakpoint, we estimate that the inversion polymorphism is approximately 1.63 N generations old, where N is the effective population size. An excess of low-frequency segregating polymorphisms is detected; mostly in the ancestral 2st arrangement and probably indicating a population expansion that predates the coalescent time of inversion 2j. Heterogeneity in mutation rates between the markers linked to the inversions may be sufficient to explain the different levels of nucleotide diversity observed. When considered in the context of other studies on patterns of variation relative to physical distance to inversion breakpoints, our data appear to be consistent with the conclusion that inversions are unlikely to be "long-lived" balanced polymorphisms.     

Toward a physical map of Drosophila buzzatii. Use of randomly amplified polymorphic dna polymorphisms and sequence-tagged site landmarks

We present a physical map based on RAPD polymorphic fragments and sequence-tagged sites (STSs) for the repleta group species Drosophila buzzatii. One hundred forty-four RAPD markers have been used as probes for in situ hybridization to the polytene chromosomes, and positive results allowing the precise localization of 108 RAPDs were obtained. Of these, 73 behave as effectively unique markers for physical map construction, and in 9 additional cases the probes gave two hybridization signals, each on a different chromosome. Most markers (68%) are located on chromosomes 2 and 4, which partially agree with previous estimates on the distribution of genetic variation over chromosomes. One RAPD maps close to the proximal breakpoint of inversion 2z(3) but is not included within the inverted fragment. However, it was possible to conclude from this RAPD that the distal breakpoint of 2z(3) had previously been wrongly assigned. A total of 39 cytologically mapped RAPDs were converted to STSs and yielded an aggregate sequence of 28,431 bp. Thirty-six RAPDs (25%) did not produce any detectable hybridization signal, and we obtained the DNA sequence from three of them. Further prospects toward obtaining a more developed genetic map than the one currently available for D. buzzatii are discussed.     

Chromosomal Homology and Molecular Organization of Muller''s Elements D and E in the Drosophila Repleta Species Group

Thirty-three DNA clones containing protein-coding genes have been used for in situ hybridization to the polytene chromosomes of two Drosophila repleta group species, D. repleta and D. buzzatii. Twenty-six clones gave positive results allowing the precise localization of 26 genes and the tentative identification of another nine. The results were fully consistent with the currently accepted chromosomal homologies and in no case was evidence for reciprocal translocations or pericentric inversions found. Most of the genes mapped to chromosomes 2 and 4 that are homologous, respectively, to chromosome arms 3R and 3L of D. melanogaster (Muller's elements E and D). The comparison of the molecular organization of these two elements between D. melanogaster and D. repleta (two species that belong to different subgenera and diverged some 62 million years ago) showed an extensive reorganization via paracentric inversions. Using a maximum likelihood procedure, we estimated that 130 paracentric inversions have become fixed in element E after the divergence of the two lineages. Therefore, the evolution rate for element E is approximately one inversion per million years. This value is comparable to previous estimates of the rate of evolution of chromosome X and yields an estimate of 4.5 inversions per million years for the whole Drosophila genome.