Cactophilic Drosophila in South America: a model for evolutionary studies

Evidences of differential methylation in the genome during development in the cactophilic Drosophila species

DNA methylation with 5-methylcytosine (5mC) has been reported in the genome of several eukaryotes, with marked differences between vertebrates and invertebrates. DNA methylation is poorly understood as its role in evolution in insects. Drosophila gouveai (cluster Drosophila buzzatii) presents larvae that develop obligatorily in necrotic tissues of cacti in nature, with the distribution of populations in South America, and plasticity of phenotypes in insect-plant interaction. We characterize organisms at developmental stages and analyze variations at multiple methylation-sensitive loci in pupae, and adult flies using methylation sensitive amplification polymorphism. We obtained 326 loci with CCGG targets in the genome of D. gouveai. Genomic regions with molecular lengths from 100 to 700 pb were most informative about methylation states. Multiple loci show differences in methylation-sensitive sites (MSL) concerning developmental stages, such as in pupae (MSL = 40), female reproductive tissue (MSL = 76), and male reproductive tissues (MSL = 58). Our results are the first evidence of genome-wide methylation in D. gouveai organisms.

Systematics and spatio-temporal evolutionary patterns of the flavopilosa group of Drosophila (Diptera, Drosophilidae)

The Drosophila flavopilosa group comprises morphologically cryptic species that are ecologically restricted to feeding, breeding and ovipositing on flowers of Cestrum and Sessea (Solanaceae). Previous studies confirmed the monophyly of the group and the success of DNA barcoding in identifying a subset of its species, but several others remain yet to be evaluated. Furthemore, the taxonomy of the group remains incomplete, with only nine of the 17 species assigned to subgroups. Here, we accessed the phylogenetic relationships and spatio-temporal evolutionary patterns of the flavopilosa group based on a mitochondrial and two nuclear genes, providing the first molecular support to the subdivision of the group and suggesting a new taxonomic scheme for its species. Barcoding proved to be an effective tool, as all species were reciprocally monophyletic and different analyses of species delimitation yielded congruent results. The close relationship of D. flavopilosa with D. cestri and D. cordeiroi was strongly supported, suggesting that the latter should be placed in the flavopilosa subgroup together with the first. Furthermore, D. mariaehelenae was positioned as sister to D. incompta, supporting its inclusion in the nesiota subgroup. Despite new taxonomic assignments, the synapomorphic status of the diagnostic characters proposed for both subgroups was supported. Based on them, each of the remaining species were placed into one of both subgroups. Divergence time estimates suggest that their diversification coincided with the divergence of Sessea and Cestrum, providing an interesting case of coevolution.

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.

Evolution of the odorant-binding protein gene family in Drosophila

Odorant-binding proteins (OBPs) are encoded by a gene family involved in the perception of olfactory signals in insects. This chemosensory gene family has been advocated as a candidate to mediate host preference and host shifts in insects, although it also participates in other physiological processes. Remarkable differences in the OBP gene repertoire have been described across insect groups, suggesting an accelerated gene turnover rate. The genus Drosophila, is a valuable resource for ecological genomics studies since it comprises groups of ecologically diverse species and there are genome data for many of them. Here, we investigate the molecular evolution of this chemosensory gene family across 19 Drosophila genomes, including the melanogaster and repleta species groups, which are mostly associated with rotting fruit and cacti, respectively. We also compared the OBP repertoire among the closely related species of the repleta group, associated with different subfamilies of Cactaceae that represent disparate chemical challenges for the flies. We found that the gene family size varies widely between species, ranging from 39 to 54 candidate OBPs. Indeed, more than 54% of these genes are organized in clusters and located on chromosomes X, 2, and 5, with a distribution conserved throughout the genus. The family sizes in the repleta group and D. virilis (virilis-repleta radiation) were smaller than in the melanogaster group. We tested alternative evolutionary models for OBP family size and turnover rates based on different ecological scenarios. We found heterogeneous gene turnover rates (GR) in comparisons involving columnar cactus specialists, prickly pear specialists, and fruit dwellers lineages, and signals of rapid molecular evolution compatible with positive selection in specific OBP genes. Taking ours and previous results together, we propose that this chemosensory gene family is involved in host adaptation and hypothesize that the adoption of the cactophilic lifestyle in the repleta group accelerated the evolution of members of the family.

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.

A holocenic and dynamic hybrid zone between two cactophilic Drosophila species in a coastal lowland plain of the Brazilian Atlantic Forest

Hybridization and introgression are processes that contribute to shaping biological diversity. The factors promoting the formation of these processes are multiples but poorly explored in a biogeographical and ecological context. In the southeast coastal plain of the Brazilian Atlantic Forest, a hybrid zone was described between two closely related cactophilic species, Drosophila antonietae and D. serido. Here, we revisited and analysed specimens from this hybrid zone to evaluate its temporal and spatial dynamic. We examined allopatric and sympatric populations of the flies using independent sources of data such as mitochondrial and nuclear sequences, microsatellite loci, morphometrics of wings and male genitalia, and climatic niche models. We also verified the emergence of the flies from necrotic tissues of collected cacti to verify the role of host association for the population dynamics. Our results support the existence of a hybrid zone due to secondary contact and limited to the localities where the two species are currently in contact. Furthermore, we detected asymmetric bidirectional introgression and the maintenance of the species integrity, ecological association and morphological characters, suggesting selection and limited introgression. Considering our paleomodels, probably this hybrid zone is recent and the contact occurred during the Holocene to the present day, favoured by range expansion of their populations due to expansion of open and dry areas in eastern South America during palaeoclimatic and geomorphological events.

Evolution of male genitalia in the Drosophila repleta species group (Diptera: Drosophilidae)

The Drosophila repleta group comprises more than one hundred species that inhabit several environments in the Neotropics and use different hosts as rearing and feeding resources. Rather homogeneous in their external morphology, they are generally distinguished by the male genitalia, seemingly their fastest evolving morphological trait, constituting an excellent model to study patterns of genital evolution in the context of a continental adaptive radiation. Although much is known about the evolution of animal genitalia at population level, surveys on macroevolutionary scale of this phenomenon are scarce. This study used a suite of phylogenetic comparative methods to elucidate the macroevolutionary patterns of genital evolution through deep time and large continental scales. Our results indicate that male genital size and some aspects of shape have been evolving by speciational evolution, probably due to the microevolutionary processes involved in species mate recognition. In contrast, several features of the aedeagus shape seemed to have evolved in a gradual fashion, with heterogeneous evolutionary phenotypic rates among clades. In general, the tempo of the evolution of aedeagus morphology was constant from the origin of the group until the Pliocene, when it accelerated in some clades that diversified mainly in this period. The incidence of novel ecological conditions in the tempo of aedeagus evolution and the relationship between species mate recognition and speciation in the Drosophila repleta group are discussed.

Plastic Variation in the Phyletic Lineages of Cactophilic Drosophila meridionalis and Relation to Hosts as Potential for Diversification

The insect/plant interaction is known to be a trigger for diversification and even speciation. Experimental analyses on fitness traits and phenotypic variation using alternative host sites have been performed to understand the process of diversification relative to insect/plant interactions. For cactophilic species of Drosophila, the speciation process is considered an adaptive radiation in response to the exploration of species of the Cactaceae as breeding and feeding sites. In this work, we analyzed life history and morphological traits in individuals from two phyletic lineages (Evolutionarily Significant Units ESU) of the cactophilic species Drosophila meridionalis (Wasserman 1962) (Diptera: Drosophilidae) raised from media prepare. The characters analyzed corresponded to viability, developmental time, and four morphological measurements. The experiments were performed in a semi-natural medium prepared with fermenting tissues of the natural hosts, Cereus hildmaniannus and Opuntia monacantha. Viability, development time, and three morphological measurements were influenced by lineage, suggesting differentiation between the lineages. However, in O. monacantha, the mean viability was greater (~15%) and development time was longer (~336 h) than in C. hildmaniannus (~11% and ~301 h, respectively). Only the developmental time was significantly affected by the host cactus. In general, ESU group A had better values than ESU group BC for the evaluated traits. This finding suggested differentiation between the two lineages and different plastic responsiveness to the contrasting environments of the hosts, and that C. hildmaniannus may be a relatively stressful environment for the larvae, as for other Drosophila species.

Phylogenetic Relationships and Evolutionary Trends in the Cactus Family

Members of the cactus family are keystone species of arid and semiarid biomes in the Americas, as they provide shelter and resources to support other members of ecosystems. Extraordinary examples are the several species of flies of the genus Drosophila that lay eggs and feed in their rotting stems, which provide a model system for studying evolutionary processes. Although there is significant progress in understanding the evolution of Drosophila species, there are gaps in our knowledge about the cactus lineages hosting them. Here, we review the current knowledge about the evolution of Cactaceae, focusing on phylogenetic relationships and trends revealed by the study of DNA sequence data. During the last several decades, the availability of molecular phylogenies has considerably increased our understanding of the relationships, biogeography, and evolution of traits in the family. Remarkably, although succulent cacti have very low growth rates and long generation times, they underwent some of the fastest diversifications observed in the plant kingdom, possibly fostered by strong ecological interactions. We have a better understanding of the reproductive biology, population structure and speciation mechanisms in different clades. The recent publication of complete genomes for some species has revealed the importance of phenomena such as incomplete lineage sorting. Hybridization and polyploidization are common in the family, and have been studied using a variety of phylogenetic methods. We discuss potential future avenues for research in Cactaceae, emphasizing the need of a concerted effort among scientists in the Americas, together with the analyses of data from novel sequencing techniques.

Genetics and Distribution Modeling: The Demographic History of the Cactophilic Drosophila buzzatii Species Cluster in Open Areas of South America

Biodiversity is the result of historical and recurrent events acting on populations and species. The Drosophila buzzatii species cluster is distributed along a diagonal of open areas in South America. Combining genetic analyses with species distribution models we evaluated the influence of climatic changes in the demography history of this cluster. We performed a Bayesian Skyline analysis and reconstructed the ancestral areas based on mitochondrial cytochrome oxidase I (mtCOI) gene. We modeled the current distribution and projected it to past (mid-Holocene and Last Glacial Maximum) and future. Our results demonstrate that climate change plays a critical role in historical demography and in defining the current and future geographic ranges of these species. The inter-Andean dry valleys and the Chiquitano Seasonally Dry Tropical Forests (SDTF) in Bolivia are considered the ancestral area for the D. buzzatii cluster. From this area, the migration route was through a west-east corridor from central Andes, throughout Bolivia, Paraguay and Argentina toward eastern and northeastern Brazil, along the edges of rain forest. The responses of the species to the climatic changes differ from the dominant assumptions of expansion during dry/cold weather and contraction during wet/warm weather that characterized the glacial cycles. We suggest that the influence of ecology and the potential responsiveness of each taxon to the environmental dynamics should be considered as well. Predictions for the future suggest a decline of suitable areas for the cluster, threatening biodiversity of these habitats. This work showed the importance of an integrative analysis of genetics and geography information to improve the inferences about demographic history hypotheses for the cluster.

What does mitogenomics tell us about the evolutionary history of the Drosophila buzzatii cluster (repleta group)?

The Drosophila repleta group is an array of more than 100 species endemic to the “New World”, many of which are cactophilic. The ability to utilize decaying cactus tissues as breeding and feeding sites is a key aspect that allowed the successful diversification of the repleta group in American deserts and arid lands. Within this group, the Drosophila buzzatii cluster is a South American clade of seven closely related species in different stages of divergence, making them a valuable model system for evolutionary research. Substantial effort has been devoted to elucidating the phylogenetic relationships among members of the D. buzzatii cluster, including molecular phylogenetic studies that have generated ambiguous results where different tree topologies have resulted dependent on the kinds of molecular marker used. Even though mitochondrial DNA regions have become useful markers in evolutionary biology and population genetics, none of the more than twenty Drosophila mitogenomes assembled so far includes this cluster. Here, we report the assembly of six complete mitogenomes of five species: D. antonietae, D. borborema, D. buzzatii, two strains of D. koepferae and D. seriema, with the aim of revisiting phylogenetic relationships and divergence times by means of mitogenomic analyses. Our recovered topology using complete mitogenomes supports the hypothesis of monophyly of the D. buzzatii cluster and shows two main clades, one including D. buzzatii and D. koepferae (both strains), and the other containing the remaining species. These results are in agreement with previous reports based on a few mitochondrial and/or nuclear genes, but conflict with the results of a recent large-scale nuclear phylogeny, indicating that nuclear and mitochondrial genomes depict different evolutionary histories.

Pupal emergence pattern in cactophilic Drosophila and the effect of host plants

Drosophila buzzatii and D. koepferae are sibling cactophilic species. The former breeds primarily on prickly pears (genus Opuntia) whereas the latter breeds on columnar cacti of the genera Cereus and Trichocereus, although with certain degree of niche overlapping. We examined the interspecific differences in diurnal temporal patterns of adult emergence from puparia and evaluated whether this behavior is affected by rearing in the different cactus hosts available in nature. We detected important host-dependent genetic variation for this trait differentially affecting the emergence schedule of these species. Diurnal pattern of emergence time was directly correlated with developmental time and negatively correlated with adult wing size, suggesting that early emergences are at least indirectly correlated with increased fitness. We discussed our results in terms of their putative effects on fitness and the genetic-metabolic pathways that would be presumably affected by host's nutritional-chemical differences.

The Calibrated Phylogeny of the Drosophila fasciola Subgroup (D. repleta Group Wasserman) Indicates Neogene Diversification of Its Internal Branches

The species of the Drosophila fasciola subgroup Wasserman represent the dominant section of the Drosophila repleta group Wasserman in the American rainforests and have a broad geographical distribution in the New World. However, despite of its wide range, the D. fasciola subgroup is one of the most overlooked D. repleta subgroups. Here, we report a molecular phylogenetic analysis focused on the D. fasciola subgroup using two mitochondrial [cytochrome oxidase subunit I (COI), cytochrome oxidase subunit II (COII)] and two nuclear [elongation factor-1alpha F1 (EF-alphaF1) and transformer (tra)] genes. Overall, we found that this subgroup is a monophyletic taxon, subdivided into two main internal branches: named Fas1 and Fas2 clades. The diversification of these clades is estimated to have begun in the middle Miocene, around 12 Ma [95% high posterior density (HPD) 9.0-15 Ma], and might be associated with the colonization of South America by Central America populations after the closure of Isthmus of Panama due to the temporal congruence between these events. The terminal branches had their origins estimated to be in the Pliocene or the Plio-Pleistocene transition. For the later estimates, both the geomorphological influences and the climatic oscillations of the Pleistocene may have played a role in shaping the diversification of the D. fasciola group.

Dissecting the Satellite DNA Landscape in Three Cactophilic Drosophila Sequenced Genomes

Eukaryote genomes are replete with repetitive DNAs. This class includes tandemly repeated satellite DNAs (satDNA) which are among the most abundant, fast evolving (yet poorly studied) genomic components. Here, we used high-throughput sequencing data from three cactophilic Drosophila species, D. buzzatii, D. seriema, and D. mojavensis, to access and study their whole satDNA landscape. In total, the RepeatExplorer software identified five satDNAs, three previously described (pBuM, DBC-150 and CDSTR198) and two novel ones (CDSTR138 and CDSTR130). Only pBuM is shared among all three species. The satDNA repeat length falls within only two classes, between 130 and 200 bp or between 340 and 390 bp. FISH on metaphase and polytene chromosomes revealed the presence of satDNA arrays in at least one of the following genomic compartments: centromeric, telomeric, subtelomeric, or dispersed along euchromatin. The chromosomal distribution ranges from a single chromosome to almost all chromosomes of the complement. Fiber-FISH and sequence analysis of contigs revealed interspersion between pBuM and CDSTR130 in the microchromosomes of D. mojavensis Phylogenetic analyses showed that the pBuM satDNA underwent concerted evolution at both interspecific and intraspecific levels. Based on RNA-seq data, we found transcription activity for pBuM (in D. mojavensis) and CDSTR198 (in D. buzzatii) in all five analyzed developmental stages, most notably in pupae and adult males. Our data revealed that cactophilic Drosophila present the lowest amount of satDNAs (1.9-2.9%) within the Drosophila genus reported so far. We discuss how our findings on the satDNA location, abundance, organization, and transcription activity may be related to functional aspects.

The role of courtship song in female mate choice in South American Cactophilic Drosophila

Courtship songs have undergone a spectacular diversification in the Drosophila buzzatii cluster. Accordingly, it has been suggested that sexual selection has played a significant role in promoting rapid diversification, reproductive isolation and speciation. However, there is no direct evidence (i.e., song playback experiments with wingless males) supporting this tenet. Moreover, several studies have showed that the courtship song in the genus Drosophila is not always used in female mate choice decisions, nor plays the same role when it is taken into account. In this vein, we use an approach that combines manipulative and playback experiments to explore the importance and the role of courtship song in female mate choice in four species of the D. buzzatii cluster and one species of the closely related D. martensis cluster for outgroup comparison. We also investigate the importance of courtship song in sexual isolation in sympatry between the only semi-cosmopolitan species, D. buzzatii, and the other species of the D. buzzatii cluster. Our study revealed that the courtship song is used by females of the D. buzzatii cluster as a criterion for male acceptance or influences the speed with which males are chosen. In contrast, we showed that this characteristic is not shared by D. venezolana, the representative species of the D. martensis cluster. We also found that the studied species of the D. buzzatii cluster differ in the role that conspecific and heterospecific songs have in female mate choice and in sexual isolation. Our findings support the hypothesis that divergence in female preferences for courtship songs has played a significant role in promoting rapid diversification and reproductive isolation in the D. buzzatii cluster. However, evidence from D. venezolana suggests that the use of the courtship song in female mate choice is not a conserved feature in the D. buzzatii complex.

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.

Demographic Structure and Evolutionary History of Drosophila ornatifrons (Diptera, Drosophilidae) from Atlantic Forest of Southern Brazil

Drosoph1la ornatifrons of the guarani group (Diptera: Drosophilidae) is found mainly in humid areas of the Atlantic Forest biome, especially in the southern region of Brazil. Historical and contemporary fragmentation events influenced species diversity and distribution in this biome, although the role of paleoclimatic and paleogeographic events remain to be verified. The objective of the present study was to evaluate the demographic structure of D. ornatifrons from collection sites that are remnants of Atlantic Forest in southern Brazil, in order to contribute to the understanding of the processes that affected the patterns of genetic variability in this species. To achieve this goal, we sequenced 51 individuals from nine localities and 64 individuals from six localities for the mitochondrial genes Cytochrome Oxidase I and II, respectively. Our results indicate that D. ornatifrons may have experienced a demographic expansion event from the southernmost locations of its distribution, most likely from those located next to the coast and in fragments of Atlantic Forest inserted in the Pampa biome (South 2 group), towards the interior (South 1 group). This expansion probably started after the last glacial maximum, between 20,000 and 18,000 years ago, and was intensified near the Pleistocene-Holocene transition, around 12,000 years ago, when temperature started to rise. In this work we discuss how the haplotypes found barriers to gene flow and dispersal, influenced by the biogeographic pattern of Atlantic Forest.

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.

Volatile codes: Correlation of olfactory signals and reception in Drosophila-yeast chemical communication

Drosophila have evolved strong mutualistic associations with yeast communities that best support their growth and survival, resulting in the development of novel niches. It has been suggested that flies recognize their cognate yeasts primarily based on the rich repertoire of volatile organic compounds (VOCs) derived from the yeasts. Thus, it remained an exciting avenue to study whether fly spp. detect and discriminate yeast strains based on odor alone, and if so, how such resolution is achieved by the olfactory system in flies. We used two fly species known to exploit different niches and harboring different yeasts, D. suzukii (a pest of fresh fruit) and D. melanogaster (a saprophytic fly and a neurogenetic model organism). We initially established the behavioral preference of both fly species to six Drosophila-associated yeasts; then chemically analyzed the VOC profile of each yeast which revealed quantitative and qualitative differences; and finally isolated and identified the physiologically active constituents from yeast VOCs for each drosophilid that potentially define attraction. By employing chemical, behavioral, and electrophysiological analyses, we provide a comprehensive portrait of the olfactory neuroethological correlates underlying fly-yeast coadaptation in two drosophilids with distinct habitats.

Genomics of ecological adaptation in cactophilic Drosophila

Cactophilic Drosophila species provide a valuable model to study gene–environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii–D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii–D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

Identification of differentially expressed genes in female Drosophila antonietae and Drosophila meridionalis in response to host cactus odor

Background

Studies of insect-plant interactions have provided critical insights into the ecology and evolution of adaptive processes within and among species. Cactophilic Drosophila species have received much attention because larval development occurs in the necrotic tissues of cacti, and both larvae and adults feed on these tissues. Such Drosophila-cactus interactions include effects of the host plant on the physiology and behavior of the flies, especially so their nutritional status, mating condition and reproduction. The aim of this work was to compare the transcriptional responses of two species, Drosophila antonietae and Drosophila meridionalis, and identify genes potentially related to responses to odors released by their host cactus, Cereus hildmannianus. The two fly species are sympatric in most of their populations and use this same host cactus in nature.

Results

We obtained 47 unique sequences (USs) for D. antonietae in a suppression subtractive hybridization screen, 30 of these USs had matches with genes predicted for other Drosophila species. For D. meridionalis we obtained 81 USs, 46 of which were orthologous with genes from other Drosophila species. Functional information (Gene Ontology) revealed that these differentially expressed genes are related to metabolic processes, detoxification mechanisms, signaling, response to stimuli, and reproduction. The expression of 13 genes from D. meridionalis and 12 from D. antonietae were further analyzed by quantitative real time-PCR, showing that four genes were significantly overexpressed in D. antonietae and six in D. meridionalis.

Conclusions

Our results revealed the differential expression of genes related to responses to odor stimuli by a cactus, in two associated fly species. Although the majority of activated genes were similar between the two species, we also observed that certain metabolic pathways were specifically activated, especially those related to signaling pathways and detoxification mechanisms. The activation of these genes may reflect different metabolic pathways used by these flies in their interaction with this host cactus. Our findings provide insight into how the use of C. hildmannianus may have arisen independently in the two fly species, through genetic differentiation in metabolic pathways to effectively explore this cactus as a host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0191-2) contains supplementary material, which is available to authorized users.

Differences in tolerance to host cactus alkaloids in Drosophila koepferae and D. buzzatii

The evolution of cactophily in the genus Drosophila was a major ecological transition involving over a hundred species in the Americas that acquired the capacity to cope with a variety of toxic metabolites evolved as feeding deterrents in Cactaceae. D. buzzatii and D. koepferae are sibling cactophilic species in the D. repleta group. The former is mainly associated with the relatively toxic-free habitat offered by prickly pears (Opuntia sulphurea) and the latter has evolved the ability to use columnar cacti of the genera Trichocereus and Cereus that contain an array of alkaloid secondary compounds. We assessed the effects of cactus alkaloids on fitness-related traits and evaluated the ability of D. buzzatii and D. koepferae to exploit an artificial novel toxic host. Larvae of both species were raised in laboratory culture media to which we added increasing doses of an alkaloid fraction extracted from the columnar cactus T. terschekii. In addition, we evaluated performance on an artificial novel host by rearing larvae in a seminatural medium that combined the nutritional quality of O. sulphurea plus amounts of alkaloids found in fresh T. terschekii. Performance scores in each rearing treatment were calculated using an index that took into account viability, developmental time, and adult body size. Only D. buzzatii suffered the effects of increasing doses of alkaloids and the artificial host impaired viability in D. koepferae, but did not affect performance in D. buzzatii. These results provide the first direct evidence that alkaloids are key determinants of host plant use in these species. However, the results regarding the artificial novel host suggest that the effects of alkaloids on performance are not straightforward as D. koepferae was heavily affected. We discuss these results in the light of patterns of host plan evolution in the Drosophila repleta group.

Contrasting patterns of within-species morphological variation in two cactophilic Drosophila species (Diptera: Drosophilidae)

In this work, we investigated the morphological variation of the intromittent male copulatory organ (aedeagus) of specimens from natural populations of two cactophilic Drosophila species distributed in the southeast region of Brazil, Drosophila gouveai Tidon-Sklorz & Sene and Drosophila antonietae Tidon-Sklorz & Sene. It was explored how the within-species variability is arranged for both species, considering their historical and ecological features. Our results showed two distinct aedeagal morphologies for these species, and differences within species were observed only in D. gouveai as specimens could be distinguished by their population origin. In contrast, after size discrepancies correction, this feature was not detected in D. antonietae. The contrasting patterns of intraspecific variation, together with the other features exhibited by these two species, are most likely to be explained by differences in the historical host plant association and distribution and in demographic events, which determined the evolutionary history of these two South American cactophilic Drosophila species.

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.

Segmental duplication, microinversion, and gene loss associated with a complex inversion breakpoint region in Drosophila

Chromosomal inversions are usually portrayed as simple two-breakpoint rearrangements changing gene order but not gene number or structure. However, increasing evidence suggests that inversion breakpoints may often have a complex structure and entail gene duplications with potential functional consequences. Here, we used a combination of different techniques to investigate the breakpoint structure and the functional consequences of a complex rearrangement fixed in Drosophila buzzatii and comprising two tandemly arranged inversions sharing the middle breakpoint: 2m and 2n. By comparing the sequence in the breakpoint regions between D. buzzatii (inverted chromosome) and D. mojavensis (noninverted chromosome), we corroborate the breakpoint reuse at the molecular level and infer that inversion 2m was associated with a duplication of a ~13 kb segment and likely generated by staggered breaks plus repair by nonhomologous end joining. The duplicated segment contained the gene CG4673, involved in nuclear transport, and its two nested genes CG5071 and CG5079. Interestingly, we found that other than the inversion and the associated duplication, both breakpoints suffered additional rearrangements, that is, the proximal breakpoint experienced a microinversion event associated at both ends with a 121-bp long duplication that contains a promoter. As a consequence of all these different rearrangements, CG5079 has been lost from the genome, CG5071 is now a single copy nonnested gene, and CG4673 has a transcript ~9 kb shorter and seems to have acquired a more complex gene regulation. Our results illustrate the complex effects of chromosomal rearrangements and highlight the need of complementing genomic approaches with detailed sequence-level and functional analyses of breakpoint regions if we are to fully understand genome structure, function, and evolutionary dynamics.

Variations on a theme: diversification of cuticular hydrocarbons in a clade of cactophilic Drosophila

Background

We characterized variation and chemical composition of epicuticular hydrocarbons (CHCs) in the seven species of the Drosophila buzzatii cluster with gas chromatography/mass spectrometry. Despite the critical role of CHCs in providing resistance to desiccation and involvement in communication, such as courtship behavior, mating, and aggregation, few studies have investigated how CHC profiles evolve within and between species in a phylogenetic context. We analyzed quantitative differences in CHC profiles in populations of the D. buzzatii species cluster in order to assess the concordance of CHC differentiation with species divergence.

Results

Thirty-six CHC components were scored in single fly extracts with carbon chain lengths ranging from C₂₉ to C₃₉, including methyl-branched alkanes, n-alkenes, and alkadienes. Multivariate analysis of variance revealed that CHC amounts were significantly different among all species and canonical discriminant function (CDF) analysis resolved all species into distinct, non-overlapping groups. Significant intraspecific variation was found in different populations of D. serido suggesting that this taxon is comprised of at least two species. We summarized CHC variation using CDF analysis and mapped the first five CHC canonical variates (CVs) onto an independently derived period (per) gene + chromosome inversion + mtDNA COI gene for each sex. We found that the COI sequences were not phylogenetically informative due to introgression between some species, so only per + inversion data were used. Positive phylogenetic signal was observed mainly for CV1 when parsimony methods and the test for serial independence (TFSI) were used. These results changed when no outgroup species were included in the analysis and phylogenetic signal was then observed for female CV3 and/or CV4 and male CV4 and CV5. Finally, removal of divergent populations of D. serido significantly increased the amount of phylogenetic signal as up to four out of five CVs then displayed positive phylogenetic signal.

Conclusions

CHCs were conserved among species while quantitative differences in CHC profiles between populations and species were statistically significant. Most CHCs were species-, population-, and sex-specific. Mapping CHCs onto an independently derived phylogeny revealed that a significant portion of CHC variation was explained by species' systematic affinities indicating phylogenetic conservatism in the evolution of these hydrocarbon arrays, presumptive waterproofing compounds and courtship signals as in many other drosophilid species.

Low satellite DNA variability in natural populations of Drosophila antonietae involved in different evolutionary events

Drosophila antonietae is a cactophilic species that is found in the mesophilic forest of the Paraná-Paraguay river basin and in the dunes of the South Atlantic coast of Brazil. Although the genetic structure of the Paraná-Paraguay river basin populations has already been established, the relationship between these populations and those on the Atlantic coast is controversial. In this study, we compared 33 repetitive units of pBuM-2 satellite DNA isolated from individuals from 8 populations of D. antonietae in these geographic regions, including some populations found within a contact zone with the closely related D. serido. The pBuM-2 sequences showed low interpopulational variability. This result was interpreted as a consequence of both gene flow among the populations and unequal crossing over promoting homogenization of the tandem arrays. The results presented here, together with those of previous studies, highlight the use of pBuM-2 for solving taxonomic conflicts within the D. buzzatii species cluster.

The non-regular orbit: three satellite DNAs in Drosophila martensis (buzzatii complex, repleta group) followed three different evolutionary pathways

The genome of species from the buzzatii cluster (buzzatii complex, repleta group) is hosted by a number of satellite DNAs (satDNAs) showing contrasting structural characteristics, genomic organization and evolution, such as pBuM-alpha (~190 bp repeats), pBuM-alpha/beta (~370 bp repeats) and the DBC-150 (~150 bp repeats). In the present study, we aimed to investigate the evolution of these three satDNAs by looking for homologous sequences in the genome of the closest outgroup species: Drosophila martensis (buzzatii complex). After PCR, we isolated and sequenced 9 alpha, 8 alpha/beta and 11 DBC-150 sequences from this species. The results were compared to all pBuM and DBC-150 sequences available in literature. After D. martensis split from the buzzatii cluster some 6 Mya, the three satDNAs evolved differently in the genome of D. martensis by: (1) maintenance of a collection of major types of ancestral repeats in the genome (alpha); (2) fixation for a single major type of ancestral repeats (alpha/beta) or (3) fixation for new divergent species-specific repeat types (DBC-150). Curiously, D. seriema and D. martensis, although belonging to different and allopatric clusters, became independently fixed for the same major type of alpha/beta ancestral repeats, illustrating a rare case of parallelism in satDNA evolution. The contrasting pictures illustrate the diversity of evolutionary pathways a satDNA can follow, defining a "non-regular orbit" with outcomes difficult to predict.

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.