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

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.

Divergent evolutionary trajectories shape the postmating transcriptional profiles of conspecifically and heterospecifically mated cactophilic Drosophila females

Postmating-prezygotic (PMPZ) reproductive isolation is hypothesized to result from divergent coevolutionary trajectories of sexual selection and/or sexual conflict in isolated populations. However, the genetic basis of PMPZ incompatibilities between species is poorly understood. Here, we use a comparative framework to compare global gene expression in con- and heterospecifically mated Drosophila mojavensis and D. arizonae female reproductive tracts. We find striking divergence between the species in the female postmating transcriptional response to conspecific mating, including differences in differential expression (DE), alternative splicing (AS), and intron retention (IR). As predicted, heterospecific matings produce disrupted transcriptional profiles, but the overall patterns of misregulation are different between the reciprocal crosses. Moreover, we find a positive correlation between postmating transcriptional divergence between species and levels of transcriptional disruption in heterospecific crosses. This result indicates that mating responsive genes that have diverged more in expression also have more disrupted transcriptional profiles in heterospecifically mated females. Overall, our results provide insights into the evolution of PMPZ isolation and lay the foundation for future studies aimed at identifying specific genes involved in PMPZ incompatibilities and the evolutionary forces that have contributed to their divergence in closely related species.

Comparison of global gene expression patterns in con- and heterospecifically mated Drosophila mojavensis and Drosophila arizonae suggest that mating-responsive genes with divergent expression also exhibit more disrupted transcriptional profiles in heterospecifically mated females, providing further insight into the evolution of postmating-prezygotic reproductive isolation.

Evolution of the odorant-binding protein gene family in Drosophila

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

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

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

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

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.

Seven Questions on the Chemical Ecology and Neurogenetics of Resource-Mediated Speciation

Adaptation to different environments can result in reproductive isolation between populations and the formation of new species. Food resources are among the most important environmental factors shaping local adaptation. The chemosensory system, the most ubiquitous sensory channel in the animal kingdom, not only detects food resources and their chemical composition, but also mediates sexual communication and reproductive isolation in many taxa. Chemosensory divergence may thus play a crucial role in resource-mediated adaptation and speciation. Understanding how the chemosensory system can facilitate resource-mediated ecological speciation requires integrating mechanistic studies of the chemosensory system with ecological studies, to link the genetics and physiology of chemosensory properties to divergent adaptation. In this review, we use examples of insect research to present seven key questions that can be used to understand how the chemosensory system can facilitate resource-mediated ecological speciation in consumer populations.

Mate discrimination among subspecies through a conserved olfactory pathway

Communication mechanisms underlying the sexual isolation of species are poorly understood. Using four subspecies of Drosophila mojavensis as a model, we identify two behaviorally active, male-specific pheromones. One functions as a conserved male antiaphrodisiac in all subspecies and acts via gustation. The second induces female receptivity via olfaction exclusively in the two subspecies that produce it. Genetic analysis of the cognate receptor for the olfactory pheromone indicates an important role for this sensory pathway in promoting sexual isolation of subspecies, in combination with auditory signals. Unexpectedly, the peripheral sensory pathway detecting this pheromone is conserved molecularly, physiologically, and anatomically across subspecies. These observations imply that subspecies-specific behaviors arise from differential interpretation of the same peripheral cue, reminiscent of sexually conserved detection but dimorphic interpretation of male pheromones in Drosophila melanogaster. Our results reveal that, during incipient speciation, pheromone production, detection, and interpretation do not necessarily evolve in a coordinated manner.