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

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

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

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.

Gene expression and alternative splicing dynamics are perturbed in female head transcriptomes following heterospecific copulation

BACKGROUND

Despite the growing interest in the female side of copulatory interactions, the roles played by differential expression and alternative splicing mechanisms of pre-RNA on tissues outside of the reproductive tract have remained largely unknown. Here we addressed these questions in the context of con- vs heterospecific matings between Drosophila mojavensis and its sister species, D. arizonae. We analyzed transcriptional responses in female heads using an integrated investigation of genome-wide patterns of gene expression, including differential expression (DE), alternative splicing (AS) and intron retention (IR).

RESULTS

Our results indicated that early transcriptional responses were largely congruent between con- and heterospecific matings but are substantially perturbed over time. Conspecific matings induced functional pathways related to amino acid balance previously associated with the brain's physiology and female postmating behavior. Heterospecific matings often failed to activate regulation of some of these genes and induced expression of additional genes when compared with those of conspecifically-mated females. These mechanisms showed functional specializations with DE genes mostly linked to pathways of proteolysis and nutrient homeostasis, while AS genes were more related to photoreception and muscle assembly pathways. IR seems to play a more general role in DE regulation during the female postmating response.

CONCLUSIONS

We provide evidence showing that AS genes substantially perturbed by heterospecific matings in female heads evolve at slower evolutionary rates than the genome background. However, DE genes evolve at evolutionary rates similar, or even higher, than those of male reproductive genes, which highlights their potential role in sexual selection and the evolution of reproductive barriers.

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.

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.

Genome-Wide Patterns of Sequence Divergence of Protein-Coding Genes Between Drosophila buzzatii and D. mojavensis

Evolutionary rates for protein-coding genes are determined not only by natural selection but also by multiple genomic factors including mutation rates, recombination, gene expression levels, and chromosomal location. To investigate the joint effects of different genomic determinants on protein evolution, we compared the coding sequences of 9017 single-copy orthologs between 2 cactophilic species from the Drosophila subgenus, Drosophila mojavensis and D. buzzatii, whose genomes have been previously sequenced. We assessed the impact of 7 genomic determinants, that is, chromosome type, recombination, chromosomal inversions, expression breadth, expression level, gene length, and the number of exons, on divergence rates of protein-coding genes to understand patterns of evolutionary variation. Integrative analysis of these factors revealed that 1) X-linked and autosomal genes evolve at significantly different rates in agreement with the faster-X hypothesis, 2) genes located on the dot chromosome and pericentromeric regions have higher divergence rates, 3) genes located at chromosomes with more fixed inversions have higher pairwise divergence than those located at nearly collinear chromosomes, and 4) gene expression patterns can be considered the strongest determinant of protein evolution. In addition, the number of exons and protein length had a significant effect on pairwise divergence at synonymous sites. All in all, our results show the relative importance of each genomic factor on the rates of protein evolution and functional constraint in these 2 cactophilic Drosophila species.

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.

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.

Effects of temperature on transcriptome and cuticular hydrocarbon expression in ecologically differentiated populations of desert Drosophila

We assessed the effects of temperature differences on gene expression using whole-transcriptome microarrays and cuticular hydrocarbon variation in populations of cactophilic Drosophila mojavensis. Four populations from Baja California and mainland Mexico and Arizona were each reared on two different host cacti, reared to sexual maturity on laboratory media, and adults were exposed for 12 hr to 15, 25, or 35°C. Temperature differences influenced the expression of 3,294 genes, while population differences and host plants affected >2,400 each in adult flies. Enriched, functionally related groups of genes whose expression changed at high temperatures included heat response genes, as well as genes affecting chromatin structure. Gene expression differences between mainland and peninsular populations included genes involved in metabolism of secondary compounds, mitochondrial activity, and tRNA synthases. Flies reared on the ancestral host plant, pitaya agria cactus, showed upregulation of genes involved in metabolism, while flies reared on organ pipe cactus had higher expression of DNA repair and chromatin remodeling genes. Population × environment (G × E) interactions had widespread effects on the transcriptome where population × temperature interactions affected the expression of >5,000 orthologs, and there were >4,000 orthologs that showed temperature × host plant interactions. Adults exposed to 35°C had lower amounts of most cuticular hydrocarbons than those exposed to 15 or 25°C, including abundant unsaturated alkadienes. For insects adapted to different host plants and climatic regimes, our results suggest that temperature shifts associated with climate change have large and significant effects on transcriptomes of genetically differentiated natural populations.

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.

Genome-wide tests for introgression between cactophilic Drosophila implicate a role of inversions during speciation

Models of speciation‐with‐gene‐flow have shown that the reduction in recombination between alternative chromosome arrangements can facilitate the fixation of locally adaptive genes in the face of gene flow and contribute to speciation. However, it has proven frustratingly difficult to show empirically that inversions have reduced gene flow and arose during or shortly after the onset of species divergence rather than represent ancestral polymorphisms. Here, we present an analysis of whole genome data from a pair of cactophilic fruit flies, Drosophila mojavensis and D. arizonae, which are reproductively isolated in the wild and differ by several large inversions on three chromosomes. We found an increase in divergence at rearranged compared to colinear chromosomes. Using the density of divergent sites in short sequence blocks we fit a series of explicit models of species divergence in which gene flow is restricted to an initial period after divergence and may differ between colinear and rearranged parts of the genome. These analyses show that D. mojavensis and D. arizonae have experienced postdivergence gene flow that ceased around 270 KY ago and was significantly reduced in chromosomes with fixed inversions. Moreover, we show that these inversions most likely originated around the time of species divergence which is compatible with theoretical models that posit a role of inversions in speciation with gene flow.

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.

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.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

Species coherence in the face of karyotype diversification in holocentric organisms: the case of a cytogenetically variable sedge (Carex scoparia, Cyperaceae)

BACKGROUND AND AIMS

The sedge genus Carex, the most diversified angiosperm genus of the northern temperate zone, is renowned for its holocentric chromosomes and karyotype variability. The genus exhibits high variation in chromosome numbers both among and within species. Despite the possibility that this chromosome evolution may play a role in the high species diversity of Carex, population-level patterns of molecular and cytogenetic differentiation in the genus have not been extensively studied.

METHODS

Microsatellite variation (11 loci, 461 individuals) and chromosomal diversity (82 individuals) were investigated in 22 Midwestern populations of the North American sedge Carex scoparia and two Northeastern populations.

KEY RESULTS

Among Midwestern populations, geographic distance is the most important predictor of genetic differentiation. Within populations, inbreeding is high and chromosome variation explains a significant component of genetic differentiation. Infrequent dispersal among populations separated by >100 km explains an important component of molecular genetic and cytogenetic diversity within populations. However, karyotype variation and correlation between genetic and chromosomal variation persist within populations even when putative migrants based on genetic data are excluded.

CONCLUSIONS

These findings demonstrate dispersal and genetic connectivity among widespread populations that differ in chromosome numbers, explaining the phenomenon of genetic coherence in this karyotypically diverse sedge species. More generally, the study suggests that traditional sedge taxonomic boundaries demarcate good species even when those species encompass a high range of chromosomal diversity. This finding is important evidence as we work to document the limits and drivers of biodiversity in one of the world's largest angiosperm genera.

Density dynamics of diverse Spiroplasma strains naturally infecting different species of Drosophila

Facultative heritable bacterial endosymbionts can have dramatic effects on their hosts, ranging from mutualistic to parasitic. Within-host bacterial endosymbiont density plays a critical role in maintenance of a symbiotic relationship, as it can affect levels of vertical transmission and expression of phenotypic effects, both of which influence the infection prevalence in host populations. Species of genus Drosophila are infected with Spiroplasma, whose characterized phenotypic effects range from that of a male-killing reproductive parasite to beneficial defensive endosymbiont. For many strains of Spiroplasma infecting at least 17 species of Drosophila, however, the phenotypic effects are obscure. The infection prevalence of these Spiroplasma vary within and among Drosophila species, and little is known about the within-host density dynamics of these diverse strains. To characterize the patterns of Spiroplasma density variation among Drosophila we used quantitative PCR to assess bacterial titer at various life stages of three species of Drosophila naturally-infected with two different types of Spiroplasma. For naturally infected Drosophila species we found that non-male-killing infections had consistently lower densities than the male-killing infection. The patterns of Spiroplasma titer change during aging varied among Drosophila species infected with different Spiroplasma strains. Bacterial density varied within and among populations of Drosophila, with individuals from the population with the highest prevalence of infection having the highest density. This density variation underscores the complex interaction of Spiroplasma strain and host genetic background in determining endosymbiont density.

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.

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.

Diversity-enhancing selection acts on a female reproductive protease family in four subspecies of Drosophila mojavensis

Protein components of the Drosophila male ejaculate are critical modulators of reproductive success, several of which are known to evolve rapidly. Recent evidence of adaptive evolution in female reproductive tract proteins suggests this pattern may reflect sexual selection at the molecular level. Here we explore the evolutionary dynamics of a five-paralog gene family of female reproductive proteases within geographically isolated subspecies of Drosophila mojavensis. Remarkably, four of five paralogs show exceptionally low differentiation between subspecies and unusually structured haplotypes that suggest the retention of old polymorphisms. These gene genealogies are accompanied by deviations from neutrality consistent with diversifying selection. While diversifying selection has been observed among the reproductive molecules of mammals and marine invertebrates, our study provides the first evidence of this selective regime in any Drosophila reproductive protein, male or female.

Mitotic and polytene chromosomes analysis of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

The Oriental fruit fly, Batrocera dorsalis s.s. (Hendel) is one of the most destructive agricultural pests, belonging to a large group of difficult to distinguish morphologically species, referred as the B. dorsalis complex. We report here a cytogenetic analysis of two laboratory strains of the species and provide a photographic polytene chromosome map from larval salivary glands. The mitotic complement consists of six chromosome pairs including a heteromorphic sex (XX/XY) chromosome pair. Analysis of the polytene complement has shown a total of five polytene chromosomes (10 polytene arms) that correspond to the five autosomes. The most important landmarks of each polytene chromosome and characteristic asynapsis at a specific chromosomal region are presented and discussed. Chromosomal homology between B. dorsalis and Ceratitis capitata has been determined by comparing chromosome banding patterns. The detection of chromosome inversions in both B. dorsalis strains is shown and discussed. Our results show that the polytene maps presented here are suitable for cytogenetic analysis of this species and can be used for comparative studies among species of the Tephritidae family. They also provide a diagnostic tool that could accelerate species identification within the B. dorsalis complex and could shed light on the ongoing speciation in this complex. Polytene chromosome maps can facilitate the development of biological control methods and support the genome mapping project of the species that is currently in progress.

Islands of speciation or mirages in the desert? Examining the role of restricted recombination in maintaining species

Over the past decade, many studies documented high genetic divergence between closely related species in genomic regions experiencing restricted recombination in hybrids, such as within chromosomal rearrangements or areas adjacent to centromeres. Such regions have been called 'islands of speciation' because of their presumed role in maintaining the integrity of species despite gene flow elsewhere in the genome. Here, we review alternative explanations for such patterns. Segregation of ancestral variation or artifacts of nucleotide diversity within species can readily lead to higher F(ST) in regions of restricted recombination than other parts of the genome, even in the complete absence of interspecies gene flow, and thereby cause investigators to erroneously conclude that islands of speciation exist. We conclude by discussing strengths and weaknesses of various means for testing the role of restricted recombination in maintaining species.

Metabolic pools differ among ecologically diverse Drosophila species

Studies of the genetic mechanisms underlying metabolic storage have focused on a few model organisms. Although very fruitful, these studies have not allowed for the examination of mechanisms across a phylogenetic spectrum. The exploration of natural patterns of metabolic pool size variation across species will help us to better understand the genetics of metabolic adaptation. We examined the metabolic pools size (triglyceride, glycogen and protein) at two ages in 12 Drosophila species with distinctly different ecologies for which complete genome sequences (for 11 of the 12 species) are known. Overall, there were significant differences across species for all three pools, while age and sex appear to affect some metabolic pools more than others. After correcting for the phylogenetic relatedness of the species used, we observed no association between triglyceride and glycogen content. Although within species these two pools sometimes are correlated, at a larger phylogenetic scale control of triglyceride and glycogen contents may have been shaped independently by natural selection.

The genetics and ecology of reinforcement: implications for the evolution of prezygotic isolation in sympatry and beyond

Reinforcement, the evolution of prezygotic reproductive barriers by natural selection in response to maladaptive hybridization, is one of the most debated processes in speciation. Critics point to "fatal" conceptual flaws for sympatric evolution of prezygotic isolation, but recent theoretical and empirical work on genetics and ecology of reinforcement suggests that such criticisms can be overcome. New studies provide evidence for reinforcement in frogs, fish, insects, birds, and plants. While such evidence lays to rest the argument over reinforcement's existence, our understanding remains incomplete. We lack data on (1) the genetic basis of female preferences and the links between genetics of pre- and postzygotic isolation, (2) the ecological basis of reproductive isolation, (3) connections between prezygotic isolation between species and within-species sexual selection (potentially leading to a "cascade" of effects on reproductive isolation), (4) the role of habitat versus mate preference in reinforcement, and (5) additional detailed comparative studies. Here, we review data on these issues and highlight why they are important for understanding speciation.

Post-zygotic isolation in cactophilic Drosophila: larval viability and adult life-history traits of D. mojavensis/D. arizonae hybrids

Drosophila mojavensis and Drosophila arizonae are cactophilic flies that have been used extensively in speciation studies. Incomplete premating isolation, evidence of reinforcement, and a lack of recent introgression between these species point to a potentially important role for post-zygotic isolating barriers in this system. Other than hybrid male sterility, however, post-zygotic isolation between D. mojavensis and D. arizonae has received little attention. In this study, we examined viability and life-history traits of D. mojavensis/D. arizonae F(1) hybrids from sympatric crosses. Specifically, we reared hybrids and purebreds on the natural host cacti of each parental species and compared viability, development time, thorax length, and desiccation resistance between hybrids and purebreds. Interestingly, hybrid females from both crosses performed similarly or even better than purebred females. In contrast, hybrid sons of D. arizonae mothers, in addition to being sterile, had shorter average thorax length than males of both parental species, and hybrid males from both crosses had substantially lower desiccation resistance than D. mojavensis males. The probable cost to hybridization for D. mojavensis females resulting from reduced desiccation resistance of hybrid sons may have been an important selective factor in the history of reinforcement for crosses involving these females.

Multiple introductions of the Spiroplasma bacterial endosymbiont into Drosophila

Bacterial endosymbionts are common in insects and can have dramatic effects on their host's evolution. So far, the only heritable symbionts found in Drosophila have been Wolbachia and Spiroplasma. While the incidence and effects of Wolbachia have been studied extensively, the prevalence and significance of Spiroplasma infections in Drosophila are less clear. These small, gram-positive, helical bacteria infect a diverse array of plant and arthropod hosts, conferring a variety of fitness effects. Male-killing Spiroplasma are known from certain Drosophila species; however, in others, Spiroplasma appear not to affect sex ratio. Previous studies have identified different Spiroplasma haplotypes in Drosophila populations, although no extensive surveys have yet been reported. We used a multilocus sequence analysis to reconstruct a robust Spiroplasma endosymbiont phylogeny, assess genetic diversity, and look for evidence of recombination. Six loci were sequenced from over 65 Spiroplasma-infected individuals from nine different Drosophila species. Analysis of these sequences reveals at least five separate introductions of four phylogenetically distinct Spiroplasma haplotypes, indicating that more extensive sampling will likely reveal an even greater Spiroplasma endosymbiont diversity. Patterns of variation in Drosophila mitochondrial haplotypes in Spiroplasma-infected and uninfected flies imply imperfect vertical transmission in host populations and possible horizontal transmission.

Duplication, selection and gene conversion in a Drosophila mojavensis female reproductive protein family

Protein components of the Drosophila male ejaculate, several of which evolve rapidly, are critical modulators of reproductive success. Recent studies of female reproductive tract proteins indicate they also are extremely divergent between species, suggesting that reproductive molecules may coevolve between the sexes. Our current understanding of intersexual coevolution, however, is severely limited by the paucity of genetic and evolutionary studies on the female molecules involved. Physiological evidence of ejaculate-female coadaptation, paired with a promiscuous mating system, makes Drosophila mojavensis an exciting model system in which to study the evolution of reproductive proteins. Here we explore the evolutionary dynamics of a five-paralog gene family of female reproductive proteases within populations of D. mojavensis and throughout the repleta species group. We show that the proteins have experienced ongoing gene duplication and adaptive evolution and further exhibit dynamic patterns of pseudogenation, copy number variation, gene conversion, and selection within geographically isolated populations of D. mojavensis. The integration of these patterns in a single gene family has never before been documented in a reproductive protein.

Genetic diversification and demographic history of the cactophilic pseudoscorpion Dinocheirus arizonensis from the Sonoran Desert

Sequence data from a segment of the mitochondrial cytochrome c oxidase subunit I (COI) gene were used to examine phylogenetic relationships, estimate gene flow and infer demographic history of the cactophilic chernetid pseudoscorpion, Dinocheirus arizonensis (Banks), from the Sonoran Desert. Phylogenetic trees resolved two clades of D. arizonensis, one from mainland Sonora, Mexico and southern Arizona (clade I) and the other from the Baja California peninsula and southern Arizona (clade II). The two clades were separated by a mean genetic distance (d) of approximately 2.6%. Hierarchical analysis of molecular variance indicated highly significant population structuring in D. arizonensis (overall Phi(ST)=0.860; P<0.0001), with 80% of the genetic variation distributed among the two clades. Most pairwise comparisons of Phi(ST) among populations within each clade, however, were not significant. The results suggest that phoretic dispersal on vagile cactophilic insects such as the neriid cactus fly Odontoloxozus longicornis (Coquillett) provides sufficient gene flow to offset the accumulation of unique haplotypes within each clade of the non-vagile pseudoscorpion. Preliminary results on dispersal capability of O. longicornis were consistent with this conclusion. Tests designed to reconstruct demographic history from sequence data indicated that both clades of D. arizonensis, as well as O. longicornis, have experienced historical population expansions. Potential barriers to gene flow that may have led to genetic isolation and diversification in clades I and II of D. arizonensis are discussed.

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

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.

Molecular demographic history of the annual sunflowers Helianthus annuus and H. petiolaris--large effective population sizes and rates of long-term gene flow

Hybridization between distinct species may lead to introgression of genes across species boundaries, and this pattern can potentially persist for extended periods as long as selection at some loci or genomic regions prevents thorough mixing of gene pools. However, very few reliable estimates of long-term levels of effective migration are available between hybridizing species throughout their history. Accurate estimates of divergence dates and levels of gene flow require data from multiple unlinked loci as well as an analytical framework that can distinguish between lineage sorting and gene flow and incorporate the effects of demographic changes within each species. Here we use sequence data from 18 anonymous nuclear loci in two broadly sympatric sunflower species, Helianthus annuus and H. petiolaris, analyzed within an "isolation with migration" framework to make genome-wide estimates of the ages of these two species, long-term rates of gene flow between them, and effective population sizes and historical patterns of population growth. Our results indicate that H. annuus and H. petiolaris are approximately one million years old and have exchanged genes at a surprisingly high rate (long-term N(ef)m estimates of approximately 0.5 in each direction), with somewhat higher rates of introgression from H. annuus into H. petiolaris than vice versa. In addition, each species has undergone dramatic population expansion since divergence, and both species have among the highest levels of genetic diversity reported for flowering plants. Our results provide the most comprehensive estimate to date of long-term patterns of gene flow and historical demography in a nonmodel plant system, and they indicate that species integrity can be maintained even in the face of extensive gene flow over a prolonged period.

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.