Evolutionary relationships of Drosophila mojavensis geographic host races and their sister species Drosophila arizonae

Elevated sleep quota in a stress-resilient Drosophila species

Sleep is broadly conserved across the animal kingdom but can vary widely between species. It is currently unclear which selective pressures and regulatory mechanisms influence differences in sleep between species. The fruit fly Drosophila melanogaster has become a successful model system for examining sleep regulation and function, but little is known about the sleep patterns in many related fly species. Here, we find that fly species with adaptations to extreme desert environments, including D. mojavensis, exhibit strong increases in baseline sleep compared with D. melanogaster. Long-sleeping D. mojavensis show intact homeostasis, indicating that desert flies carry an elevated drive for sleep. In addition, D. mojavensis exhibit altered abundance or distribution of several sleep/wake-related neuromodulators and neuropeptides that are consistent with their reduced locomotor activity and increased sleep. Finally, we find that in a nutrient-deprived environment, the sleep patterns of individual D. mojavensis are strongly correlated with their survival time and that disrupting sleep via constant light stimulation renders D. mojavensis more sensitive to starvation. Our results demonstrate that D. mojavensis is a novel model for studying organisms with high sleep drive and for exploring sleep strategies that provide resilience in extreme environments.

Evolution of a fatty acyl-CoA elongase underlies desert adaptation in Drosophila

Traits that allow species to survive in extreme environments such as hot-arid deserts have independently evolved in multiple taxa. However, the genetic and evolutionary mechanisms underlying these traits have thus far not been elucidated. Here, we show that Drosophila mojavensis, a desert-adapted fruit fly species, has evolved high desiccation resistance by producing long-chain methyl-branched cuticular hydrocarbons (mbCHCs) that contribute to a cuticular lipid layer reducing water loss. We show that the ability to synthesize these longer mbCHCs is due to evolutionary changes in a fatty acyl-CoA elongase (mElo). mElo knockout in D. mojavensis led to loss of longer mbCHCs and reduction of desiccation resistance at high temperatures but did not affect mortality at either high temperatures or desiccating conditions individually. Phylogenetic analysis showed that mElo is a Drosophila-specific gene, suggesting that while the physiological mechanisms underlying desert adaptation may be similar between species, the genes involved in these mechanisms may be species or lineage specific.

Elevated sleep need in a stress-resilient Drosophila species

Sleep is broadly conserved across the animal kingdom, but can vary widely between species. It is currently unclear which types of selective pressures and sleep regulatory mechanisms influence differences in sleep between species. The fruit fly Drosophila melanogaster has become a successful model system for examining sleep regulation and function, but little is known about the sleep patterns and need for sleep in many related fly species. Here, we find that Drosophila mojavensis, a fly species that has adapted to extreme desert environments, exhibits strong increases in sleep compared to D. melanogaster. Long-sleeping D. mojavensis show intact sleep homeostasis, indicating that these flies carry an elevated need for sleep. In addition, D. mojavensis exhibit altered abundance or distribution of several sleep/wake related neuromodulators and neuropeptides that are consistent with their reduced locomotor activity, and increased sleep. Finally, we find that in a nutrient-deprived environment, the sleep responses of individual D. mojavensis are correlated with their survival time. Our results demonstrate that D. mojavensis is a novel model for studying organisms with high sleep need, and for exploring sleep strategies that provide resilience in extreme environments.

Physiological and life history changes associated with seasonal adaptation in the cactophilic Drosophila mojavensis

Many insects inhabiting temperate climates are faced with changing environmental conditions throughout the year. Depending on the species, these environmental fluctuations can be experienced within a single generation or across multiple generations. Strategies for dealing with these seasonal changes vary across populations. Drosophila mojavensis is a cactophilic Drosophila species endemic to the Sonoran Desert. The Sonoran Desert regularly reaches temperatures of 50°C in the summer months. As individuals of this population are rare to collect in the summer months, we simulated the cycling temperatures experienced by D. mojavensis in the Sonoran Desert from April to July (four generations) in a temperature- and light-controlled chamber, to understand the physiological and life history changes that allow this population to withstand these conditions. In contrast to our hypothesis of a summer aestivation, we found that D. mojavensis continue to reproduce during the summer months, albeit with lower viability, but the adult survivorship of the population is highly reduced during this period. As expected, stress resistance increased during the summer months in both the adult and the larval stages. This study examines several strategies for withstanding the Sonoran Desert summer conditions which may be informative in the study of other desert endemic species.

Unveiling the Mycodrosophila projectans (Diptera, Drosophilidae) species complex: Insights into the evolution of three Neotropical cryptic and syntopic species

The Zygothrica genus group has been shown to be speciose, with a high number of cryptic species. DNA barcoding approaches have been a valuable tool to uncover cryptic diversity in this lineage, as recently suggested for the Neotropical Mycodrosophila projectans complex, which seems to comprise at least three different species. The aim of this study was to confirm the subdivision of the M. projectans complex while shedding some light on the patterns and processes related to its diversification. In this sense, the use of single and multi-locus datasets under phylogenetic, distance, coalescence, and diagnostic nucleotide approaches confirmed the presence of at least three species under the general morphotype previously described as M. projectans. Only a few subtle morphological differences were found for the three species in terms of aedeagus morphology and abdominal color patterns. Ecologically, sympatry and syntopy seem to be recurrent for these three cryptic species, which present widely overlapping niches, implying niche conservatism. This morphological and ecological similarity has persisted though cladogenesis within the complex, which dates back to the Miocene, providing an interesting example of morphological conservation despite ancient divergence. These results, in addition to contrasting patterns of past demographic fluctuations, allowed us to hypothesize patterns of allopatric or parapatric diversification with secondary contact in Southern Brazil. Nevertheless, genetic diversity was generally high within species, suggesting that migration may encompass an adaptive response to the restrictions imposed by the ephemerality of resources.

Transposable Element Expression and Regulation Profile in Gonads of Interspecific Hybrids of Drosophila arizonae and Drosophila mojavensis wrigleyi

Interspecific hybridization may lead to sterility and/or inviability through differential expression of genes and transposable elements (TEs). In Drosophila, studies have reported massive TE mobilization in hybrids from interspecific crosses of species presenting high divergence times. However, few studies have examined the consequences of TE mobilization upon hybridization in recently diverged species, such as Drosophila arizonae and D. mojavensis. We have sequenced transcriptomes of D. arizonae and the subspecies D. m. wrigleyi and their reciprocal hybrids, as well as piRNAs, to analyze the impact of genomic stress on TE regulation. Our results revealed that the differential expression in both gonadal tissues of parental species was similar. Globally, ovaries and testes showed few deregulated TEs compared with both parental lines. Analyses of small RNA data showed that in ovaries, the TE upregulation is likely due to divergence of copies inherited from parental genomes and lack of piRNAs mapping to them. Nevertheless, in testes, the divergent expression of genes associated with chromatin state and piRNA pathway potentially indicates that TE differential expression is related to the divergence of regulatory genes that play a role in modulating transcriptional and post-transcriptional mechanisms.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F₁ hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Comparative transcriptomics between Drosophila mojavensis and D. arizonae reveals transgressive gene expression and underexpression of spermatogenesis-related genes in hybrid testes

Interspecific hybridization is a stressful condition that can lead to sterility and/or inviability through improper gene regulation in Drosophila species with a high divergence time. However, the extent of these abnormalities in hybrids of recently diverging species is not well known. Some studies have shown that in Drosophila, the mechanisms of postzygotic isolation may evolve more rapidly in males than in females and that the degree of viability and sterility is associated with the genetic distance between species. Here, we used transcriptomic comparisons between two Drosophila mojavensis subspecies and D. arizonae (repleta group, Drosophila) and identified greater differential gene expression in testes than in ovaries. We tested the hypothesis that the severity of the interspecies hybrid phenotype is associated with the degree of gene misregulation. We showed limited gene misregulation in fertile females and an increase in the amount of misregulation in males with more severe sterile phenotypes (motile vs. amotile sperm). In addition, for these hybrids, we identified candidate genes that were mostly associated with spermatogenesis dysfunction.

Host Preference and Olfaction in Drosophila mojavensis

Many organisms live in complex environments that vary geographically in resource availability. This environmental heterogeneity can lead to changes within species in their phenotypic traits. For example, in many herbivorous insects, variation in host plant availability has been shown to influence insect host preference behavior. This behavior can be mediated in part through the insect olfactory system and the odor-evoked responses of olfactory sensory neurons (OSNs), which are in turn mediated by their corresponding odorant receptor genes. The desert dwelling fly Drosophila mojavensis is a model species for understanding the mechanisms underlying host preference in a heterogeneous environment. Depending on geographic region, one to multiple host plant species are available. Here, we conducted electrophysiological studies and found variation in responses of ORNs to host plant volatiles both within and between 2 populations-particularly to the odorant 4-methylphenol. Flies from select localities within each population were found to lack a response to 4-methylphenol. Experiments then assessed the extent to which these electrophysiological differences were associated with differences in several odor-mediated behavioral responses. No association between the presence/absence of these odor-evoked responses and short range olfactory behavior or oviposition behavior was observed. However, differences in odor-induced feeding behavior in response to 4-methylphenol were found. Localities that exhibit an odor-evoked response to the odorant had increased feeding behavior in the presence of the odorant. This study sets the stage for future work examining the functional genetics underlying variation in odor perception.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

Genetic Diversity and Demographic History in the Cactophilic Drosophila repleta Species Group (Diptera: Drosophilidae) in North America Inferred from Mitochondrial DNA Barcodes

Genetic diversity in mitochondrial DNA barcodes, comprising a segment of the cytochrome c oxidase subunit I (COI) gene, was used to infer demographic histories in selected taxa of the cactophilic Drosophila repleta species group in North America. Haplotype and nucleotide diversities were determined in 16 taxa based on both previously published and new sequences. Haplotype diversity (h) differed dramatically in different taxa, varying from h = 0 in Drosophila eremophila, Drosophila hexastigma, and Drosophila bifurca to h = 0.99 in Drosophila hamatofila. Genetic diversity indices and sample sizes were sufficient to infer demographic histories from mismatch distribution analysis and Bayesian skyline plots for 9 taxa: Drosophila mojavensis baja, Drosophila mojavensis sonorensis, Drosophila arizonae, Drosophila aldrichi, D. hamatofila, Drosophila spenceri, Drosophila mainlandi, Drosophila mettleri, and Drosophila nigrospiracula. Evidence was found for both population expansions and relatively stable populations in these species. Demographic history varied dramatically in subspecies of D. mojavensis, showing a relatively stable population size over time in D. m. sonorensis from the mainland Sonoran Desert whereas a large population expansion was evident in D. m. baja from the Baja California Peninsula, providing support for the hypothesis that the split of sister species D. mojavensis and D. arizonae from a common ancestor occurred on the mainland rather than the peninsula as proposed by others. No evidence was found for a causal relationship between a stable or expanding population and host plant shifts from prickly-pear cactus to columnar cacti, which has occurred independently in many taxa of the repleta group.

Genomic analysis of the four ecologically distinct cactus host populations of Drosophila mojavensis

BACKGROUND

Relationships between an organism and its environment can be fundamental in the understanding how populations change over time and species arise. Local ecological conditions can shape variation at multiple levels, among these are the evolutionary history and trajectories of coding genes. This study examines the rate of molecular evolution at protein-coding genes throughout the genome in response to host adaptation in the cactophilic Drosophila mojavensis. These insects are intimately associated with cactus necroses, developing as larvae and feeding as adults in these necrotic tissues. Drosophila mojavensis is composed of four isolated populations across the deserts of western North America and each population has adapted to utilize different cacti that are chemically, nutritionally, and structurally distinct.

RESULTS

High coverage Illumina sequencing was performed on three previously unsequenced populations of D. mojavensis. Genomes were assembled using the previously sequenced genome of D. mojavensis from Santa Catalina Island (USA) as a template. Protein coding genes were aligned across all four populations and rates of protein evolution were determined for all loci using a several approaches.

CONCLUSIONS

Loci that exhibited elevated rates of molecular evolution tend to be shorter, have fewer exons, low expression, be transcriptionally responsive to cactus host use and have fixed expression differences across the four cactus host populations. Fast evolving genes were involved with metabolism, detoxification, chemosensory reception, reproduction and behavior. Results of this study give insight into the process and the genomic consequences of local ecological adaptation.

Assessing the Architecture of Drosophila mojavensis Locomotor Evolution with Bulk Segregant Analysis

Behavior is frequently predicted to be especially important for evolution in novel environments. If these predictions are accurate, there might be particular patterns of genetic architecture associated with recently diverged behaviors. Specifically, it has been predicted that behaviors linked to population divergence should be underpinned by a few genes of relatively large effect, compared to architectures of intrapopulation behavioral variation, which is considered to be highly polygenic. More mapping studies of behavioral variation between recently diverged populations are needed to continue assessing the generality of these predictions. Here, we used a bulk segregant mapping approach to dissect the genetic architecture of a locomotor trait that has evolved between two populations of the cactophilic fly Drosophila mojavensis We created an F8 mapping population of 1,500 individuals from advanced intercross lines and sequenced the 10% of individuals with the highest and lowest levels of locomotor activity. Using three alternative statistical approaches, we found strong evidence for two relatively large-effect QTL that is localized in a region homologous to a region of densely packed behavior loci in Drosophila melanogaster, suggesting that clustering of behavior genes may display relatively deep evolutionary conservation. Broadly, our data are most consistent with a polygenic architecture, though with several loci explaining a high proportion of variation in comparison to similar behavioral traits. We further note the presence of several antagonistic QTL linked to locomotion and discuss these results in light of theories regarding behavioral evolution and the effect size and direction of QTL for diverging traits in general.

Inverse resource allocation between vision and olfaction across the genus Drosophila

Divergent populations across different environments are exposed to critical sensory information related to locating a host or mate, as well as avoiding predators and pathogens. These sensory signals generate evolutionary changes in neuroanatomy and behavior; however, few studies have investigated patterns of neural architecture that occur between sensory systems, or that occur within large groups of closely-related organisms. Here we examine 62 species within the genus Drosophila and describe an inverse resource allocation between vision and olfaction, which we consistently observe at the periphery, within the brain, as well as during larval development. This sensory variation was noted across the entire genus and appears to represent repeated, independent evolutionary events, where one sensory modality is consistently selected for at the expense of the other. Moreover, we provide evidence of a developmental genetic constraint through the sharing of a single larval structure, the eye-antennal imaginal disc. In addition, we examine the ecological implications of visual or olfactory bias, including the potential impact on host-navigation and courtship.

Genomic changes associated with adaptation to arid environments in cactophilic Drosophila species

BACKGROUND

Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. Here we investigate the genetic changes associated with adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species.

RESULTS

We find disproportionately high rates of gene gains in internal branches in the species' phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species' differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one.

CONCLUSION

Adaptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change.

An Improved Genome Assembly for Drosophila navojoa, the Basal Species in the mojavensis Cluster

Three North American cactophilic Drosophila species, D. mojavensis, D. arizonae, and D. navojoa, are of considerable evolutionary interest owing to the shift from breeding in Opuntia cacti to columnar species. The 3 species form the "mojavensis cluster" of Drosophila. The genome of D. mojavensis was sequenced in 2007 and the genomes of D. navojoa and D. arizonae were sequenced together in 2016 using the same technology (Illumina) and assembly software (AllPaths-LG). Yet, unfortunately, the D. navojoa genome was considerably more fragmented and incomplete than its sister species, rendering it less useful for evolutionary genetic studies. The D. navojoa read dataset does not fully meet the strict insert size required by the assembler used (AllPaths-LG) and this incompatibility might explain its assembly problems. Accordingly, when we re-assembled the genome of D. navojoa with the SPAdes assembler, which does not have the strict AllPaths-LG requirements, we obtained a substantial improvement in all quality indicators such as N50 (from 84 kb to 389 kb) and BUSCO coverage (from 77% to 97%). Here we share a new, improved reference assembly for D. navojoa genome, along with a RNAseq transcriptome. Given the basal relationship of the Opuntia breeding D. navojoa to the columnar breeding D. arizonae and D. mojavensis, the improved assembly and annotation will allow researchers to address a range of questions associated with the genomics of host shifts, chromosomal rearrangements and speciation in this group.

Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

Background

Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. Here we aim to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae.

Results

We used a phylogenetic and population genetic analysis framework to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. We found strong evidence of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains.

Conclusions

The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1250-x) contains supplementary material, which is available to authorized users.

Behavioral evolution accompanying host shifts in cactophilic Drosophila larvae

For plant utilizing insects, the shift to a novel host is generally accompanied by a complex set of phenotypic adaptations. Many such adaptations arise in response to differences in plant chemistry, competitive environment, or abiotic conditions. One less well-understood factor in the evolution of phytophagous insects is the selective environment provided by plant shape and volume. Does the physical structure of a new plant host favor certain phenotypes? Here, we use cactophilic Drosophila, which have colonized the necrotic tissues of cacti with dramatically different shapes and volumes, to examine this question. Specifically, we analyzed two behavioral traits in larvae, pupation height, and activity that we predicted might be related to the ability to utilize variably shaped hosts. We found that populations of D. mojavensis living on lengthy columnar or barrel cactus hosts have greater activity and pupate higher in a laboratory environment than populations living on small and flat prickly pear cactus cladodes. Crosses between the most phenotypically extreme populations suggest that the genetic architectures of these behaviors are distinct. A comparison of activity in additional cactophilic species that are specialized on small and large cactus hosts shows a consistent trend. Thus, we suggest that greater motility and an associated tendency to pupate higher in the laboratory are potential larval adaptations for life on a large plant where space is more abundant and resources may be more sparsely distributed.

Mutation predicts 40 million years of fly wing evolution

Mutation enables evolution, but the idea that adaptation is also shaped by mutational variation is controversial. Simple evolutionary hypotheses predict such a relationship if the supply of mutations constrains evolution, but it is not clear that constraints exist, and, even if they do, they may be overcome by long-term natural selection. Quantification of the relationship between mutation and phenotypic divergence among species will help to resolve these issues. Here we use precise data on over 50,000 Drosophilid fly wings to demonstrate unexpectedly strong positive relationships between variation produced by mutation, standing genetic variation, and the rate of evolution over the last 40 million years. Our results are inconsistent with simple constraint hypotheses because the rate of evolution is very low relative to what both mutational and standing variation could allow. In principle, the constraint hypothesis could be rescued if the vast majority of mutations are so deleterious that they cannot contribute to evolution, but this also requires the implausible assumption that deleterious mutations have the same pattern of effects as potentially advantageous ones. Our evidence for a strong relationship between mutation and divergence in a slowly evolving structure challenges the existing models of mutation in evolution.

Early events in speciation: Cryptic species of Drosophila aldrichi

Understanding the earliest events in speciation remains a major challenge in evolutionary biology. Thus identifying species whose populations are beginning to diverge can provide useful systems to study the process of speciation. Drosophila aldrichi, a cactophilic fruit fly species with a broad distribution in North America, has long been assumed to be a single species owing to its morphological uniformity. While previous reports either of genetic divergence or reproductive isolation among different D. aldrichi strains have hinted at the existence of cryptic species, the evolutionary relationships of this species across its range have not been thoroughly investigated. Here we show that D. aldrichi actually is paraphyletic with respect to its closest relative, Drosophila wheeleri, and that divergent D. aldrichi lineages show complete hybrid male sterility when crossed. Our data support the interpretation that there are at least two species of D. aldrichi, making these flies particularly attractive for studies of speciation in an ecological and geographical context.

Identification of misexpressed genetic elements in hybrids between Drosophila-related species

Crosses between close species can lead to genomic disorders, often considered to be the cause of hybrid incompatibility, one of the initial steps in the speciation process. How these incompatibilities are established and what are their causes remain unclear. To understand the initiation of hybrid incompatibility, we performed reciprocal crosses between two species of Drosophila (D. mojavensis and D. arizonae) that diverged less than 1 Mya. We performed a genome-wide transcriptomic analysis on ovaries from parental lines and on hybrids from reciprocal crosses. Using an innovative procedure of co-assembling transcriptomes, we show that parental lines differ in the expression of their genes and transposable elements. Reciprocal hybrids presented specific gene categories and few transposable element families misexpressed relative to the parental lines. Because TEs are mainly silenced by piwi-interacting RNAs (piRNAs), we hypothesize that in hybrids the deregulation of specific TE families is due to the absence of such small RNAs. Small RNA sequencing confirmed our hypothesis and we therefore propose that TEs can indeed be major players of genome differentiation and be implicated in the first steps of genomic incompatibilities through small RNA regulation.

Lipidomic profiles of Drosophila melanogaster and cactophilic fly species: models of human metabolic diseases

Untargeted metabolomics, combined with data mining, reveals different sensibility of fly species against diet changes.

Genome Evolution in Three Species of Cactophilic Drosophila

We report genomes of two species of cactophilic Drosophila: Drosophila arizonae and D. navojoa These two are the closest relatives of D. mojavensis, forming the D. mojavensis cluster. D. mojavensis and D. arizonae diverged from D. navojoa ∼5.8 Mya, while the split between D. arizonae and D. mojavensis is more recent, at 1.5 Mya. Together the three genomes provide opportunities to examine genomic changes associated with speciation and host shifts in this ecologically defined group of flies. The three species are also separated by fixed inversion differences in three of their six chromosomes. While the levels of nucleotide divergence in the colinear chromosomes are significantly lower than in the inverted chromosomes, consistent with a past role of the inversions in preventing gene flow, the patterns differ among the inverted chromosomes when the locations of nucleotides inside or outside of the inversions are considered. For Muller element E, there is greater divergence external to the inversion breakpoints. For Muller A, the divergence is slightly higher inside the inversions, while for Muller B, the breakpoints and hence the difference in substitutions in relation to the inversions could not be determined. The differences among the inverted chromosomes, especially once the breakpoints are clearly established, could aid in dating the origins of the inversions.

The changing biodiversity of Alabama Drosophila: important impacts of seasonal variation, urbanization, and invasive species

Global warming and anthropogenic disturbances significantly influence the biosphere, tremendously increasing species extinction rates. In Central Alabama, we analyzed Drosophilidae species composition change nearly 100 years after the previous survey. We found ten Drosophilid species that were not reported during the last major biodiversity studies, two of which are invasive pests. In addition, we analyzed the influence of seasonal environmental variables characteristic of the subtropical climate zone on Drosophila abundance and biodiversity. We found a significant correlation between temperature and abundance of total Drosophila as well as for six of the seven most represented species individually, with a maximum abundance at intermediate temperatures (18–26°C). In addition, temperature was positively correlated with biodiversity of Drosophila. Precipitation produced a significant effect on the abundance of five species of Drosophila, with different optima for each species, but did not affect overall biodiversity. We concluded that in the subtropical climate zone of Central Alabama, seasonal temperature and precipitation changes produce a significant effect on Drosophila abundance and biodiversity, while local land use also impacts fly abundance, contributing to an apparent shift in species composition over the last century. We expect global climate change and other anthropogenic factors to further impact Drosophila species composition in the subtropical climate zone into the future.

SNP calling from RNA-seq data without a reference genome: identification, quantification, differential analysis and impact on the protein sequence

SNPs (Single Nucleotide Polymorphisms) are genetic markers whose precise identification is a prerequisite for association studies. Methods to identify them are currently well developed for model species, but rely on the availability of a (good) reference genome, and therefore cannot be applied to non-model species. They are also mostly tailored for whole genome (re-)sequencing experiments, whereas in many cases, transcriptome sequencing can be used as a cheaper alternative which already enables to identify SNPs located in transcribed regions. In this paper, we propose a method that identifies, quantifies and annotates SNPs without any reference genome, using RNA-seq data only. Individuals can be pooled prior to sequencing, if not enough material is available from one individual. Using pooled human RNA-seq data, we clarify the precision and recall of our method and discuss them with respect to other methods which use a reference genome or an assembled transcriptome. We then validate experimentally the predictions of our method using RNA-seq data from two non-model species. The method can be used for any species to annotate SNPs and predict their impact on the protein sequence. We further enable to test for the association of the identified SNPs with a phenotype of interest.

Molecular evolution of candidate genes involved in post-mating-prezygotic reproductive isolation

Traits involved in post-copulatory interactions between the sexes may evolve rapidly as a result of sexual selection and/or sexual conflict, leading to post-mating-prezygotic (PMPZ) reproductive isolating barriers between diverging lineages. Although the importance of PMPZ isolation is recognized, the molecular basis of such incompatibilities is not well understood. Here, we investigate molecular evolution of a subset of Drosophila mojavensis and Drosophila arizonae reproductive tract genes. These include genes that are transcriptionally regulated by conspecific mating in females, many of which are misregulated in heterospecific crosses, and a set of male genes whose transcripts are transferred to females during mating. As a group, misregulated female genes are not more divergent and do not appear to evolve under different selection pressures than other female reproductive genes. Male transferred genes evolve at a higher rate than testis-expressed genes, and at a similar rate compared to accessory gland protein genes, which are known to evolve rapidly. Four of the individual male transferred genes show patterns of divergent positive selection between D. mojavensis and D. arizonae. Three of the four genes belong to the sperm-coating protein-like family, including an ortholog of antares, which influences female fertility and receptivity in Drosophila melanogaster. Synthesis of these molecular evolutionary analyses with transcriptomics and predicted functional information makes these genes candidates for involvement in PMPZ reproductive incompatibilities between D. mojavensis and D. arizonae.

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.

The role of aedeagus size and shape in failed mating interactions among recently diverged taxa in the Drosophila mojavensis species cluster

Background

Investigating the evolution of species-specific insect genitalia is central to understanding how morphological diversification contributes to reproductive isolation and lineage divergence. While many studies evoke some form of sexual selection to explain genitalia diversity, the basis of selection and the mechanism of heterospecific mate exclusion remains vague. I conducted reciprocal mate pair trials in the Drosophila mojavensis species cluster to quantify the frequency of failed insemination attempts, historically referred to as pseudocopulation, between lineages with discrete size and shape differences of the male aedeagus.

Results

In cross-taxon matings aedeagus size had a significant effect on pseudocopulation frequencies, while aedeagus shape and genetic distance did not. The direction of the size difference was an important factor for successful mating. When females were mated to a cross-taxon male with a larger aedeagus than males from her own species, the pair could not establish a successful mating interaction. Females mated to cross-taxon males with a smaller aedeagus than conspecific males were able to establish the mating interaction but had issues disengaging at the end of the interaction.

Conclusions

The results of this study support a role for aedeagus size in the male-female mating interaction, with a secondary role for aedeagus shape. In natural populations, mating failure based on aedeagus size could serve as an important reproductive isolating mechanism resulting in failed insemination attempts after both the male and female show a willingness to mate.

Electronic supplementary material

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

Sex-specific triacylglycerides are widely conserved in Drosophila and mediate mating behavior

Pheromones play an important role in the behavior, ecology, and evolution of many organisms. The structure of many insect pheromones typically consists of a hydrocarbon backbone, occasionally modified with various functional oxygen groups. Here we show that sex-specific triacylclyerides (TAGs) are broadly conserved across the subgenus Drosophila in 11 species and represent a novel class of pheromones that has been largely overlooked. In desert-adapted drosophilids, 13 different TAGs are secreted exclusively by males from the ejaculatory bulb, transferred to females during mating, and function synergistically to inhibit courtship from other males. Sex-specific TAGs are comprised of at least one short branched tiglic acid and a long linear fatty acyl component, an unusual structural motif that has not been reported before in other natural products. The diversification of chemical cues used by desert-adapted Drosophila as pheromones may be related to their specialized diet of fermenting cacti.

DOI:http://dx.doi.org/10.7554/eLife.01751.001

For animals, the ultimate purpose of life is to have sex, as nothing is more important than passing down your genes to future generations. A wide range of strategies are therefore employed throughout nature to maximize the chances of sexual success, from ostentatious courtship rituals to the subtle subliminal signals sent out using chemicals called pheromones. Plants and animals release pheromones to influence the behavior of other plants and animals, often without the recipient being aware of it.

Hundreds of different insect pheromones have been discovered. Fruit flies release a number of different pheromones, all with similar chemical structures. Now, Chin et al. have discovered that male flies belonging to several species of fruit fly that live in the desert release chemicals called triacylglycerides (TAGs), which are commonly used for energy storage by many organisms as pheromones. During sex, the male fly rubs the TAGs onto the body of the female, which makes her less attractive to other male flies for several hours, thus increasing his chances of parenthood and passing his genes to future generations.

TAGs are also found in other insect species, but have been largely overlooked as pheromones. Moreover, the TAGs discovered by Chin et al. have an unusual structure, not previously seen in nature, which may result from the diet of fermenting cacti the desert-dwelling fruit flies enjoy.

DOI:http://dx.doi.org/10.7554/eLife.01751.002

Genome-wide dissection of hybrid sterility in Drosophila confirms a polygenic threshold architecture

To date, different studies about the genetic basis of hybrid male sterility (HMS), a postzygotic reproductive barrier thoroughly investigated using Drosophila species, have demonstrated that no single major gene can produce hybrid sterility without the cooperation of several genetic factors. Early work using hybrids between Drosophila koepferae (Dk) and Drosophila buzzatii (Db) was consistent with the idea that HMS requires the cooperation of several genetic factors, supporting a polygenic threshold (PT) model. Here we present a genome-wide mapping strategy to test the PT model, analyzing serially backcrossed fertile and sterile males in which the Dk genome was introgressed into the Db background. We identified 32 Dk-specific markers significantly associated with hybrid sterility. Our results demonstrate 1) a strong correlation between the number of segregated sterility markers and males' degree of sterility, 2) the exchangeability among markers, 3) their tendency to cluster into low-recombining chromosomal regions, and 4) the requirement for a minimum number (threshold) of markers to elicit sterility. Although our findings do not contradict a role for occasional major hybrid-sterility genes, they conform more to the view that HMS primarily evolves by the cumulative action of many interacting genes of minor effect in a complex PT architecture.

Ecological genomics of host shifts in Drosophila mojavensis

Advances in next-generation sequencing technologies have liberated our dependency on model laboratory species for answering genomic and transcriptomic level questions. These new techniques have dramatically expanded our breadth of study organisms and have allowed the analysis of species from diverse ecological environments. One such species is the cactophilic Drosophila mojavensis that inhabits the deserts of western North America. These insects feed and develop in the necrotic cacti, feeding largely on the microflora of the necrotic plant tissues. Drosophila mojavensis is composed of four geographically and ecologically separated populations. Each population (Baja California peninsula, mainland Sonoran Desert, Mojave Desert and Santa Catalina Island) utilizes the necrotic tissues of distinct cactus species. The differences in the nutritional and chemical composition of the necroses include a set of toxic compounds to which resident population must adapt. These ecological differences have facilitated many of the life history, behavior, physiological and genetic differences between the cactus host populations. Genomic resources have allowed investigators to examine the genomic and transcriptional level changes associated with the local adaptation of the four D. mojavensis populations, thereby providing further understanding of the genetic mechanism of adaptation and its role in the divergence of ecologically distinct populations.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.

Genetic diversity and population genetics of mosquitoes (Diptera: Culicidae: Culex spp.) from the Sonoran Desert of North America

The population genetics and phylogenetic relationships of Culex mosquitoes inhabiting the Sonoran Desert region of North America were studied using mitochondrial DNA and microsatellite molecular markers. Phylogenetic analyses of mitochondrial cytochrome c oxidase subunit I (COI) from mosquitoes collected over a wide geographic area, including the Baja California peninsula, and mainland localities in southern Arizona, USA and Sonora, Mexico, showed several well-supported partitions corresponding to Cx. quinquefasciatus, Cx. tarsalis, and two unidentified species, Culex sp. 1 and sp. 2. Culex quinquefasciatus was found at all localities and was the most abundant species collected. Culex tarsalis was collected only at Tucson, Arizona and Guaymas, Sonora. The two unidentified species of Culex were most abundant at Navojoa in southern Sonora. Haplotype and nucleotide diversities in the COI gene segment were substantially lower in Cx. quinquefasciatus compared with the other three species. Analysis of molecular variance revealed little structure among seven populations of Cx. quinquefasciatus, whereas significant structure was found between the two populations of Cx. tarsalis. Evidence for an historical population expansion beginning in the Pleistocene was found for Cx. tarsalis. Possible explanations for the large differences in genetic diversity between Cx. quinquefasciatus and the other species of Culex are presented.

Divergence in olfactory host plant preference in D. mojavensis in response to cactus host use

Divergence in host adaptive traits has been well studied from an ecological and evolutionary perspective, but identification of the proximate mechanisms underlying such divergence is less well understood. Behavioral preferences for host plants are often mediated by olfaction and shifts in preference may be accompanied by changes in the olfactory system. In this study, we examine the evolution of host plant preferences in cactophilic Drosophila mojavensis that feeds and breeds on different cacti throughout its range. We show divergence in electrophysiological responses and olfactory behavior among populations with host plant shifts. Specifically, significant divergence was observed in the Mojave Desert population that specializes on barrel cactus. Differences were observed in electrophysiological responses of the olfactory organs and in behavioral responses to barrel cactus volatiles. Together our results suggest that the peripheral nervous system has changed in response to different ecological environments and that these changes likely contribute to divergence among D. mojavensis populations.

Functional genomic and phenotypic responses to desiccation in natural populations of a desert drosophilid

We used whole-transcriptome microarrays to assess changes in gene expression and monitored mortality rates and epicuticular hydrocarbons (CHCs) in response to desiccation stress in four natural populations of Drosophila mojavensis from Baja California and mainland Mexico. Desiccation had the greatest effect on gene expression, followed by biogeographical variation at regional and population levels. Genes involved in environmental sensing and cuticular structure were up-regulated in dry conditions, while genes involved in transcription itself were down-regulated. Flies from Baja California had higher expression of reproductive and mitochondrial genes, suggesting that these populations have greater fecundity and higher metabolic rates. Host plant differences had a surprisingly minor effect on the transcriptome. In most cases, desiccation-caused mortality was greater in flies reared on fermenting cactus tissues than that on laboratory media. Water content of adult females and males was significantly different and was lower in Baja California males. Different groups of CHCs simultaneously increased and decreased in amounts due to desiccation exposure of 9 and 18 h and were population-specific and dependent on larval rearing substrates. Overall, we observed that changes in gene expression involved a coordinated response of behavioural, cuticular and metabolic genes. Together with differential expression of cuticular hydrocarbons, this study revealed some of the mechanisms that have allowed D. mojavensis to exploit its harsh desert conditions. Certainly, for D. mojavensis that uses different host plants, population-level understanding of responses to stressors associated with future climate change in desert regions must be evaluated across geographical and local ecological scales.

Transcriptome-wide expression variation associated with environmental plasticity and mating success in cactophilic Drosophila mojavensis

Ecological speciation occurs with the adaptation of populations to different environments and concurrent evolution of reproductive isolation. Phenotypic plasticity might influence both ecological adaptation and reproductive traits. We examined environment-specific gene expression and male mating success in cactophilic Drosophila mojavensis using transcriptome sequencing. This species exhibits cactus-dependent mating success across different species of host plants, with genotype-by-environment interactions for numerous traits. We cultured flies from egg to eclosion on two natural cactus hosts and surveyed gene expression in adult males that were either successful or unsuccessful in achieving copulation in courtship trials. We identified gene expression differences that included functions involved with metabolism, most likely related to chemical differences between host cactus species. Several epigenetic-related functions were identified that might play a role in modulating gene expression in adults due to host cactus effects on larvae, and mating success. Cactus-dependent mating success involved expression differences of genes implicated in translation, transcription, and nervous system development. This suggests a role of neurological function genes in the mating success of D. mojavensis males. Together, these results suggest that the influence of environmental variation on mating success via regulation of gene expression might be an important aspect of ecological speciation.

Striking structural dynamism and nucleotide sequence variation of the transposon Galileo in the genome of Drosophila mojavensis

Background

Galileo is a transposable element responsible for the generation of three chromosomal inversions in natural populations of Drosophila buzzatii. Although the most characteristic feature of Galileo is the long internally-repetitive terminal inverted repeats (TIRs), which resemble the Drosophila Foldback element, its transposase-coding sequence has led to its classification as a member of the P-element superfamily (Class II, subclass 1, TIR order). Furthermore, Galileo has a wide distribution in the genus Drosophila, since it has been found in 6 of the 12 Drosophila sequenced genomes. Among these species, D. mojavensis, the one closest to D. buzzatii, presented the highest diversity in sequence and structure of Galileo elements.

Results

In the present work, we carried out a thorough search and annotation of all the Galileo copies present in the D. mojavensis sequenced genome. In our set of 170 Galileo copies we have detected 5 Galileo subfamilies (C, D, E, F, and X) with different structures ranging from nearly complete, to only 2 TIR or solo TIR copies. Finally, we have explored the structural and length variation of the Galileo copies that point out the relatively frequent rearrangements within and between Galileo elements. Different mechanisms responsible for these rearrangements are discussed.

Conclusions

Although Galileo is a transposable element with an ancient history in the D. mojavensis genome, our data indicate a recent transpositional activity. Furthermore, the dynamism in sequence and structure, mainly affecting the TIRs, suggests an active exchange of sequences among the copies. This exchange could lead to new subfamilies of the transposon, which could be crucial for the long-term survival of the element in the genome.

Phylogenetic patterns of geographical and ecological diversification in the subgenus Drosophila

Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. Members of the subgenus Drosophila are distributed across the globe and show a large diversity of ecological niches. Furthermore, taxonomic classification of Drosophila includes the rank radiation, which refers to closely related species groups. Nevertheless, it has never been tested if these taxonomic radiations correspond to evolutionary radiations. Here we present a study of the patterns of diversification of Drosophila to test for increased diversification rates in relation to the geographic and ecological diversification processes. For this, we have estimated and dated a phylogeny of 218 species belonging to the major species groups of the subgenus. The obtained phylogenies are largely consistent with previous studies and indicate that the major groups appeared during the Oligocene/Miocene transition or early Miocene, characterized by a trend of climate warming with brief periods of glaciation. Ancestral reconstruction of geographic ranges and ecological resource use suggest at least two dispersals to the Neotropics from the ancestral Asiatic tropical disribution, and several transitions to specialized ecological resource use (mycophagous and cactophilic). Colonisation of new geographic regions and/or of new ecological resources can result in rapid species diversification into the new ecological niches available. However, diversification analyses show no significant support for adaptive radiations as a result of geographic dispersal or ecological resource shift. Also, cactophily has not resulted in an increase in the diversification rate of the repleta and related groups. It is thus concluded that the taxonomic radiations do not correspond to adaptive radiations.

Model-based comparisons of phylogeographic scenarios resolve the intraspecific divergence of cactophilic Drosophila mojavensis

The cactophilic fly Drosophila mojavensis exhibits considerable intraspecific genetic structure across allopatric geographic regions and shows associations with different host cactus species across its range. The divergence between these populations has been studied for more than 60years, yet their exact historical relationships have not been resolved. We analysed sequence data from 15 intronic X-linked loci across populations from Baja California, mainland Sonora-Arizona and Mojave Desert regions under an isolation-with-migration model to assess multiple scenarios of divergence. We also compared the results with a pre-existing sequence data set of eight autosomal loci. We derived a population tree with Baja California placed at its base and link their isolation to Pleistocene climatic oscillations. Our estimates suggest the Baja California population diverged from an ancestral Mojave Desert/mainland Sonora-Arizona group around 230,000-270,000years ago, while the split between the Mojave Desert and mainland Sonora-Arizona populations occurred one glacial cycle later, 117,000-135,000years ago. Although we found these three populations to be effectively allopatric, model ranking could not rule out the possibility of a low level of gene flow between two of them. Finally, the Mojave Desert population showed a small effective population size, consistent with a historical population bottleneck. We show that model-based inference from multiple loci can provide accurate information on the historical relationships of closely related groups allowing us to set into historical context a classic system of incipient ecological speciation.

Population transcriptomics of cactus host shifts in Drosophila mojavensis

In the presence of environmental change, natural selection can shape the transcriptome. Under a scenario of environmental change, genotypes that are better able to modulate gene expression to maximize fitness will tend to be favoured. Therefore, it is important to examine gene expression at the population level to distinguish random or neutral gene expression variation from the pattern produced by natural selection. This study investigates the natural variation in transcriptional response to a cactus host shift utilizing the mainland Sonora population of Drosophila mojavensis. Drosophila mojavensis is a cactophilic species composed of four cactus host populations endemic to the deserts of North America. Overall, the change in cactus host was associated with a significant reduction in larval viability as well as the differential expression of 21% of the genome (3109 genes). Among the genes identified were a set of genes previously known to be involved in xenobiotic metabolism, as well as genes involved in cellular energy production, oxidoreductase/carbohydrate metabolism, structural components and mRNA binding. Interestingly, of the 3109 genes whose expression was affected by host use, there was a significant overrepresentation of genes that lacked an orthologous call to the D. melanogaster genome, suggesting the possibility of an accelerated rate of evolution in these genes. Of the genes with a significant cactus effect, the majority, 2264 genes, did not exhibit a significant cactus-by-line interaction. This population-level approach facilitated the identification of genes involved in past cactus host shifts.

Evolution of reproductive morphology among recently diverged taxa in the Drosophila mojavensis species cluster

The morphological evolution of sexual traits informs studies of speciation due to the potential role of these characters in reproductive isolation. In the current study, we quantified and compared genitalic variation within the Drosophila mojavensis species cluster to infer the mode of evolution of the male aedeagus. This system is ideal for such studies due to the opportunity to test and compare levels of variation along a divergence continuum at various taxonomic levels within the group. Shape variation was quantified using elliptic Fourier descriptors and compared among the four D. mojavensis host races, and between D. mojavensis and its sister species Drosophila arizonae. Aedeagus shape was diagnostic for D. arizonae, and among three of the four D. mojavensis subspecies. In each of these cases, there was less variation within subspecies than among subspecies, which is consistent with the pattern predicted if genitalia are evolving according to a punctuated change model, and are involved with mate recognition. However, aedeagus shape in Drosophila mojavensis sonorensis was highly variable and broadly overlapping with the other three subspecies, suggesting aedeagus evolution in this subspecies is more complex and subject to additional evolutionary factors. These results are interpreted and discussed in the context of selection on mate recognition systems and the potential for failed copulation.

Phylogeography of the Cactophilic Drosophila and Other Arthropods Associated with Cactus Necroses in the Sonoran Desert

Studies on the population genetics, phylogenetic relationships, systematics and evolution of arthropods that inhabit necrotic tissue of cacti in the Sonoran Desert of North America are reviewed. These studies have focused upon several species of insects (orders Diptera and Coleoptera) and arachnids (order Pseudoscorpiones). For most taxa studied, little genetic structure and high dispersal ability are found in populations inhabiting the mainland and Baja California peninsula regions of the Sonoran Desert, consistent with the availability of the rotting cactus microhabitat which is patchily distributed and ephemeral. There is evidence, however, that the Gulf of California, which bisects the Sonoran Desert, has played a role in limiting gene flow and promoting speciation in several taxa, including histerid beetles, whereas other taxa, especially Drosophila nigrospiracula and D. mettleri, apparently are able to freely cross the Gulf, probably by taking advantage of the Midriff Islands in the northern Gulf as dispersal “stepping stones”. Genetic evidence has also been found for historical population expansions dating to the Pleistocene and late Pliocene in several taxa. Overall, these studies have provided important insights into how arthropods with different life history traits, but generally restricted to a necrotic cactus microhabitat, have evolved in an environmentally harsh and tectonically active region. In addition, they suggest some taxa for further, and more detailed, hypothesis driven studies of speciation.