Evolutionary genomics of host plant adaptation: insights from Drosophila

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

Diverse signatures of convergent evolution in cactus-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently approximately 17 times. Using a machine learning-based approach, we further found that cactophily can be predicted with 76% accuracy from both functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which we found to be likely associated with altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved independently through disparate molecular mechanisms. Notably, we found that multiple cactophilic species and their close relatives have been reported as emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-might preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high-throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Characterization of the chemoreceptor repertoire of a highly specialized fly with comparisons to other Drosophila species

To explore the diversity of scenarios in nature, animals have evolved tools to interact with different environmental conditions. Chemoreceptors are an important interface component and among them, olfactory receptors (ORs) and gustatory receptors (GRs) can be used to find food and detect healthy resources. Drosophila is a model organism in many scientific fields, in part due to the diversity of species and niches they occupy. The contrast between generalists and specialists Drosophila species provides an important model for studying the evolution of chemoreception. Here, we compare the repertoire of chemoreceptors of different species of Drosophila with that of D. incompta, a highly specialized species whose ecology is restricted to Cestrum flowers, after reporting the preferences of D. incompta to the odor of Cestrum flowers in olfactory tests. We found evidence that the chemoreceptor repertoire in D. incompta is smaller than that presented by species in the Sophophora subgenus. Similar patterns were found in other non-Sophophora species, suggesting the presence of underlying phylogenetic trends. Nevertheless, we also found autapomorphic gene losses and detected some genes that appear to be under positive selection in D. incompta, suggesting that the specific lifestyle of these flies may have shaped the evolution of individual genes in each of these gene families.

Elevated sleep quota in a stress-resilient Drosophila species

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

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

Diverse signatures of convergent evolution in cacti-associated yeasts

Many distantly related organisms have convergently evolved traits and lifestyles that enable them to live in similar ecological environments. However, the extent of phenotypic convergence evolving through the same or distinct genetic trajectories remains an open question. Here, we leverage a comprehensive dataset of genomic and phenotypic data from 1,049 yeast species in the subphylum Saccharomycotina (Kingdom Fungi, Phylum Ascomycota) to explore signatures of convergent evolution in cactophilic yeasts, ecological specialists associated with cacti. We inferred that the ecological association of yeasts with cacti arose independently ~17 times. Using machine-learning, we further found that cactophily can be predicted with 76% accuracy from functional genomic and phenotypic data. The most informative feature for predicting cactophily was thermotolerance, which is likely associated with duplication and altered evolutionary rates of genes impacting the cell envelope in several cactophilic lineages. We also identified horizontal gene transfer and duplication events of plant cell wall-degrading enzymes in distantly related cactophilic clades, suggesting that putatively adaptive traits evolved through disparate molecular mechanisms. Remarkably, multiple cactophilic lineages and their close relatives are emerging human opportunistic pathogens, suggesting that the cactophilic lifestyle-and perhaps more generally lifestyles favoring thermotolerance-may preadapt yeasts to cause human disease. This work underscores the potential of a multifaceted approach involving high throughput genomic and phenotypic data to shed light onto ecological adaptation and highlights how convergent evolution to wild environments could facilitate the transition to human pathogenicity.

Elevated sleep need in a stress-resilient Drosophila species

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

Host Preferences Shown by Drosophilids (Diptera) in a Commercial Fruit and Vegetable Distribution Center Follow the Wild Neotropical Pattern

Although drosophilids have been extensively studied in laboratories worldwide, their ecology is still relatively poorly understood. This is unfortunate because some species are currently expanding their geographic distribution and infesting fruit crops. Here, we investigated the relationship between drosophilids and potential plant hosts in a commercial fruit and vegetable distribution center in the Neotropical region. We collected discarded fruits and vegetables from this commercial center during two time periods (2007-2008 and 2017-2018). Resources were weighted and individually monitored in the laboratory. The drosophilids that emerged were identified, and the relationship between them and their resources was explored. From the 99,478 kg of potential hosts collected, we identified 48 plant taxa, from which 48,894 drosophilids of 16 species emerged. On both collecting occasions, drosophilid assemblages were strongly dominated by basically the same exotic species, which explore a broader range of resources, especially those of exotic origin, when compared to neotropical drosophilids. These results are concerning because the studied site, Along with other urban markets around the world, might be acting as sources of generalist widespread species that disperse to surrounding natural vegetation and contribute to biotic homogenization.

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.

Ortholog genes from cactophilic Drosophila provide insight into human adaptation to hallucinogenic cacti

Cultural transformations of lifestyles and dietary practices have been key drivers of human evolution. However, while most of the evidence of genomic adaptations is related to the hunter-gatherer transition to agricultural societies, little is known on the influence of other major cultural manifestations. Shamanism is considered the oldest religion that predominated throughout most of human prehistory and still prevails in many indigenous populations. Several lines of evidence from ethno-archeological studies have demonstrated the continuity and importance of psychoactive plants in South American cultures. However, despite the well-known importance of secondary metabolites in human health, little is known about its role in the evolution of ethnic differences. Herein, we identified candidate genes of adaptation to hallucinogenic cactus in Native Andean populations with a long history of shamanic practices. We used genome-wide expression data from the cactophilic fly Drosophila buzzatii exposed to a hallucinogenic columnar cactus, also consumed by humans, to identify ortholog genes exhibiting adaptive footprints of alkaloid tolerance. Genomic analyses in human populations revealed a suite of ortholog genes evolving under recent positive selection in indigenous populations of the Central Andes. Our results provide evidence of selection in genetic variants related to alkaloids toxicity, xenobiotic metabolism, and neuronal plasticity in Aymara and Quechua populations, suggesting a possible process of gene-culture coevolution driven by religious practices.

Transcriptional responses are oriented towards different components of the rearing environment in two Drosophila sibling species

BACKGROUND

The chance to compare patterns of differential gene expression in related ecologically distinct species can be particularly fruitful to investigate the genetics of adaptation and phenotypic plasticity. In this regard, a powerful technique such as RNA-Seq applied to ecologically amenable taxa allows to address issues that are not possible in classic model species. Here, we study gene expression profiles and larval performance of the cactophilic siblings Drosophila buzzatii and D. koepferae reared in media that approximate natural conditions and evaluate both chemical and nutritional components of the diet. These closely related species are complementary in terms of host-plant use since the primary host of one is the secondary of the other. D. koepferae is mainly a columnar cactus dweller while D. buzzatii prefers Opuntia hosts.

RESULTS

Our comparative study shows that D. buzzatii and D. koepferae have different transcriptional strategies to face the challenges posed by their natural resources. The former has greater transcriptional plasticity, and its response is mainly modulated by alkaloids of its secondary host, while the latter has a more canalized genetic response, and its transcriptional plasticity is associated with the cactus species.

CONCLUSIONS

Our study unveils a complex pleiotropic genetic landscape in both species, with functional links that relate detox responses and redox mechanisms with developmental and neurobiological processes. These results contribute to deepen our understanding of the role of host plant shifts and natural stress driving ecological specialization.

Soft Selective Sweep on Chemosensory Genes Correlates with Ancestral Preference for Toxic Noni in a Specialist Drosophila Population

Understanding how organisms adapt to environmental changes is a major question in evolution and ecology. In particular, the role of ancestral variation in rapid adaptation remains unclear because its trace on genetic variation, known as soft selective sweep, is often hardly recognizable from genome-wide selection scans. Here, we investigate the evolution of chemosensory genes in Drosophila yakuba mayottensis, a specialist subspecies on toxic noni (Morinda citrifolia) fruits on the island of Mayotte. We combine population genomics analyses and behavioral assays to evaluate the level of divergence in chemosensory genes and perception of noni chemicals between specialist and generalist subspecies of D. yakuba. We identify a signal of soft selective sweep on a handful of genes, with the most diverging ones involving a cluster of gustatory receptors expressed in bitter-sensing neurons. Our results highlight the potential role of ancestral genetic variation in promoting host plant specialization in herbivorous insects and identify a number of candidate genes underlying behavioral adaptation.

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.

Combining Experimental Evolution and Genomics to Understand How Seed Beetles Adapt to a Marginal Host Plant

Genes that affect adaptive traits have been identified, but our knowledge of the genetic basis of adaptation in a more general sense (across multiple traits) remains limited. We combined population-genomic analyses of evolve-and-resequence experiments, genome-wide association mapping of performance traits, and analyses of gene expression to fill this knowledge gap and shed light on the genomics of adaptation to a marginal host (lentil) by the seed beetle Callosobruchus maculatus. Using population-genomic approaches, we detected modest parallelism in allele frequency change across replicate lines during adaptation to lentil. Mapping populations derived from each lentil-adapted line revealed a polygenic basis for two host-specific performance traits (weight and development time), which had low to modest heritabilities. We found less evidence of parallelism in genotype-phenotype associations across these lines than in allele frequency changes during the experiments. Differential gene expression caused by differences in recent evolutionary history exceeded that caused by immediate rearing host. Together, the three genomic datasets suggest that genes affecting traits other than weight and development time are likely to be the main causes of parallel evolution and that detoxification genes (especially cytochrome P450s and beta-glucosidase) could be especially important for colonization of lentil by C. maculatus.