Transcriptional responses are oriented towards different components of the rearing environment in two Drosophila sibling species.

BACKGROUND: The chance to compare patterns of differential gene expression in related ecologically distinct species can be particularly fruitful to investigate the genetics of adaptation and phenotypic plasticity. In this regard, a powerful technique such as RNA-Seq applied to ecologically amenable taxa allows to address issues that are not possible in classic model species. Here, we study gene expression profiles and larval performance of the cactophilic siblings Drosophila buzzatii and D. koepferae reared in media that approximate natural conditions and evaluate both chemical and nutritional components of the diet. These closely related species are complementary in terms of host-plant use since the primary host of one is the secondary of the other. D. koepferae is mainly a columnar cactus dweller while D. buzzatii prefers Opuntia hosts. RESULTS: Our comparative study shows that D. buzzatii and D. koepferae have different transcriptional strategies to face the challenges posed by their natural resources. The former has greater transcriptional plasticity, and its response is mainly modulated by alkaloids of its secondary host, while the latter has a more canalized genetic response, and its transcriptional plasticity is associated with the cactus species. CONCLUSIONS: Our study unveils a complex pleiotropic genetic landscape in both species, with functional links that relate detox responses and redox mechanisms with developmental and neurobiological processes. These results contribute to deepen our understanding of the role of host plant shifts and natural stress driving ecological specialization.

Studies on the Biology of Hypogeococcus pungens (sensu stricto) (Hemiptera: Pseudococcidae) in Argentina to Aid the Identification of the Mealybug Pest of Cactaceae in Puerto Rico.

Hypogeococcus pungens Granara de Willink, sensu stricto, is a serious pest of cacti in Puerto Rico threating many Caribbean islands. A classical biological control program for H. pungens was initiated for Puerto Rico in 2010 with a survey for natural enemies of H. pungens in its native range of Argentina. Biological differences were observed between populations of H. pungens sampled on Amaranthaceae and Cactaceae. Molecular studies suggested that H. pungens populations from different host plant families are likely a complex of species. Our objective was to study the biology of H. pungens sensu stricto on specimens collected in the same locality and host plant as the holotype [Tucumán Province, Argentina; Alternanthera pungens Kunth (Amaranthaceae)]. We were interested in the reproductive biology of females, longevity and survival of adults, the effect of temperature on the development, and nymph performance (survival and development) on five Cactaceae species. We found that H. pungens s.s showed marked biological differences from the populations collected on Cactaceae and exported to Australia for the biological control of the cactus Harrisia spp. The main differences were the presence of deuterotoky parthenogenesis and the fact that H. pungens did not attack Cactaceae in the laboratory. Our results provide biological evidence that H. pungens is a species complex. We propose that the population introduced to Australia is neither Hypogeococcus festerianus Lizer y Trelles nor H. pungens, but an undescribed species with three circuli, and that the Hypogeococcus pest of cacti in Puerto Rico is not H. pungens.

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.

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.

Tracking the origin of an invasive species: Drosophila subobscura in Argentina.

Biological invasions are excellent opportunities to study the evolutionary forces leading to the adaptation of a species to a new habitat. Knowledge of the introduction history of colonizing species helps tracking colonizing routes and assists in defining management strategies for invasive species. The Palearctic species Drosophila subobscura is a good model organism for tracking colonizations since it was detected in Chile and western North America three decades ago and later on in the Atlantic coast of Argentina. To unravel the origin of the Argentinean colonizers two populations have been analysed with several genetic markers. Chromosomal arrangements and microsatellite alleles found in Argentina are almost similar to those observed in Chile and USA. The lethal allelism test demonstrates that the lethal gene associated with the O(5) inversions in Argentina is identical to that found in Chile and USA, strongly supporting the hypothesis that all the American colonizing populations originated from the same colonization event. A secondary bottleneck is detected in the Argentinean populations and the genetic markers suggest that these populations originated from the invasion of 80-150 founding individuals from Chile.

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.

First Record of Drosophila buzzatii (Patterson & Wheeler) (Diptera: Drosophilidae) Emerging from a Non-Cactus Host.

Drosophila buzzatii (Patterson & Wheeler), a typical cactophilic species of the repleta group, is registered for the first time emerging from Melon (Cucumis melo) in western Argentina. The analysis of inversion polymorphism and genetic diversity of mitochondrial cytochrome oxidase subunit I gene (mtCOI) provided additional evidence that corroborated the presence of a high proportion of D. buzzatii among the flies emerged from melon. This finding set the scenario for a broader range of possible hosts and host-related distribution and dispersion for this widespread species.

Contrasting histories of three gene regions associated with In(3L)Payne of Drosophila melanogaster.

In the present report, we studied nucleotide variation in three gene regions of Drosophila melanogaster, spanning > 5 kb and showing different degrees of association with the cosmopolitan inversion In(3L)Payne. The analysis of sequence variation in the regions surrounding the breakpoints and the heat shock 83 (Hsp83) gene locus, located close to the distal breakpoint, revealed the absence of shared polymorphisms and the presence of a number of fixed differences between arrangements, indicating absence of genetic exchange. In contrast, for the esterase-6 gene region, located in the center of the inversion, we observed the presence of shared polymorphisms between arrangements suggesting genetic exchange. In the regions close to the breakpoints, the common St arrangement is 10 times more polymorphic than inverted chromosomes. We propose that the lack of recombination between arrangements in these regions coupled with genetic hitchhiking is the best explanation for the low heterozygosity observed in inverted lines. Using the data for the breakpoints, we estimate that this inversion polymorphism is around 0.36 million yr old. Although it is widely accepted that inversions are examples of balanced polymorphisms, none of the current neutrality tests including our Monte Carlo simulations showed significant departure from neutral expectations.

Genetic variance for body size in a natural population of Drosophila buzzatii.

Previous work has shown thorax length to be under directional selection in the Drosophila buzzatii population of Carboneras. In order to predict the genetic consequences of natural selection, genetic variation for this trait was investigated in two ways. First, narrow sense heritability was estimated in the laboratory F2 generation of a sample of wild flies by means of the offspring-parent regression. A relatively high value, 0.59, was obtained. Because the phenotypic variance of wild flies was 7-9 times that of the flies raised in the laboratory, "natural" heritability may be estimated as one-seventh to one-ninth that value. Second, the contribution of the second and fourth chromosomes, which are polymorphic for paracentric inversions, to the genetic variance of thorax length was estimated in the field and in the laboratory. This was done with the assistance of a simple genetic model which shows that the variance among chromosome arrangements and the variance among karyotypes provide minimum estimates of the chromosome's contribution to the additive and genetic variances of the trait, respectively. In males raised under optimal conditions in the laboratory, the variance among second-chromosome karyotypes accounted for 11.43% of the total phenotypic variance and most of this variance was additive; by contrast, the contribution of the fourth chromosome was nonsignificant. The variance among second-chromosome karyotypes accounted for 1.56-1.78% of the total phenotypic variance in wild males and was nonsignificant in wild females. The variance among fourth chromosome karyotypes accounted for 0.14-3.48% of the total phenotypic variance in wild flies. At both chromosomes, the proportion of additive variance was higher in mating flies than in nonmating flies.

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

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 2jq7 had longer DTs and larger TLs. Arrangements 2st and 2jq7 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.

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.

Behavioral differentiation in oviposition activity in Drosophila buzzatii from highland and lowland populations in Argentina: plasticity or thermal adaptation?

Highland populations of several Drosophila species in Argentina were active early in the afternoon in the field as opposed to populations from a much warmer lowland site, where flies were mainly active in the early evening prior to sunset. For one of these species, Drosophila buzzatii, we tested for a genetic component of activity differences by carrying out crosses within and between populations and measuring oviposition activity of the progeny in the laboratory. We found that activity in the highland population exceeded that in the lowland one during the midafternoon, whereas activity in the lowland population exceeded that in the highland one prior to the beginning of the dark period. Oviposition activity for the period corresponding to the field observations was regressed on the proportion of the genome derived from the highland population. This variable significantly predicted oviposition activity between 1400 and 1600 and between 2000 and 2200 h. Activity of both reciprocal crosses was intermediate and not significantly different from each other, suggesting that nuclear genetic, rather than cytoplasmic factors contribute to differences in oviposition activity between the populations. Two morphological, one genetic, and one stress resistance trait were also scored to examine whether temperature differences between environments were associated with other differences between populations. Wing length of wild-caught and laboratory-reared flies from the highland population significantly exceeded that in the lowland. Thorax length of laboratory-reared flies from the highland population also significantly exceeded that from the lowland. Chromosomal inversion frequencies differed significantly between the two populations with a fivefold reduction in the frequency of arrangement 2st in the highland as compared to the lowland population. This arrangement is known for its negative dose effect on size, and thus, the highland population has experienced a genetic change, perhaps as a result of adaptation to the colder environment, where body size and the frequency of arrangement 2st have changed in concert. Finally, a heat knockdown test revealed that the lowland population was significantly more resistant to high temperature than the highland one. In conclusion, we suggest that temperature has been an important selective agent causing adaptive differentiation between these two populations. We also suggest that the activity rhythms of the two populations have diverged as a consequence of behavioral evolution, that is, through avoidance of stressful temperatures as a mean of thermal adaptation.

Does salinity induce early ageing of pea root tissue?

The effect of salinity on ageing of pea roots was studied. The distance from the apex at which differentiation of xylem elements occurred and the relative increase in the function of pentose phosphate pathways were taken as parameters for maturation or ageing.Pea seeds (Pisum stivum L.) of the varieties Alaska and Dan were used in these experiments. The seeds were germinated and grown in vermiculite moistened with Hoagland's solution or Hoagland's solution containing either 96 or 120 mM NaCl. In Alaska roots salinity induced differentiation in a lower section of the root than in controls, and the increase in the function of the pentose phosphate pathway paralleled the advance of maturation. Salinity apparently induces earlier ageing in Alaska roots. This is not the case in Dan roots which tolerate slightly higher salinity levels than Alaska.

Spontaneous second-chromosome recessive lethals in a Drosophila melanogaster repair-deficient (mei-9a) mutant.

The frequency of spontaneous second-chromosome recessive lethal mutations was determined in a Drosophila melanogaster strain kept homozygous for the X-linked mei-9 mutant for 50 generations. The data show that the rate of lethals of independent origin is 9 times higher and that of semilethals 4 times higher than the control values. The mei-9 mutant, originally isolated on the ground of its abnormal behaviour during meiosis, was later found to be defective in excision repair in somatic cells and to yield, after exposure to different physical and chemical agents, higher frequencies of genetic damage than mei-9+ flies. Our results are taken as evidence that the normal functioning of the mei-9 locus plays an important role in the repair of spontaneous lesions in the second chromosome and give support to positive findings in the X chromosome. Monitoring the accumulation of mutations in an experimental population could be a useful methodological approach to detect slight, though relevant, genetic damage induced by long-term exposure to environmental mutagens, whose detection is beyond the resolving power of conventional tests in Drosophila.

Male genital and wing morphology in the cactophilic sibling species Drosophila gouveai and Drosophila antonietae and their hybrids reared in different host plants.

Cactophilic Drosophila flies are excellent models to study adaptation to a relatively narrow spectrum of potential host plants and host-driven evolutionary diversification. Previous studies suggested a complex genetic architecture of wing and male genital morphology in phylogenetically basal species of the D. buzzatii cluster. In this work, we investigate the effect of experimental hybridization and host plant shifts on male genital and wing morphology in D. gouveai Tidon-Sklorz and Sene and D. antonietae Tidon-Sklorz and Sene, a pair of more recently derived species. We explicitly tested the hypotheses that wing and male genital morphology in interspecific hybrids depend on the host plant in which flies were grown. Our study shows that cactus hosts exert a strong effect on genital and wing morphology and that hybrids can be clearly differentiated on the basis of wing and genital morphology from both parental species. However, the extent of morphological differentiation between hybrids and pure species as well as plasticity patterns varied across organs, suggesting a complex genetic architecture for the studied traits.

Molecular population genetics of the alpha-esterase5 gene locus in original and colonized populations of Drosophila buzzatii and its sibling Drosophila koepferae.

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 alpha-esterase5 (alphaE5), 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 2jz3 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 alphaE5 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.

Mitochondrial sequence diversity of the southernmost extant New World monkey, Alouatta caraya.

Variability in mitochondrial DNA sequences was analyzed in the howler monkey, Alouatta caraya, in order to delineate evolutionary relationships among populations in the most southerly distributed New World monkey. Based on new and previously published sequence data, fourteen cytochrome b haplotypes were observed among 33 howlers sampled in Argentina, Paraguay and Brazil, and grouped in two main haplogroups. In northeastern Argentina and southern Paraguay, new sequence data on 73 specimens sampled from six localities gave 34 control region haplotypes that also clustered in two main haplogroups. At this southern distribution, both mitochondrial markers revealed the presence of two sympatric and differentiated clades that we interpret to be the consequence of a secondary contact between previously allopatric populations. Given evidence for a demographic expansion at the beginning of the Holocene 15,500-7000 years ago (Fu's test, F(S)=-12.137; P<0.001), we suggest that atleast two populations of A. caraya have colonized the southernmost range since the Holocene employing forested corridors on the Paraná and Paraguay Rivers.

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.