Genetic basis and repeatability for desiccation resistance in Drosophila melanogaster (Diptera: Drosophilidae).

Dehydration is a stress factor for organisms inhabiting natural habitats where water is scarce. Thus, it may be expected that species facing arid environments will develop mechanisms that maximize resistance to desiccation. Insects are excellent models for studying the effects of dehydration as well as the mechanisms and processes that prevent water loss since the effect of desiccation is greater due to the higher area/volume ratio than larger animals. Even though physiological and behavioral mechanisms to cope with desiccation are being understood, the genetic basis underlying the mechanisms related to variation in desiccation resistance and the context-dependent effect remain unsolved. Here we analyze the genetic bases of desiccation resistance in Drosophila melanogaster and identify candidate genes that underlie trait variation. Our quantitative genetic analysis of desiccation resistance revealed sexual dimorphism and extensive genetic variation. The phenotype-genotype association analyses (GWAS) identified 71 candidate genes responsible for total phenotypic variation in desiccation resistance. Half of these candidate genes were sex-specific suggesting that the genetic architecture underlying this adaptive trait differs between males and females. Moreover, the public availability of desiccation data analyzed on the same lines but in a different lab allows us to investigate the reliability and repeatability of results obtained in independent screens. Our survey indicates a pervasive micro-environment lab-dependent effect since we did not detect overlap in the sets of genes affecting desiccation resistance identified between labs.

Phylogenomics provides insights into the evolution of cactophily and host plant shifts in Drosophila.

Cactophilic species of the Drosophila buzzatii cluster (repleta group) comprise an excellent model group to investigate genomic changes underlying adaptation to extreme climate conditions and host plants. In particular, these species form a tractable system to study the transition from chemically simpler breeding sites (like prickly pears of the genus Opuntia) to chemically more complex hosts (columnar cacti). Here, we report four highly contiguous genome assemblies of three species of the buzzatii cluster. Based on this genomic data and inferred phylogenetic relationships, we identified candidate taxonomically restricted genes (TRGs) likely involved in the evolution of cactophily and cactus host specialization. Functional enrichment analyses of TRGs within the buzzatii cluster identified genes involved in detoxification, water preservation, immune system response, anatomical structure development, and morphogenesis. In contrast, processes that regulate responses to stress, as well as the metabolism of nitrogen compounds, transport, and secretion were found in the set of species that are columnar cacti dwellers. These findings are in line with the hypothesis that those genomic changes brought about key mechanisms underlying the adaptation of the buzzatii cluster species to arid regions in South America.

Oviposition behaviour in Drosophila melanogaster: Genetic and behavioural decoupling between oviposition acceptance and preference for natural fruits.

In phytophagous insects, oviposition behaviour is an important component of habitat selection and, given the multiplicity of genetic and environmental factors affecting its expression, is defined as a complex character resulting from the sum of interdependent traits. Here, we study two components of egg-laying behaviour: oviposition acceptance (OA) and oviposition preference (OP) in Drosophila melanogaster using three natural fruits as resources (grape, tomato and orange) by means of no-choice and two-choice experiments, respectively. This experimental design allowed us to show that the results obtained in two-choice assays (OP) cannot be accounted for by those resulting from no-choice assays (OA). Since the genomes of all lines used are completely sequenced, we perform a genome-wide association study to identify and characterize the genetic underpinnings of these oviposition behaviour traits. The analyses revealed different candidate genes affecting natural genetic variation of both OA and OP traits. Moreover, our results suggest behavioural and genetic decoupling between OA and OP and that egg-laying behaviour is plastic and context-dependent. Such independence in the genetic architectures of OA and OP variation may influence different aspects of oviposition behaviour, including plasticity, canalization, host shift and maintenance of genetic variability, which contributes to the adoption of adaptive strategies during habitat selection.

Extended lifespan and sex-specific fertility loss in cold-acclimated flies of the sibling species Drosophila buzzatii and Drosophila koepferae.

Survival and reproduction are the core elements of Darwinian fitness. In the context of a fixed energy budget, organisms tend to allocate resources in order to maximize one at the expense of the other, in what has been called the lifespan-reproduction trade-off. Reproductive arrest and extended lifespan are common responses to low temperatures in many insects including fruit flies. In this study, we aim to understand the overwintering strategy of two closely-related Drosophila species with contrasting distribution ranges. We compared survival, lifespan, ovarian maturation, and reproductive output (fecundity and fertility) of virgin and mated adults of both Drosophila buzzatii and Drosophila koepferae after long-term cold exposure at dormancy-inducing conditions (10 °C, 10:14 L:D) and controls (25 °C, 12:12 L:D). Virgin flies of D. buzzatii showed the longest lifespan (averaging 102 days) under dormancy-inducing conditions. Cold-induced reproductive arrest preserves reproductive capacity mainly in virgin females that mated after reproductive dormancy, indicating that males were much more susceptible to fertility loss than females, in both species. Notably, females of D. buzzatii were capable of protecting stored sperm from cold damage and produced viable progeny. Even if, in D. buzzatii, fertility of flies mated after the cold-exposure was extremely low, cold temperature likely sterilized D. koepferae males, indicating that cold carry-over effects are stronger for the species with the shorter lifespan. Such species-specific effects of low temperature over fitness likely contributed to the divergence of these closely-related species and to the spread of D. buzzatii into cooler environments.

Spatial and host related genomic variation in partially sympatric cactophagous moth species.

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.

Research gaps and new insights in the evolution of Drosophila seminal fluid proteins.

While the striking effects of seminal fluid proteins (SFPs) on females are fairly conserved among Diptera, most SFPs lack detectable homologues among the SFP repertoires of phylogenetically distant species. How such a rapidly changing proteome conserves functions across taxa is a fascinating question. However, this and other pivotal aspects of SFPs' evolution remain elusive because discoveries on these proteins have been mainly restricted to the model Drosophila melanogaster. Here, we provide an overview of the current knowledge on the inter-specific divergence of the SFP repertoire in Drosophila and compile the increasing amount of relevant genomic information from multiple species. Capitalizing on the accumulated knowledge in D. melanogaster, we present novel sets of high-confidence SFP candidates and transcription factors presumptively involved in regulating the expression of SFPs. We also address open questions by performing comparative genomic analyses that failed to support the existence of many conserved SFPs shared by most dipterans and indicated that gene co-option is the most frequent mechanism accounting for the origin of Drosophila SFP-coding genes. We hope our update establishes a starting point to integrate further data and thus widen the understanding of the intricate evolution of these proteins.

On how to identify a seminal fluid protein: A response to Wigby et al.

The choice of criteria to delimit a group or class is a subjective matter, even though the reasoning, the objectives and the criteria themselves should always be clearly stated. This paper is part of a discussion about the criteria used to identify seminal fluid proteins (SFPs) in Drosophila species. SFPs are proteins that are transferred to females during copulation together with sperm. The only way to ascertain that a protein is an SFP is to prove that it is produced in a male reproductive organ and is found in the female reproductive tract after insemination. Nevertheless, the required methodology is labour-intensive and expensive, and therefore this kind of data is unlikely to be available for many species, precluding comparative and evolutionary studies on the subject. To conduct evolutionary analyses, in a previous study, we capitalized on the accumulated knowledge we have in the model species D. melanogaster to recommend a set of criteria for identifying candidate SFPs in other Drosophila species. Those criteria, based on transcriptomic evidence and in silico predictions from sequences, would allow a good balance between sensitivity (the inclusion of true SFPs) and specificity (the exclusion of false positives). In view of the criticism raised by another group, here we defend our criteria on one hand while accepting there is room for improvement on the other. The results are updated sets of criteria and SFPs that we believe can be useful in future evolutionary studies.

Divergent metabolomic profiles of cold-exposed mature and immature females of tropical versus temperate Drosophila species.

Temperate species, contrary to their tropical counterparts, are exposed not only to thermally variable environments with low temperatures but also to long winters. Different selective pressures may have driven divergent physiological adaptations in closely related species with different biogeographic origins. To survive unfavourable winter conditions, Drosophila species in temperate areas generally undergo a period of reproductive dormancy, associated with a cold-induced cessation of oogenesis and metabolic reorganization. This work aims to compare cold tolerance and metabolic signatures of cold-exposed females exhibiting different reproductive maturity status (mature and immature females) of four Drosophila species from tropical vs. temperate origins. We expected that the capacity for delayed reproduction of immature females could result in the redirection of the energy-related metabolites to be utilized for surviving the cold season. To do so, we studied an array of 45 metabolites using quantitative target GC-MS profiling. Reproductively immature females of temperate species showed the lower CTmin and the faster chill coma recovery time (i.e. the most cold-tolerant group). Principal component analysis captured differences across species, but also between reproductive maturity states. Notably, temperate species exhibited significantly higher levels of glucose, alanine, and gluconolactone than tropical ones. As proline and glycerol showed higher abundances in immature females of temperate species compared to the levels exhibited by the rest of the groups, we reasoned that glucose and alanine could serve as intermediates in the synthesis of these compatible solutes. All in all, our findings suggest that cold-exposed females of temperate species accumulate energy-related and protective metabolites (e.g. glycerol and proline) while delaying reproduction, and that these metabolites are relevant to cold tolerance even at modest concentrations.

Host Plant Adaptation in Cactophilic Species of the Drosophila buzzatii Cluster: Fitness and Transcriptomics.

Host plant shifts in herbivorous insects often involve facing new environments that may speed up the evolution of oviposition behavior, performance-related traits, morphology, and, incidentally, reproductive isolation. In the genus Drosophila, cactophilic species of the repleta group include emblematic species in the study of the evolution of host plant utilization. The South American D. buzzatii and its sibling D. koepferae are a model system for the study of differential host plant use. Although these species exhibit a certain degree of niche overlap, the former breeds primarily on decaying cladodes of Opuntia cacti while D. koepferae main hosts are columnar cacti of the genus Trichocereus. Opuntia sulphurea and Trichocereus terscheckii are among the main hosts in nature. These cacti differ in ecological (spatial and temporal predictability) and chemical characteristics. Particularly relevant is the presence of toxic alkaloids in T. terscheckii. Studies of the effects of these cacti and alkaloids revealed the remarkable impact on oviposition behavior, viability, developmental time, wing morphology, mating success, and developmental stability in both species. Recent whole-genome expression studies showed that expression profiles are massively affected by the rearing cactus, and that the presence of alkaloids is the main factor modulating gene expression in D. buzzatii. Functional enrichment analysis indicated that differentially expressed genes are related to detoxification processes and stress response-though genes involved in development are an important part of the transcriptomic response. The implications of our studies in the evolution of host plant use in the repleta group are discussed.

The influence of developmental environment on courtship song in cactophilic Drosophila.

Closely related species often differ in the signals involved in sexual communication and mate recognition. Determining the factors influencing signal quality (i.e. signal's content and conspicuousness) provides an important insight into the potential pathways by which these interspecific differences evolve. Host specificity could bias the direction of the evolution of sexual communication and the mate recognition system, favouring sensory channels that work best in the different host conditions. In this study, we focus on the cactophilic sibling species Drosophila buzzatii and D. koepferae that have diverged not only in the sensory channel used for sexual communication and mate recognition but also in the cactus species that use as primary hosts. We evaluate the role of the developmental environment in generating courtship song variation using an isofemale line design. Our results show that host environment during development induces changes in the courtship song of D. koepferae males, but not in D. buzzatii males. Moreover, we report for the first time that host rearing environment affects the conspicuousness of courtship song (i.e. song volume). Our results are mainly discussed in the context of the sensory drive hypothesis.

Enhanced fertility and chill tolerance after cold-induced reproductive arrest in females of temperate species of the Drosophila buzzatii complex.

Long-term exposure to low temperatures during adult maturation might decrease fertility after cold recovery as a consequence of carry-over effects on reproductive tissues. This pattern should be more pronounced in tropical than in temperate species as protective mechanisms against chilling injuries are expected to be more effective in the latter. We initially determined the lower thermal thresholds to induce ovarian maturation in four closely related Drosophila species, two inhabiting temperate regions and the other two tropical areas of South America. As expected, only temperate species regularly experience cold-inducing conditions for reproductive arrest during winter in their natural environment. Subsequently, we exposed reproductively arrested and mature females to cold-inducing conditions for reproductive arrest over a long period. Following cold exposure, tropical species exhibited a dramatic fertility decline, irrespective of reproductive status. In contrast, not only were temperate females fecund and fertile but also fertility was superior in females that underwent cold-induced reproductive arrest, suggesting that it might act as a protecting mechanism ensuring fertility after cold recovery. Based on these findings, we decided to evaluate the extent to which reproductive status affects cold tolerance and energy metabolism at low temperature. We found a lower metabolic rate and a higher cold tolerance in reproductively arrested females, although only temperate species attained high levels of chill tolerance. These findings highlight the role of cold-induced reproductive arrest as part of an integrated mechanism of cold adaptation that could potentially contribute to the spread of temperate species into higher latitudes or altitudes.

Transcriptome modulation during host shift is driven by secondary metabolites in desert Drosophila.

High-throughput transcriptome studies are breaking new ground to investigate the responses that organisms deploy in alternative environments. Nevertheless, much remains to be understood about the genetic basis of host plant adaptation. Here, we investigate genome-wide expression in the fly Drosophila buzzatii raised in different conditions. This species uses decaying tissues of cactus of the genus Opuntia as primary rearing substrate and secondarily, the necrotic tissues of the columnar cactus Trichocereus terscheckii. The latter constitutes a harmful host, rich in mescaline and other related phenylethylamine alkaloids. We assessed the transcriptomic responses of larvae reared in Opuntia sulphurea and T. terscheckii, with and without the addition of alkaloids extracted from the latter. Whole-genome expression profiles were massively modulated by the rearing environment, mainly by the presence of T. terscheckii alkaloids. Differentially expressed genes were mainly related to detoxification, oxidation-reduction and stress response; however, we also found genes involved in development and neurobiological processes. In conclusion, our study contributes new data onto the role of transcriptional plasticity in response to alternative rearing environments.

Parental age influences developmental stability of the progeny in Drosophila.

The stochastic nature of biochemical processes is a source of variability that influences developmental stability. Developmental instability (DI) is often estimated through fluctuating asymmetry (FA), a parameter that deals with within-individual variation in bilateral structures. A relevant goal is to shed light on how environment, physiology and genotype relate to DI, thus providing a more comprehensive view of organismal development. Using Drosophila melanogaster isogenic lines, we investigated the effect of parental age, parental diet and offspring heterozygosity on DI. In this work, we have uncovered a clear relationship between parental age and offspring asymmetry. We show that asymmetry of the progeny increases concomitantly with parental age. Moreover, we demonstrate that enriching the diet of parents mitigates the effect of age on offspring symmetry. We show as well that increasing the heterozygosity of the progeny eliminates the effect of parental age on offspring symmetry. Taken together, our results suggest that diet, genotype and age of the parents interact to determine offspring DI in wild populations. These findings provide us with an avenue to understand the mechanisms underlying DI.

City puzzles: Does urban land scape affect genetic population structure in Aedes aegypti?

Cities usually offer a suitable environment for the dengue vector Aedes aegypti, providing oviposition sites, accessibility to human hosts and nectar meals. However, large urban centres are highly heterogeneous environments, forming a patched landscape that could affect Ae. aegypti population dynamics and dispersal. Here, we performed a genome-wide analysis using Rad-seq data from 99 Ae. aegypti specimens collected in three areas within Buenos Aires city with varying levels of urbanization/land use: highly urbanized Area 1, intermediate Area 2 and poorly urbanized Area 3. We found an inverse association between urbanization levels and spatial genetic structure. Populations from highly urbanized Area 1 did not present genetic structure whereas two and three clusters were detected in Areas 2 and 3, respectively. In the case of Area 3, initial analyses showed separation in clusters was mostly due to elevated consanguinity within sites although three clusters were still detected after closely related individuals were discarded. Mosquitoes around each site displayed a high degree of isolation, evidencing a close dependence between the vector and human dwellings. Interestingly, specimens from distant boroughs (within the limits of the city) and the city's outskirts formed a single cluster with inner city sites (Area 1), highlighting the role of passive transport in shaping population structure. Genetic distances were poorly correlated with geographic distances in Buenos Aires, suggesting a stronger influence of passive than active dispersal on population structure. Only Area 2 displayed a significant isolation-by-distance pattern (p = 0.046), with males dispersing more than females (p = 0.004 and p = 0.016, respectively). Kinship analyses allowed us to detect full-siblings located 1.5 km apart in Area 1, which could be due to an extreme event of active female dispersal. Effective population size was higher in Area 2 confirming that cemeteries represent highly favourable environments for Ae. aegypti and need to be specifically targeted. Our results suggest that control programs should take into account urban landscape heterogeneity in order to improve vector control.

Ortholog genes from cactophilic Drosophila provide insight into human adaptation to hallucinogenic cacti.

Cultural transformations of lifestyles and dietary practices have been key drivers of human evolution. However, while most of the evidence of genomic adaptations is related to the hunter-gatherer transition to agricultural societies, little is known on the influence of other major cultural manifestations. Shamanism is considered the oldest religion that predominated throughout most of human prehistory and still prevails in many indigenous populations. Several lines of evidence from ethno-archeological studies have demonstrated the continuity and importance of psychoactive plants in South American cultures. However, despite the well-known importance of secondary metabolites in human health, little is known about its role in the evolution of ethnic differences. Herein, we identified candidate genes of adaptation to hallucinogenic cactus in Native Andean populations with a long history of shamanic practices. We used genome-wide expression data from the cactophilic fly Drosophila buzzatii exposed to a hallucinogenic columnar cactus, also consumed by humans, to identify ortholog genes exhibiting adaptive footprints of alkaloid tolerance. Genomic analyses in human populations revealed a suite of ortholog genes evolving under recent positive selection in indigenous populations of the Central Andes. Our results provide evidence of selection in genetic variants related to alkaloids toxicity, xenobiotic metabolism, and neuronal plasticity in Aymara and Quechua populations, suggesting a possible process of gene-culture coevolution driven by religious practices.

Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects.

Hypogeococcus pungens is a species complex native to southern South America that is composed of at least five putative species, each one specialized in the use of different host plants. Two of these undescribed species were registered as invasive in Central and North America: Hyp-C is a cactophagous mealybug that became an important pest that threatens endemic cactus species in Puerto Rico, and Hyp-AP feeds on Amaranthaceae and Portulacaceae hosts, but does not produce severe damage to the host plants. We quantified genomic variation and investigated the demographic history of both invasive species by means of coalescent-based simulations using high throughput sequencing data. We also evaluated the incidence of host plant infestation produced by both species and used an ecological niche modeling approach to assess potential distribution under current and future climatic scenarios. Our genetic survey evinced the footprints of strong effective population size reduction and signals of genetic differentiation among populations within each species. Incidence of plant attacks varied between species and among populations within species, with some host plant species preferred over others. Ecological niche modeling suggested that under future climatic scenarios both species would expand their distribution ranges in Puerto Rico. These results provide valuable information for the design of efficient management and control strategies of the Puerto Rican cactus pest and shed light on the evolutionary pathways of biological invasions.

Geographic host use variability and host range evolutionary dynamics in the phytophagous insect Apagomerella versicolor (Cerambycidae).

The high diversity of phytophagous insects has been explained by the tendency of the group towards specialization; however, generalism may be advantageous in some environments. The cerambycid Apagomerella versicolor exhibits intraspecific geographical variation in host use. In northern Argentina it is highly specialized on the herb Pluchea sagittalis (Asteraceae), while in central and southern areas it uses seven Asteraceae species. To study host species geographical variation from ecological and evolutionary perspectives, we investigated field host availability and use across a wide latitudinal range, and performed laboratory studies on insect oviposition preference and larval performance and mitochondrial DNA (mtDNA) variation in a phylogeographical framework. Geographic variation in host use was unrelated to host availability but was highly associated with laboratory oviposition preference, larval performance, and mtDNA variation. Genetic studies revealed three geographic races of A. versicolor with gene flow restriction and recent geographic expansion. Trophic generalism and oligophagy within A. versicolor seem to have evolved as adaptations to seasonal and spatial unavailability of the preferred host P. sagittalis in cooler areas of the species' geographic range. No single genotype is successful in all environments; specialization may be advantageous in environments with uniform temporal and spatial host availability, while being a trophic generalist may provide an adaptive advantage in host-constrained environments.

Cis-regulatory variation in the shavenbaby gene underlies intraspecific phenotypic variation, mirroring interspecific divergence in the same trait.

Despite considerable progress in recent decades in dissecting the genetic causes of natural morphological variation, there is limited understanding of how variation within species ultimately contributes to species differences. We have studied patterning of the non-sensory hairs, commonly known as "trichomes," on the dorsal cuticle of first-instar larvae of Drosophila. Most Drosophila species produce a dense lawn of dorsal trichomes, but a subset of these trichomes were lost in D. sechellia and D. ezoana due entirely to regulatory evolution of the shavenbaby (svb) gene. Here, we describe intraspecific variation in dorsal trichome patterns of first-instar larvae of D. virilis that is similar to the trichome pattern variation identified previously between species. We found that a single large effect QTL, which includes svb, explains most of the trichome number difference between two D. virilis strains and that svb expression correlates with the trichome difference between strains. This QTL does not explain the entire difference between strains, implying that additional loci contribute to variation in trichome numbers. Thus, the genetic architecture of intraspecific variation exhibits similarities and differences with interspecific variation that may reflect differences in long-term and short-term evolutionary processes.

A quantitative genetic study of starvation resistance at different geographic scales in natural populations of Drosophila melanogaster.

Food shortage is a stress factor that commonly affects organisms in nature. Resistance to food shortage or starvation resistance (SR) is a complex quantitative trait with direct implications on fitness. However, surveys of natural genetic variation in SR at different geographic scales are scarce. Here, we have measured variation in SR in sets of lines derived from nine natural populations of Drosophila melanogaster collected in western Argentina. Our study shows that within population variation explained a larger proportion of overall phenotypic variance (80%) than among populations (7·2%). We also noticed that an important fraction of variation was sex-specific. Overall females were more resistant to starvation than males; however, the magnitude of the sexual dimorphism (SD) in SR varied among lines and explained a significant fraction of phenotypic variance in all populations. Estimates of cross-sex genetic correlations suggest that the genetic architecture of SR is only partially shared between sexes in the populations examined, thus, facilitating further evolution of the SD.

Geographic patterns of inversion polymorphism in the second chromosome of the cactophilic Drosophila buzzatii from northeastern Argentina.

The inversion polymorphisms of the cactophilic Drosophila buzzatti Patterson and Wheeler (Diptera: Drosophilidae) were studied in new areas of its distribution in Argentina. A total of thirty-eight natural populations, including 29 from previous studies, were analyzed using multiple regression analyses. The results showed that about 23% of total variation was accounted for by a multiple regression model in which only altitude contributed significantly to population variation, despite the fact that latitude and longitude were also included in the model. Also, inversion frequencies exhibited significant associations with mean annual temperature, precipitation, and atmospheric pressure. In addition, expected heterozygosity exhibited a negative association with temperature and precipitation and a positive association with atmospheric pressure. The close similarity of the patterns detected in this larger dataset to previous reports is an indication of the stability of the clines. Also, the concurrence of the clines detected in Argentina with those reported for colonizing populations of Australia suggests the involvement of natural selection as the main mechanism shaping inversion frequencies in D. buzzatii.