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.     

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.     

Genotype by environment interactions in viability and developmental time in populations of cactophilic Drosophila

The genetic and ecological basis of viability and developmental time differences between Drosophila buzzatii and D. koepferae were analysed using the isofemale line technique. Several isofemale lines were sampled from pairs of allopatric/sympatric populations of each species. Flies were reared in media prepared with decaying tissues of two of the main natural cactus hosts of each species. This experimental design enabled us to evaluate the relative contribution of phenotypic plasticity, genetic variation and genotype by environment interaction (G x E) to total phenotypic variation for two fitness traits, viability and developmental time. Our results revealed significant G x E in both traits, suggesting that the maintenance of genetic variation can be explained, at least in part, by diversifying selection in different patches of a heterogeneous environment in both species. However, the relative importance of the factors involved in the G x E varied between traits and populations within species. For viability, the G x E can be mainly attributed to changes in the rank order of lines across cacti. However, the pattern was different for developmental time. In D. buzzatii the G x E can be mainly accounted for by changes in among line variance across cacti, whereas changes in the rank order of lines across cacti was the main component in D. koepferae. These dissimilar patterns of variation between traits and species suggest that the evolutionary forces shaping genetic variation for developmental time and viability vary between populations within species and between species.     

Nucleotide polymorphism at the xanthine dehydrogenase locus in Drosophila pseudoobscura.

Sequential polyacrylamide electrophoresis has revealed 20 allozymes of xanthine dehydrogenase (XDH) in Drosophila pseudoobscura. DNA sequence determination of seven isolates of the Xdh locus that represent six allozyme classes are presented here. Of the 5,456 sites examined, 180 are polymorphic, with 27 polymorphisms occurring at nonsynonymous, or replacement, sites. An average of nine amino acids differ between XDH allozyme classes, with 85% of the polymorphic amino acids singly represented. The level and pattern of variation observed at Xdh argue that the effective population size of the species is quite large--i.e., on the order of 2 x 10(6)--and that the populations sampled are quite ancient. In addition, as judged by two statistical tests, the levels of nucleotide polymorphism observed at Xdh are compatible with predictions from the neutral theory of molecular evolution.     

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.     

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.     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

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.     

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.     

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.     

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 variation in original and colonizing Drosophila buzzatii populations analysed by microsatellite loci isolated with a new PCR screening method

A new polymerase chain reaction-based screening method for microsatellites is presented. Using this method, we isolated 12 microsatellite loci from Drosophila buzzatii, two of which were X-linked. We applied the other 10 microsatellite loci to the analysis of genetic variation in five natural populations of D. buzzatii. Two populations were from the species' original distribution in Argentina, whereas the other three were from Europe (two) and Australia that were colonized 200 and 65 years ago, respectively. Allelic variation was much larger in the original populations than in the colonizing ones and there was a tendency to decreased heterozygosity in the colonizing populations. We used three different statistical procedures for detecting population bottlenecks. All procedures suggested that the low variability in the populations in the Old World was not the result of the recent population decline, but was due to a founder effect followed by a population expansion. In fact, one procedure which detects population expansions and declines based on the genealogical history of microsatellite data suggested that an expansion had taken place in all the colonized populations.     

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.     

Synonymous Substitutions in the Xdh Gene of Drosophila: Heterogeneous Distribution along the Coding Region

The Xdh (rosy) region of Drosophila subobscura has been sequenced and compared to the homologous region of D. pseudoobscura and D. melanogaster. Estimates of the numbers of synonymous substitutions per site (Ks) confirm that Xdh has a high synonymous substitution rate. The distributions of both nonsynonymous and synonymous substitutions along the coding region were found to be heterogeneous. Also, no relationship has been detected between Ks estimates and codon usage bias along the gene, in contrast with the generally observed relationship among genes. This heterogeneous distribution of synonymous substitutions along the Xdh gene, which is expression-level independent, could be explained by a differential selection pressure on synonymous sites along the coding region acting on mRNA secondary structure. The synonymous rate in the Xdh coding region is lower in the D. subobscura than in the D. pseudoobscura lineage, whereas the reverse is true for the Adh gene.     

Postglacial population growth of Cunninghamia konishii (Cupressaceae) inferred from phylogeographical and mismatch analysis of chloroplast DNA variation

Phylogeographical and mismatch analysis of chloroplast DNA (cpDNA) variation were used to infer the temporal dynamics of distributional and demographic history of Taiwan fir (Cunninghamia konishii). We examined 64 and 52 trees from 17 populations of C. konishii and 14 provenances of C. lanceolata, respectively, by sequencing three intergenic spacers and one intron using cpDNA universal primers. Of the aligned 1888 base pairs (bp) sequence, 30 varied among 28 haplotypes, which consisted of three transitions, 14 transversions and 13 indels. One ancestral haplotype was found in 86 individuals across the surveyed range of both species, C. konishii and C. lanceolata, which was distributed in all populations and provenances. The 28 haplotypes also included 15 C. konishii specific and 12 C. lanceolata-specific haplotypes. Ancestral haplotype was found fixed in five populations of C. konishii and five provenances of C. lanceolata. Other haplotypes occurred mainly as singletons. The levels of population differentiation studied are relatively low in both Cunninghamia species. The nucleotide diversity (theta) of chloroplast DNA sequences within C. konishii was slightly higher than that of C. lanceolata. Excess in singletons as well as star-like phylogeny of haplotypes suggested no clearcut migration patterns of C. konishii after glacial maximum. One probable demographic history of C. konishii is the postglacial population growth of C. konishii after a glacial bottleneck event. This inference is supported by the combined results of fossil pollen record, low nucleotide diversity, significant Tajima's d-value, phylogeographical analysis and unimodal mismatch distribution. Similarities and discrepancies between our results and those of Lu et al. (2001) are discussed.     

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]     

基于线粒体COI基因分析钩手水母的群体遗传结构

钩手水母(Gonionemus vertens)为大西洋和太平洋广布种, 是我国习见的有毒水母种类之一。本文对采自黄渤海海域4个地理群体的104个钩手水母线粒体COI基因序列进行扩增, 并结合GenBank上其他182个钩手水母同源序列进行序列变异分析。在286个基因序列中共检测出52个多态位点, 定义了14种单倍型。总群体的单倍型多样性和核苷酸多样性分别为0.743 ± 0.012和1.046% ± 0.097%, 与其他几种大型水母相比, 钩手水母总群体的遗传多样性处于较高水平。AMOVA结果显示, 60.17%的分子变异源于群组间, 13.37%的分子变异源于群体内, 26.46%的分子变异源于组内群体间, 群组间、群体内和组内群体间的遗传分化均极显著。F_st值统计检验表明, 中国厦门群体与乐亭、东营、烟台、大连群体间存在显著的遗传分化, 大连与东营、烟台群体间也存在显著的遗传分化。系统分析结果显示, 钩手水母群体间存在2个明显的单倍型谱系分支。不同的钩手水母地理群体间具有复杂的遗传模式, 钩手水母复杂的生活史、扩散能力、地理隔离和海流分布可能是影响钩手水母遗传结构的重要因素。     

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.