Population genetics and geographic variation of alcohol dehydrogenase (Adh) paralogs and glucose-6-phosphate dehydrogenase (G6pd) in Drosophila mojavensis

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.

A Bayesian Mutation-Selection Framework for Detecting Site-Specific Adaptive Evolution in Protein-Coding Genes

In recent years, codon substitution models based on the mutation–selection principle have been extended for the purpose of detecting signatures of adaptive evolution in protein-coding genes. However, the approaches used to date have either focused on detecting global signals of adaptive regimes—across the entire gene—or on contexts where experimentally derived, site-specific amino acid fitness profiles are available. Here, we present a Bayesian site-heterogeneous mutation–selection framework for site-specific detection of adaptive substitution regimes given a protein-coding DNA alignment. We offer implementations, briefly present simulation results, and apply the approach on a few real data sets. Our analyses suggest that the new approach shows greater sensitivity than traditional methods. However, more study is required to assess the impact of potential model violations on the method, and gain a greater empirical sense its behavior on a broader range of real data sets. We propose an outline of such a research program.

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.

Computational Sequence Analysis of Inversion Breakpoint Regions in the Cactophilic Drosophila mojavensis Lineage

We investigated rates of chromosomal evolution in Drosophila mojavensis using whole-genome sequence information from D. mojavensis, Drosophila buzzatii, and Drosophila virilis. Drosophila mojavensis is a cactophilic species of the repleta group living under extreme ecological conditions in the deserts of the Southwestern United States and Northwestern México. The genome of D. buzzatii, another member of the repleta group, was recently sequenced and the largest scaffolds anchored to all chromosomes using diverse procedures. Chromosome organization between D. mojavensis and D. buzzatii was compared using MUMmer and GRIMM software. Our results corroborate previous cytological analyses that indicated chromosome 2 differed between these 2 species by 10 inversions, chromosomes X and 5 differed by one inversion each, and chromosome 4 was homosequential. In contrast, we found that chromosome 3 differed by 5 inversions instead of the expected 2 that were previously inferred by cytological analyses. Thirteen of these inversions occurred in the D. mojavensis lineage: 12 are fixed and one of them is a polymorphic inversion previously described in populations from Sonora and Baja California, México. We previously investigated the breakpoints of chromosome 2 inversions fixed in D. mojavensis. Here we characterized the breakpoint regions of the 5 inversions found in chromosome 3 in order to infer the molecular mechanism that generated each inversion and its putative functional consequences. Overall, our results reveal a number of gene alterations at the inversion breakpoints with putative adaptive consequences that point to natural selection as the cause for fast chromosomal evolution in D. mojavensis.

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.

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.

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.

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.

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.

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.

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.

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.

Gene alterations at Drosophila inversion breakpoints provide prima facie evidence for natural selection as an explanation for rapid chromosomal evolution

Background

Chromosomal inversions have been pervasive during the evolution of the genus Drosophila, but there is significant variation between lineages in the rate of rearrangement fixation. D. mojavensis, an ecological specialist adapted to a cactophilic niche under extreme desert conditions, is a chromosomally derived species with ten fixed inversions, five of them not present in any other species.

Results

In order to explore the causes of the rapid chromosomal evolution in D. mojavensis, we identified and characterized all breakpoints of seven inversions fixed in chromosome 2, the most dynamic one. One of the inversions presents unequivocal evidence for its generation by ectopic recombination between transposon copies and another two harbor inverted duplications of non-repetitive DNA at the two breakpoints and were likely generated by staggered single-strand breaks and repair by non-homologous end joining. Four out of 14 breakpoints lay in the intergenic region between preexisting duplicated genes, suggesting an adaptive advantage of separating previously tightly linked duplicates. Four out of 14 breakpoints are associated with transposed genes, suggesting these breakpoints are fragile regions. Finally two inversions contain novel genes at their breakpoints and another three show alterations of genes at breakpoints with potential adaptive significance.

Conclusions

D. mojavensis chromosomal inversions were generated by multiple mechanisms, an observation that does not provide support for increased mutation rate as explanation for rapid chromosomal evolution. On the other hand, we have found a number of gene alterations at the breakpoints with putative adaptive consequences that directly point to natural selection as the cause of D. mojavensis rapid chromosomal evolution.

Resolving the phylogenetic relationships and evolutionary history of the Drosophila virilis group using multilocus data

The Drosophila virilis group is one of the major lineages of Drosophila previously recognised and it has been used as a model for different types of studies. It comprises 13 species whose phylogenetic relationships are not well resolved. In the present study, six nuclear genes (Adh, fused, Gpdh, NonA, CG9631 and CG7219) and the mitochondrial ribosomal RNA genes (12S-16S) have been used to estimate the evolutionary tree of the group using different methods of phylogenetic reconstruction. Different competing evolutionary hypotheses have also been compared using the Approximately Unbiased test to further evaluate the robustness of the inferred trees. Results are, in general, consistent with previous studies in recovering the four major lineages of the group (D. virilis phylad, Drosophila montana subphylad, Drosophila kanekoi subphylad and Drosophila littoralis subphylad), although D. kanekoi, D. littoralis and Drosophila ezoana are here inferred to be more closely related to the D. virilis phylad than to the D. montana subphylad. The age of the crown group, estimated with a Bayesian method that assumes a relaxed molecular clock, is placed in the late Miocene (∼ 10 Mya). The oldest lineages also appeared during this period (∼ 7.5 to ∼ 8.9 Mya), while the ages of the basal nodes of the montana subphylad and the virilis phylad are located in the early Pliocene (∼ 4.9 and ∼ 4.1 Mya). Major cladogenesis events correlate to geological and palaeoclimatic occurrences that most likely affected the freshwater and deciduous forests where these species are found. The inferred biogeographical history of the group, based on the statistical dispersal-vicariance analysis, indicates that the last common ancestor of the group had a Holarctic distribution from which the North American and the Eurasian lineages evolved as a result of a vicariant event.

Distinguishing among evolutionary models for the maintenance of gene duplicates

Determining the evolutionary forces responsible for the maintenance of gene duplicates is key to understanding the processes leading to evolutionary adaptation and novelty. In his highly prescient book, Susumu Ohno recognized that duplicate genes are fixed and maintained within a population with 3 distinct outcomes: neofunctionalization, subfunctionalization, and conservation of function. Subsequent researchers have proposed a multitude of population genetic models that lead to these outcomes, each differing largely in the role played by adaptive natural selection. In this paper, I present a nonmathematical review of these models, their predictions, and the evidence collected in support of each of them. Though the various outcomes of gene duplication are often strictly associated with the presence or absence of adaptive natural selection, I argue that determining the outcome of duplication is orthogonal to determining whether natural selection has acted. Despite an ever-growing field of research into the fate of gene duplicates, there is not yet clear evidence for the preponderance of one outcome over the others, much less evidence for the importance of adaptive or nonadaptive forces in maintaining these duplicates.

Protease gene duplication and proteolytic activity in Drosophila female reproductive tracts

Secreted proteases play integral roles in sexual reproduction in a broad range of taxa. In the genetic model Drosophila melanogaster, these molecules are thought to process peptides and activate enzymes inside female reproductive tracts, mediating critical postmating responses. A recent study of female reproductive tract proteins in the cactophilic fruit fly Drosophila arizonae, identified pervasive, lineage-specific gene duplication amongst secreted proteases. Here, we compare the evolutionary dynamics, biochemical nature, and physiological significance of secreted female reproductive serine endoproteases between D. arizonae and its congener D. melanogaster. We show that D. arizonae lower female reproductive tract (LFRT) proteins are significantly enriched for recently duplicated secreted proteases, particularly serine endoproteases, relative to D. melanogaster. Isolated lumen from D. arizonae LFRTs, furthermore, exhibits significant trypsin-like and elastase-like serine endoprotease activity, whereas no such activity is seen in D. melanogaster. Finally, trypsin- and elastase-like activity in D. arizonae female reproductive tracts is negatively regulated by mating. We propose that the intense proteolytic environment of the D. arizonae female reproductive tract relates to the extraordinary reproductive physiology of this species and that ongoing gene duplication amongst these proteases is an evolutionary consequence of sexual conflict.

Genetic architecture of hybrid male sterility in Drosophila: analysis of intraspecies variation for interspecies isolation

BACKGROUND

The genetic basis of postzygotic isolation is a central puzzle in evolutionary biology. Evolutionary forces causing hybrid sterility or inviability act on the responsible genes while they still are polymorphic, thus we have to study these traits as they arise, before isolation is complete.

METHODOLOGY/PRINCIPAL FINDINGS

Isofemale strains of D. mojavensis vary significantly in their production of sterile F(1) sons when females are crossed to D. arizonae males. We took advantage of the intraspecific polymorphism, in a novel design, to perform quantitative trait locus (QTL) mapping analyses directly on F(1) hybrid male sterility itself. We found that the genetic architecture of the polymorphism for hybrid male sterility (HMS) in the F(1) is complex, involving multiple QTL, epistasis, and cytoplasmic effects.

CONCLUSIONS/SIGNIFICANCE

The role of extensive intraspecific polymorphism, multiple QTL, and epistatic interactions in HMS in this young species pair shows that HMS is arising as a complex trait in this system. Directional selection alone would be unlikely to maintain polymorphism at multiple loci, thus we hypothesize that directional selection is unlikely to be the only evolutionary force influencing postzygotic isolation.

Molecular evolution and population genetics of two Drosophila mettleri cytochrome P450 genes involved in host plant utilization

Understanding the genetic basis of adaptation is one of the primary goals of evolutionary biology. The evolution of xenobiotic resistance in insects has proven to be an especially suitable arena for studying the genetics of adaptation, and resistant phenotypes are known to result from both coding and regulatory changes. In this study, we examine the evolutionary history and population genetics of two Drosophila mettleri cytochrome P450 genes that are putatively involved in the detoxification of alkaloids present in two of its cactus hosts: saguaro (Carnegiea gigantea) and senita (Lophocereus schottii). Previous studies demonstrated that Cyp28A1 was highly up-regulated following exposure to rotting senita tissue while Cyp4D10 was highly up-regulated following exposure to rotting saguaro tissue. Here, we show that a subset of sites in Cyp28A1 experienced adaptive evolution specifically in the D. mettleri lineage. Moreover, neutrality tests in several populations were also consistent with a history of selection on Cyp28A1. In contrast, we did not find evidence for positive selection on Cyp4D10, although this certainly does not preclude its involvement in host plant use. A surprising result that emerged from our population genetic analyses was the presence of significant genetic differentiation between flies collected from different host plant species (saguaro and senita) at Organ Pipe National Monument, Arizona, USA. This preliminary evidence suggests that D. mettleri may have evolved into distinctive host races that specialize on different hosts, a possibility that warrants further investigation.

The molecular basis of host adaptation in cactophilic Drosophila: molecular evolution of a glutathione S-transferase gene (GstD1) in Drosophila mojavensis

Drosophila mojavensis is a cactophilic fly endemic to the northwestern deserts of North America. This species includes four genetically isolated cactus host races each individually specializing on the necrotic tissues of a different cactus species. The necrosis of each cactus species provides the resident D. mojavensis populations with a distinct chemical environment. A previous investigation of the role of transcriptional variation in the adaptation of D. mojavensis to its hosts produced a set of candidate loci that are differentially expressed in response to host shifts, and among them was glutathione S-transferase D1 (GstD1). In both D. melanogaster and Anopheles gambiae, GstD1 has been implicated in the resistance of these species to the insecticide dichloro-diphenyl-trichloroethane (DDT). The pattern of sequence variation of the GstD1 locus from all four D. mojavensis populations, D. arizonae (sister species), and D. navojoa (outgroup) has been examined. The data suggest that in two populations of D. mojavensis GstD1 has gone through a period of adaptive amino acid evolution. Further analyses indicate that of the seven amino acid fixations that occurred in the D. mojavensis lineage, two of them occur in the active site pocket, potentially having a significant effect on substrate specificity and in the adaptation to alternative cactus hosts.

Multilocus nuclear sequences reveal intra- and interspecific relationships among chromosomally polymorphic species of cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae, a pair of sibling species endemic to North America, constitute an important model system to study ecological genetics and the evolution of reproductive isolation. This species pair can produce fertile hybrids in some crosses and are sympatric in a large part of their ranges. Despite the potential for hybridization in nature, however, evidence of introgression has not been rigorously sought. Further, the evolutionary relationships within and among the geographically distant populations of the two species have not been characterized in detail using high-resolution molecular studies. Both species have six chromosomes: five large acrocentrics and one 'dot' chromosome. Fixed inversion differences between the species exist in three chromosomes (X, 2 and 3) while three are colinear (4, 5 and 6), suggesting that were introgression to occur, it would be most likely in the colinear chromosomes. We utilized nucleotide sequence variation at multiple loci on five chromosomes to test for evidence of introgression, and to test various scenarios for the evolutionary relationships of these two species and their populations. While we do not find evidence of recent introgression, loci in the colinear chromosomes appear to have participated in exchange in the past. We also found considerable population structure within both species. The level of differentiation discovered among D. arizonae populations was unexpectedly high and suggests that its populations, as well as those of D. mojavensis, may be themselves undergoing incipient speciation and merit further attention.

Gene duplication and adaptive evolution of digestive proteases in Drosophila arizonae female reproductive tracts

It frequently has been postulated that intersexual coevolution between the male ejaculate and the female reproductive tract is a driving force in the rapid evolution of reproductive proteins. The dearth of research on female tracts, however, presents a major obstacle to empirical tests of this hypothesis. Here, we employ a comparative EST approach to identify 241 candidate female reproductive proteins in Drosophila arizonae, a repleta group species in which physiological ejaculate–female coevolution has been documented. Thirty-one of these proteins exhibit elevated amino acid substitution rates, making them candidates for molecular coevolution with the male ejaculate. Strikingly, we also discovered 12 unique digestive proteases whose expression is specific to the D. arizonae lower female reproductive tract. These enzymes belong to classes most commonly found in the gastrointestinal tracts of a diverse array of organisms. We show that these proteases are associated with recent, lineage-specific gene duplications in the Drosophila repleta species group, and exhibit strong signatures of positive selection. Observation of adaptive evolution in several female reproductive tract proteins indicates they are active players in the evolution of reproductive tract interactions. Additionally, pervasive gene duplication, adaptive evolution, and rapid acquisition of a novel digestive function by the female reproductive tract points to a novel coevolutionary mechanism of ejaculate–female interaction.

In a broad range of organisms, including humans, molecular interactions between the male ejaculate and the female reproductive tract play integral roles in sexual reproduction. Although these interactions are essential, the biochemical composition of the male ejaculate can change rapidly over short evolutionary time periods. It is often hypothesized that this rapid evolution reflects a coevolutionary relationship with the female reproductive tract. The paucity of research on females, however, presents a formidable challenge to empirical tests of this hypothesis. In this study, we sought to identify proteins in the female reproductive tracts of D. arizonae that may be interacting or coevolving with the male ejaculate. Unexpectedly, we discovered that D. arizonae females produce an array of “digestive” enzymes in their reproductive tracts. These classes of enzymes are normally found in the gut, where they degrade ingested food for nutritional uptake. In D. arizonae, these enzymes have resulted from recent gene duplications, and natural selection has caused rapid and radical changes in their amino acid sequences. We propose that this pattern of duplication and diversification reflects the “female side” of a coevolutionary relationship with the male ejaculate. Exploring the “male side” of this relationship is an important avenue for future research.

Functional genomics of cactus host shifts in Drosophila mojavensis

Understanding the genetic basis of adaptation to novel environments remains one of the major challenges confronting evolutionary biologists. While newly developed genomic approaches hold considerable promise for addressing this overall question, the relevant tools have not often been available in the most ecologically interesting organisms. Our study organism, Drosophila mojavensis, is a cactophilic Sonoran Desert endemic utilizing four different cactus hosts across its geographical range. Its well-known ecology makes it an attractive system in which to study the evolution of gene expression during adaptation. As a cactophile, D. mojavensis oviposits in the necrotic tissues of cacti, therefore exposing larvae and even adults to the varied and toxic compounds of rotting cacti. We have developed a cDNA microarray of D. mojavensis to examine gene expression associated with cactus host use. Using a population from the Baja California population we examined gene expression differences of third instar larvae when reared in two chemically distinct cactus hosts, agria (Stenocereus gummosus, native host) vs. organpipe (Stenocereus thurberi, alternative host). We have observed differential gene expression associated with cactus host use in genes involved in metabolism and detoxification.

Microsatellite variation among diverging populations of Drosophila mojavensis

Divergence and speciation may occur by various means, depending on the particular history, selective environments, and genetic composition of populations. In Drosophila mojavensis, a good model of incipient speciation, understanding the population genetic structure within this group facilitates our ability to understand the context in which reproductive isolation among populations is developing. Here we report the genetic structure and relationships of D. mojavensis populations at nuclear loci. We surveyed 29 populations throughout the distribution of D. mojavensis for four microsatellite loci to differentiation among populations of this species. These loci reveal four distinct geographical regions of D. mojavensis populations in the south-western United States and north-western Mexico--(i) Baja California peninsula (Baja), (ii) Sonora, Mexico-southern Arizona, United States (Sonora), (iii) Mojave Desert and Grand Canyon (Mojave), and (iv) Santa Catalina Island (Catalina). While all regions show strong isolation, Mojave and Catalina are highly diverged from other regions. Within any region, populations are largely homogenous over broad geographical distances. Based on the population structure, we find clear geographical barriers to gene flow appear to have a strong effect in isolating populations across regions for this species.

Sympatry, allopatry and sexual isolation between Drosophila mojavensis and D. arizonae

Populations of the North American cactophilic fruitfly Drosophila mojavensis and its sibling species D. arizonae exist both in sympatry and in allopatry. Females of D. arizonae, regardless of their population of origin, are effectively completely isolated behaviorally from D. mojavensis males. On the other hand, females of D. mojavensis from the sympatric populations in Sonora, Mexico exhibit significantly stronger premating isolation from D. arizonae males than do D. mojavensis females from allopatric populations from the Baja California peninsula. Earlier studies interpreted these limited observations as support for reinforcement. Since the time of those studies, additional allopatric populations of D. mojavensis have been collected from southern California and from Santa Catalina Island, off the coast of southern California. Here, we tested the prediction that if sympatry is in fact associated with increased isolation in D. mojavensis, these additional allopatric populations also should show, relative to the sympatric ones, less isolation from D. arizonae. Our results are consistent with this prediction and suggest that isolation is in fact stronger in sympatry.

Ribosomal intergenic spacer (IGS) length variation across the Drosophilinae (Diptera: Drosophilidae)

Background

The intergenic spacer of the ribosomal genes in eukaryotes (IGS) contains duplications of the core transcription promoter. The number of these duplicated promoters, as measured by the IGS length, appears to be correlated with growth rate and development time in several distantly related taxa. In the present study, we examined IGS length variation across a number of species of Drosophila to determine the amount of variation in this trait across different evolutionary time scales. Furthermore, we compared the usefulness of two methods commonly used to determine IGS length: Southern Blot Hybridization (SB) and Polymerase Chain Reaction (PCR).

Results

Our results show broad variation in IGS length across the genus Drosophila, but closely related species had similar IGS lengths. Our results also suggest that PCR tends to underestimate the true IGS size when the size is greater than 5 kb, and that this degree of underestimation is greater as the IGS size increases.

Conclusion

Broad variation in IGS length occurs across large evolutionary divergences in the subfamily Drosophilinae. Although average IGS length has been shown to evolve rapidly under artificial selection, closely related taxa generally have similar average IGS lengths. Our comparison of methods suggests that without previous knowledge of the DNA sequence of the IGS and flanking regions, both methods be used to accurately measure IGS length.

Activity variation in alcohol dehydrogenase paralogs is associated with adaptation to cactus host use in cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae are species of cactophilic flies that share a recent duplication of the alcohol dehydrogenase (Adh) locus. One paralog (Adh-2) is expressed in adult tissues and the other (Adh-1) in larvae and ovaries. Enzyme activity measurements of the ADH-2 amino acid polymorphism in D. mojavensis suggest that the Fast allozyme allele has a higher activity on 2-propanol than 1-propanol. The Fast allele was found at highest frequency in populations that utilize hosts with high proportions of 2-propanol, while the Slow allele is most frequent in populations that utilize hosts with high proportions of 1-propanol. This suggests that selection for ADH-2 allozyme alleles with higher activity on the most abundant alcohols is occurring in each D. mojavensis population. In the other paralog, ADH-1, significant differences between D. mojavensis and D. arizonae are associated with a previously shown pattern of adaptive protein evolution in D. mojavensis. Examination of protein sequences showed that a large number of amino acid fixations between the paralogs have occurred in catalytic residues. These changes are potentially responsible for the significant difference in substrate specificity between the paralogs. Both functional and sequence variation within and between paralogs suggests that Adh has played an important role in the adaptation of D. mojavensis and D. arizonae to their cactophilic life.