Phylogenetic utility of the mitochondrial cytochrome oxidase gene: Molecular evolution of the Drosophila buzzatii species complex

Phylogenetic relationships among eight species of the Drosophila buzzatii species complex (D. mulleri subgroup; D. repleta species group) and D. hamatofila were determined by sequencing the mitochondrial cytochrome oxidase subunit 1, II, and III genes. The species examined included members of the martensis cluster (D. martensis, D. starmeri, D. venezolana), the buzzatii cluster (D. buzzatii, D. serido, D. borborema), and the stalkeri cluster (D. stalkeri, D. richardsoni). The molecular phylogeny was found to be congruent with the chromosomal inversion phylogeny. Analyzing the cytochrome oxidase subunits separately revealed that not all the subunits seem to have the same phylogenetic information content. Parameters are discussed that might explain these differences. Correspondence to: G. Spicer     

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.     

Intra‐ and interspecific divergence in the nuclear sequences of the clock gene period in species of the Drosophila buzzatii cluster

We have measured nucleotide variation in the CLOCK/CYCLE heterodimer inhibition domain (CCID) of the clock X‐linked gene period in seven species belonging to the Drosophila buzzatii cluster, namely D. buzzatii, Drosophila koepferae, Drosophila antonietae, Drosophila serido, Drosophila gouveai, Drosophila seriema and Drosophila borborema. We detected that the purifying selection is the main force driving the sequence evolution in period, in agreement with the important role of CCID in clock machinery. Our survey revealed that period provides valuable phylogenetic information that allowed to resolve phylogenetic relationships among D. gouveai, D. borborema and D. seriema, which composed a polytomic clade in preliminary studies. The analysis of patterns of intraspecific variation revealed two different lineages of period in D. koepferae, probably reflecting introgressive hybridization from D. buzzatii, in concordance with previous molecular data.     

Reaction norms across and genetic parameters at different temperatures for thorax and wing size traits in Drosophila aldrichi and D. buzzatii

Reaction norms across seven constant and one fluctuating temperature of development were measured for thorax length and several wing size traits for up to 10 isofemale lines of each of the cactophilic Drosophila species, D. aldrichi and D. buzzatii, originating from the same locality. Maximum thorax length was reached at different low to intermediate temperatures for the two species, whereas wing length was highest at the lowest temperature in both species. Various ratio parameters showed pronounced species differences. The reaction norm for the wing loading index (wing length/thorax length) decreased monotonically with temperature in both species, but was much steeper and spanned a wider range in D. aldrichi than in D. buzzatii, suggesting either that wing loading is not a good characterization of flight capacity or, more likely, that flight optimization does not occur in the same manner in both species. The vein ratio (distal length/proximal length of the third vein) increased with temperature in D. buzzatii but decreased in D. aldrichi. Wing development in the two species thus is very different, with the proximal part of the wing in D. buzzatii more closely allied to the thorax than to the distal part. Among line variation was significant for all traits in both species, and most pronounced for thorax length and the ratio parameters. Coefficients of variation were significantly different between the species for all traits, with those in D. aldrichi higher than in D. buzzatii. Genetic variance in plasticity was significant for all traits in D. buzzatii, but only for seven out of 12 in D. aldrichi. Additive genetic variances for all traits in both species were significantly larger than zero. Genetic correlations between thorax length and several wing length parameters, and between these and wing area, were positive and generally significant in both species. The genetic correlation between the distal and the proximal length of the third vein was not significantly different from zero in D. aldrichi, but negative and significant in D. buzzatii. Heritabilites varied significantly among temperatures for almost all traits in both species. Phenotypic variances were generally higher in D. aldrichi than in D. buzzatii, and commonly highest at the extreme temperatures in the former species. At the high temperature the genetic variances also were usually highest in D. aldrichi. The data clearly suggest that the process of thermal adaptation is species specific and caution against generalizations based on the study of single species.     

A study of wing morphology and fluctuating asymmetry in interspecific hybrids between Drosophila buzzatii and D. koepferae

In this work we investigate the effect of interspecific hybridization on wing morphology using geometric morphometrics in the cactophilic sibling species D. buzzatii and D. koepferae. Wing morphology in F₁ hybrids exhibited an important degree of phenotypic plasticity and differs significantly from both parental species. However, the pattern of morphological variation between hybrids and the parental strains varied between wing size and wing shape, across rearing media, sexes, and crosses, suggesting a complex genetic architecture underlying divergence in wing morphology. Even though there was significant fluctuating asymmetry for both, wing size and shape in F₁ hybrids and both parental species, there was no evidence of an increased degree of fluctuating asymmetry in hybrids as compared to parental species. These results are interpreted in terms of developmental stability as a function of a balance between levels of heterozygosity and the disruption of coadaptation as an indirect consequence of genomic divergence.     

Variation of life-history and morphometrical traits in Drosophila buzzatii and Drosophila simulans collected along an altitudinal gradient from a Canary island

Variation in three life‐history traits (developmental time, preadult viability and daily female productivity) and five morphometrical traits (thorax length, wing length, wing width, wing/thorax ratio and wing‐aspect ratio) was studied at three developmental temperatures (20, 25 and 30 °C) in Drosophila buzzatii and Drosophila simulans collected on the island of La Gomera (Canary Archipelago). The flies originated from five closely situated localities, representing different altitudes (from 20 to 886 m above sea level) and a range of climatic conditions. We found statistically significant population effects for all traits in D. buzzatii and for most of the traits in D. simulans. Although no correlations of trait values with altitude were detected, geographical patterns for three life‐history traits and body size in D. buzzatii indicated that short‐range geographical variation in this species could be maintained by local climatic selection. Five of eight traits showed population‐by‐temperature interactions either in D. buzzatii or in D. simulans, but in all cases except wing width in D. buzzatii this could not be interpreted as adaptive responses to thermal conditions in the localities. The range of plastic changes across temperatures for particular traits differed between species, indicating a possibility for different levels of environmental stress experienced by the natural populations. The reaction norm curves and the response of within‐population variability to thermal treatments suggested better adaptations to higher and lower temperatures for D. buzzatii and D. simulans, respectively. The levels of among‐population differentiation depended on developmental temperature, implying environmental effects on the expression of the genetic variance. At 20 and 25 °C, interpopulation variability in D. buzzatii was higher than in D. simulans, while at 30 °C the opposite trend was observed. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 84 , 119–136.     

Differential Rates of Male Genital Evolution in Sibling Species of <Emphasis Type="Italic">Drosophila</Emphasis>

Genital morphology in animals with internal fertilization is considered to be among the fastest evolving traits. Sexual selection is often proposed as the main driver of genital diversification but the exact selection mechanisms involved are usually unclear. In addition, the mechanisms operating may differ even between pairs of sibling species. We investigated patterns of male genital variation within and between natural populations of the cactophilic fly Drosophila koepferae ranging its entire geographic distribution and compared them with those previously observed in its sibling species, D. buzzatii. Using both mtDNA and nDNA markers we found that genital shape variation in D. koepferae is more restricted than expected for neutral evolution, suggesting the predominance of stabilizing selection. We also detected dissimilar patterns of divergence between populations of D. koepferae that were allopatric and sympatric with D. buzzatii. The constrained evolution inferred for D. koepferae’s genitalia clearly contrasts with the rapid divergence and higher morphological disparity observed in the populations of D. buzzatii. Finally, different possible scenarios of male genital evolution in each species and within the radiation of D. buzzatii cluster are discussed.     

Genetic variation and plasticity of thorax length and wing length in Drosophila aldrichi and D. buzzatii

Reaction norms across three temperatures of development were measured for thorax length, wing length and wing length/thorax length ratio for ten isofemale lines from each of two populations of Drosophila aldrichi and D. buzzatii. Means for thorax and wing length in both species were larger at 24 °C than at either 18 °C or 31 °C, with the reduction in size at 18 °C most likely due to a nutritional constraint. Although females were larger than males, the sexes were not different for wing length/thorax length ratio. The plasticity of the traits differed between species and between populations of each species, with genetic variation in plasticity similar for the two species from one locality, but much higher for D. aldrichi from the other. Estimates of heritabilities for D. aldrichi generally were higher at 18 °C and 24 °C than at 31 °C, but for D. buzzatii they were highest at 31 °C, although heritabilities were not significantly different between species at any temperature. Additive genetic variances for D. aldrichi showed trends similar to that for heritability, being highest at 18 °C and decreasing as temperature increased. For D. buzzatii, however, additive genetic variances were lowest at 24 °C. These results are suggestive that genetic variation for body size characters is increased in more stressful environments. Thorax and wing lengths showed significant genetic correlations that were not different between the species, but the genetic correlations between each of these traits and their ratio were significantly different. For D. aldrichi, genetic variation in the wing length/thorax length ratio was due primarily to variation in thorax length, while for D. buzzatii, it was due primarily to variation in wing length. The wing length/thorax length ratio, which is the inverse of wing loading, decreased linearly as temperature increased, and it is suggested that this ratio may be of greater adaptive significance than either of its components.     

Genetic and phenotypic variation for flight ability in the cactophilicDrosophila species,D. aldrichi andD. buzzatii

The flight ability ofDrosophila aldrichi (Patterson & Crow) andD. buzzatii (Patterson & Wheeler) using tethered flights, was measured with respect to age-related changes, genetic variation and adult body size variation induced by rearing at different larval densities.Drosophila buzzatii flew for much longer thanD. aldrichi, especially females, but age-related changes in flight duration were significant only forD. aldrichi. Effects of body size on flight ability were significant inD. buzzatii, but not inD. aldrichi. InD. buzzatii, there was a significant genotype-environment interaction (larval density × line) for flight duration, with short and average flight duration isofemale lines showing longer flights, but a long flight duration line shorter flights as body size decreased (i.e., as larval density increased). Heritability estimates for flight duration were similar in the two species, but flight duration showed no significant genetic correlations with developmental time, body size or wing dimensions (except for one wing dimension inD. buzzatii). Although not significantly different between the species, heritabilities for life-history traits (adult size and developmental time) showed contrasting patterns — with higher heritability for body size (body weight and thorax length) inD. buzzatii, and higher for developmental time inD. aldrichi. In agreement with limited previous field evidence,D. buzzatii is better adapted for colonization than isD. aldrichi.     

A three-locus system of interspecific incompatibility underlies male inviability in hybrids between Drosophila buzzatii and D. koepferae

In hybrids between the sibling species D. buzzatii and D. koepferae, both sexes are more or less equally viable in the F₁: However, backcross males to D. buzzatii are frequently inviable, apparently because of interspecific genetic incompatibilities that are cryptic in the F₁. We have performed a genetic dissection of the effects of the X chromosome from D. koepferae. We found only two cytological regions, termed hmi-1 and hmi-2, altogether representing 9% of the whole chromosome, which when introgressed into D. buzzatii cause inviability of hybrid males. Observation of the pattern of asynapsis of polytene chromosomes (incomplete pairing, marking introgressed material) in females and segregation analyses were the technique used to infer the X chromosome regions responsible for this hybrid male inviability. The comparison of these results with those previously obtained with the same technique for hybrid male sterility in this same species pair indicate that in the X chromosome of D. koepferae there are at least seven times more regions that produce hybrid male sterility than hybrid male inviability. We have also found that the inviability brought about by the introgression of hmi-1 is suppressed by the cointrogression of two autosomal sections from D. koepferae. Apparently, these three regions conform to a system of species-specific complementary factors involved in an X-autosome interaction that, when disrupted in backcross hybrids by recombination with the genome of its sibling D. buzzatii, brings about hybrid male inviability.     

Multiple paternity and sperm competition in the sibling species Drosophila buzzatii and Drosophila koepferae

Sperm competition (SC) is a major component of sexual selection that enhances intra‐ and intersexual conflicts and may trigger rapid adaptive evolution of sexual characters. The actual role of SC on rapid evolution, however, is poorly understood. Besides, the relative contribution of distinctive features of the mating system to among species variation in the strength of SC remains unclear. Here, we assessed the strength of SC and mating system factors that may account for it in the closely related species Drosophila buzzatii and Drosophila koepferae. Our analyses reveal higher incidence of multiple paternity and SC risk in D. buzzatii wild‐inseminated females. The estimated number of fathers per brood was 3.57 in D. buzzatii and 1.95 in D. koepferae. In turn, the expected proportion of females inseminated by more than one male was 0.89 in D. buzzatii and 0.58 in D. koepferae. Laboratory experiments show that this pattern may be accounted for by the faster rate of stored sperm usage observed in D. koepferae and by the greater female remating rate exhibited by D. buzzatii. We also found that the male reproductive cost of SC is also higher in D. buzzatii. After a female mated with a second male, first‐mating male fertility was reduced by 71.4% in D. buzzatii and only 33.3% in D. koepferae. Therefore, we may conclude that postmating sexual selection via SC is a stronger evolutionary force in D. buzzatii than in its sibling.     

Molecular Population Genetics of the α<Emphasis Type="Italic">-Esterase5</Emphasis> Gene Locus in Original and Colonized Populations of <Emphasis Type="Italic">Drosophila buzzatii</Emphasis> and Its Sibling <Emphasis Type="Italic">Drosophila koepferae</Emphasis>

Several studies have suggested that esterase-2 (EST-2) may be the target of natural selection in the cactophilic fly Drosophila buzzatii. In this work, we analyzed nucleotide variation in a fragment of α-esterase5 (αE5), the gene encoding EST-2, in original (Argentinian) and colonized (Australian) populations of D. buzzatii and in its sibling D. koepferae. Estimates of nucleotide heterozygosity in D. buzzatii were similar in Australia and Argentina, although we detected a loss of singletons in colonized populations, suggesting a moderate founder effect. Interspecific comparisons revealed that D. buzzatii was more polymorphic for nonsynonymous variation, whereas D. koepferae was more variable for synonymous and noncoding sites. The two major chromosomal arrangements (2st and 2j) in D. buzzatii displayed similar levels of nucleotide variation, whereas 2jz ³ was monomorphic. The sequenced region allowed the discrimination of a greater number of EST-2 protein variants in the Australian sample than in the Argentinean sample. In D. koepferae, nucleotide variation in αE5 does not depart from neutral expectations, although tests of population structure were significant for silent variation. In contrast, D. buzzatii has probably undergone a recent population expansion in its South American range. In addition, the McDonald and Kreitman test revealed an excess of nonsynonymous polymorphism in both original and colonized populations of this species. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Richard Kliman]     

Chromosomal Evolution of Sibling Species of the Drosophila willistoni Group. I. Chromosomal Arm IIR (Muller’s Element B)

The phylogenetic relationships among nine entities of Drosophila belonging to the D. willistoni subgroup were investigated by establishing the homologous chromosomal segments of IIR chromosome, Muller’s element B (equivalent to chromosome 2L of D. melanogaster). The sibling species of the D. willistoni group investigated include D. willistoni, D. tropicalis tropicalis, D. tropicalis cubana, D. equinoxialis, D. insularis and four semispecies of the D. paulistorum complex. The phylogenetic relationships were based on the existence of segments in different triads of species, which could only be produced by overlapping inversions. Polytene banding similarity maps and break points of inversions between species are presented. The implications of the chromosomal data for the phylogeny of the species and comparisons with molecular data are discussed. The aim of this study is to produce phylogenetic trees depicting accurately the sequence of natural events that have occurred in the evolution of these sibling species. Claudia Rohde, Ana Cristina Lauer Garcia: These authors contributed equally to this work     

High transposition rates of Osvaldo, a new Drosophila buzzatii retrotransposon

Transposition of a new Drosophila retrotransposon was investigated. Total genomic Southern analysis and polytene in situ hybridizations in D. buzzatii strains and other related species using a 6 kb D. buzzatii clone (cDb314) showed a dispersed, repetitive DNA pattern, suggesting that this clone contains a transposable element (TE). We have sequenced the cDb314 clone and demonstrated that it contains all the conserved protein sequences and motifs typical of retrovirus-related sequences. Although cDb314 does not include the complete TE, the protein sequence alignment demonstrates that it includes a defective copy of a new long terminal repeat (LTR) retrotransposon, related to the gypsy family, which we have named Osvaldo. Using a D. buzzatii inbred line in which all insertion sites are known, we have measured Osvaldo transposition rates in hybrids between this D. buzzatii line and its sibling species D. koepferae. The results show that Osvaldo transposes in bursts at high rate, both in the D. buzzatii inbred line and in species hybrids.This paper is dedicated posthumously to Osvaldo A. Reig in recognition of his contributions to evolutionary biology and his early appreciation of the role of transposable elements in evolution     

<Emphasis Type="Italic">Drosophila</Emphasis> pupation behavior in the wild

We investigated pupa distributions of D. simulans, D. buzzatii, D. melanogaster, D. immigrans and D. hydei on a number of natural breeding sites. Pupae of all five species showed aggregated distributions, which prompted us to examine these aggregations in a more detail for two species that commonly co-occur in breeding sites, D. simulans and D. buzzatii. We found that pupae of both species tend to be aggregated in conspecific clusters. Subsequent experiments revealed that both species are attracted to the odors of other larvae, though only D. buzzatii differentiated between conspecifics and heterospecifics (they preferred conspecific). Furthermore, third instar larvae of both species preferred more alkaline substrates. Altogether, our results demonstrate that Drosophila species form conspecific pupa aggregations in natural breeding sites, and that pupation site selection depends on interactions among conspecific and heterospecific larvae and on chemical characteristics of the breeding sites.     

Aedeagal Divergence in Sympatric Populations of Two Sibling Species of Cactophilic Drosophila (Diptera: Drosophilidae): Evidence of Character Displacement?

Aedeagal morphology of two sibling cactophilic species, Drosophila buzzatii Patterson & Wheeler and Drosophila koepferae Fontdevila & Wasserman, was analyzed in nine allopatric and three sympatric locations throughout South America. Morphological differences were detected for both aedeagus size and shape between sympatric and allopatric populations of D. buzzatii, despite the significant variability within both groups. Populations of D. buzzatii sympatric with D. koepferae displayed smaller aedeagus than the allopatric ones as well as more differentiated aedeagus shape. The shape differences were non-allometric and mainly consisted in a change of curvature of the dorsal margin of the aedeagus being more pronounced in males from populations sympatric with D. koepferae. It is concluded that aedeagal morphology presented some degree of character displacement in both size and shape in populations of D. buzzatii in sympatry with D. koepferae. These results might suggest the existence of mechanisms of interspecific recognition and hybridization prevention between these species that include the morphology of the male genitalia.     

Host alkaloids differentially affect developmental stability and wing vein canalization in cactophilic Drosophila buzzatii

Host shifts cause drastic consequences on fitness in cactophilic species of Drosophila. It has been argued that changes in the nutritional values accompanying host shifts may elicit these fitness responses, but they may also reflect the presence of potentially toxic secondary compounds that affect resource quality. Recent studies reported that alkaloids extracted from the columnar cactus Trichocereus terscheckii are toxic for the developing larvae of Drosophila buzzatii. In this study, we tested the effect of artificial diets including increasing doses of host alkaloids on developmental stability and wing morphology in D. buzzatii. We found that alkaloids disrupt normal wing venation patterning and affect viability, wing size and fluctuating asymmetry, suggesting the involvement of stress–response mechanisms. Theoretical implications are discussed in the context of developmental stability, stress, fitness and their relationship with robustness, canalization and phenotypic plasticity.     

Host‐dependent phenotypic plasticity of aedeagus morphology in a pair of cactophilic sibling Drosophila species of the repleta group (Diptera,Drosophilidae)

The rapid evolution of male genital morphology is a characteristic feature of several animal groups. Such rapid divergence makes this trait a useful key for species identification. The aedeagus, the intromittent organ of male genitalia, is considered the main diagnostic trait in the Drosophila repleta group. In this study we analysed phenotypic plasticity and genetic variations associated with aedeagus size and shape in the cactophilic sibling species Drosophila gouveai Tidon‐Sklorz and Sene, 2001 and Drosophila antonietae Tidon‐Sklorz and Sene, 2001. Phenotypic plasticity in aedeagus morphology was evaluated in terms of the response to rearing media prepared with each species’ natural host plant, Pilosocereus machrisii Dawson, 1957 and Cereus hildmannianus Schum, 1890 respectively. Our results show that aedeagal shape differed significantly between species and that both shape and size presented host‐related phenotypic plasticity in both species. Flies reared on P. machrisii had, on average, larger aedeagi than those grown in C. hildmannianus. The general shape of aedeagus also differed significantly between flies that emerged in different host cactus. Patterns of variation in aedeagus morphology are discussed in the light of the current knowledge of evolutionary relationships and host plant use, in the D. buzzatii cluster, an assemblage of species in active cladogenesis.     

Evolution of male genitalia: environmental and genetic factors affect genital morphology in two <Emphasis Type="Italic">Drosophila</Emphasis> sibling species and their hybrids

Background  

The rapid evolution of genital morphology is a fascinating feature that accompanies many speciation events. However, the underlying patterns and explanatory processes remain to be settled. In this work we investigate the patterns of intraspecific variation and interspecific divergence in male genitalic morphology (size and shape) in the cactophilic sibling species Drosophila buzzatii and D. koepferae. Genital morphology in interspecific hybrids was examined and compared to the corresponding parental lines.     

Cactophilic Drosophila in South America: A Model for Evolutionary Studies

The Drosophila buzzatii cluster is composed of seven cactophilic species and their known geographical distribution encompasses the open vegetation diagonal, which includes the morphoclimatic Domains of the Caatinga, Chaco and Cerrado, which are situated between the Amazon and the Atlantic forests. Besides these areas, these cactophilic species are also found in a narrow strip along the Atlantic coast from northeastern Brazil to the southern tip of the country. The hypothesis of vicariant events, defining the core areas of each species, is proposed to explain the historical diversification for the cluster. The intraspecific analysis for the cluster shows a population structure with gene flow restricted by distance, range expansion with secondary contact resulting in introgression and simpatry, especially in the limits of the species distribution, polytypic populations and assortative mating in inter population experiments. There is a variation related to these events that depends on the species and geographic origin of the population analyzed. These events are, hypothetically, described as the results of expansion and retraction of the population ranges, as a consequence of their association with cacti, which theoretically follow the expansion and retraction of dry areas during the paleoclimatic oscillations in South America, as that promoted by the glacial cycles of the Quaternary. The Drosophila buzzatii cluster is divided into two groups. The first one is composed of D. buzzatii, a species that has a broad geographic distribution and no significant differentiation between its populations. The second is the Drosophila serido sibling set, which encompasses the others species and is characterized by a significant potential for differentiation.