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

Sexual Selection May Not Often Reduce Gene Flow Between Locally Adapted Populations. A Review of Some Evidence, and Suggestions for Better Tests

Sexually selected traits often depend on an individual’s physical condition, or otherwise indirectly reflect the ecological performance of individuals. When individuals disperse between populations that are locally adapted to different environments, their ecological performance may decline. This in turn may result in more poorly expressed sexual traits, and therefore in a lower reproductive success. Hence, sexual selection may reduce the effective gene flow between populations, and thereby maintain or even enhance population divergence. This hypothesis was published in a highly visible journal (van Doorn et al., 2009, Science). Here I review the subsequently published empirical tests of this hypothesis. I downloaded all metadata (incl. abstracts) of papers citing van Doorn et al. (2009) and read those papers that undertook relevant tests. To my surprise, only very few papers provided explicit tests of the hypothesis, this never involved plants, and only one study found support for it. While sexual selection may therefore not often reduce gene flow between locally adapted populations, some improvements to experimental design and choice of study system are noted. I therefore also provide a detailed list of suggestions for high quality tests of this hypothesis. This hopefully acts as a catalyst for more and better studies to test whether sexual and natural selection can work in synergy to reduce effective dispersal, and thereby protect and promote adaptive population divergence.

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.

Evolutionary genomics of host plant adaptation: insights from Drosophila

Variation in gene expression in response to the use of alternate host plants can reveal genetic and physiological mechanisms explaining why insect-host relationships vary from host specialism to generalism. Interpreting transcriptome variation relies on well-annotated genomes, making drosophilids valuable model systems, particularly those species with tractable ecological associations. Patterns of whole genome expression and alternate gene splicing in response to growth on different hosts have revealed expression of gene networks of known detoxification genes as well as novel functionally enriched genes of diverse metabolic and structural functions. Integrating trancriptomic responses with fitness differences and levels of phenotypic plasticity in response to alternate hosts will help to reveal the general nature of genotype-phenotype relationships.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Identification and Expression Profiling of Peripheral Olfactory Genes in the Parasitoid Wasp Aphidius ervi (Hymenoptera: Braconidae) Reared on Different Aphid Hosts

Generalist parasitoids of aphids, such as the wasp Aphidius ervi, display significant differences in terms of host preference and host acceptance, depending on the host on which they developed (natal host), which is preferred over a non-natal host, a trait known as host fidelity. This trait allows females to quickly find hosts in heterogeneous environments, a process mediated by chemosensory/olfactory mechanisms, as parasitoids rely on olfaction and chemical cues during host selection. Thus, it is expected that proteins participating in chemosensory recognition, such as odorant-binding proteins (OBPs) and odorant receptors (ORs) would play a key role in host preference. In this study, we addressed the effect of parasitoid reciprocal host switching between two aphid hosts (Sitobion avenae and Acyrthosiphon pisum) on the expression patterns of chemosensory genes in the wasp A. ervi. First, by using a transcriptomic approach based on RNAseq of A. ervi females reared on S. avenae and A. pisum, we were able to annotate a total of 91 transcripts related to chemoperception. We also performed an in-silico expression analysis and found three OBPs and five ORs displaying different expression levels. Then, by using qRT-PCR amplification, we found significant differences in the expression levels of these eight genes when the parasitoids were reciprocally transplanted from S. avenae onto A. pisum and vice versa. This suggests that the expression levels of genes coding for odorant receptors and odorant-binding proteins would be regulated by the specific plant–aphid host complex where the parasitoids develop (maternal previous experience) and that chemosensory genes coding for olfactory mechanisms would play a crucial role on host preference and host acceptance, ultimately leading to the establishment of host fidelity in A. ervi parasitoids.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Mating Changes Sexually Dimorphic Gene Expression in the Seed Beetle Callosobruchus maculatus

Sexually dimorphic phenotypes arise largely from sex-specific gene expression, which has mainly been characterized in sexually naïve adults. However, we expect sexual dimorphism in transcription to be dynamic and dependent on factors such as reproductive status. Mating induces many behavioral and physiological changes distinct to each sex and is therefore expected to activate regulatory changes in many sex-biased genes. Here, we first characterized sexual dimorphism in gene expression in Callosobruchus maculatus seed beetles. We then examined how females and males respond to mating and how it affects sex-biased expression, both in sex-limited (abdomen) and sex-shared (head and thorax) tissues. Mating responses were largely sex-specific and, as expected, females showed more genes responding compared with males (∼2,000 vs. ∼300 genes in the abdomen, ∼500 vs. ∼400 in the head and thorax, respectively). Of the sex-biased genes present in virgins, 16% (1,041 genes) in the abdomen and 17% (243 genes) in the head and thorax altered their relative expression between the sexes as a result of mating. Sex-bias status changed in 2% of the genes in the abdomen and 4% in the head and thorax following mating. Mating responses involved de-feminization of females and, to a lesser extent, de-masculinization of males relative to their virgin state: mating decreased rather than increased dimorphic expression of sex-biased genes. The fact that regulatory changes of both types of sex-biased genes occurred in both sexes suggests that male- and female-specific selection is not restricted to male- and female-biased genes, respectively, as is sometimes assumed.

Transcriptome modulation during host shift is driven by secondary metabolites in desert Drosophila

High-throughput transcriptome studies are breaking new ground to investigate the responses that organisms deploy in alternative environments. Nevertheless, much remains to be understood about the genetic basis of host plant adaptation. Here, we investigate genome-wide expression in the fly Drosophila buzzatii raised in different conditions. This species uses decaying tissues of cactus of the genus Opuntia as primary rearing substrate and secondarily, the necrotic tissues of the columnar cactus Trichocereus terscheckii. The latter constitutes a harmful host, rich in mescaline and other related phenylethylamine alkaloids. We assessed the transcriptomic responses of larvae reared in Opuntia sulphurea and T. terscheckii, with and without the addition of alkaloids extracted from the latter. Whole-genome expression profiles were massively modulated by the rearing environment, mainly by the presence of T. terscheckii alkaloids. Differentially expressed genes were mainly related to detoxification, oxidation-reduction and stress response; however, we also found genes involved in development and neurobiological processes. In conclusion, our study contributes new data onto the role of transcriptional plasticity in response to alternative rearing environments.

The Interplay of Temperature and Genotype on Patterns of Alternative Splicing in Drosophila melanogaster

Alternative splicing is the highly regulated process of variation in the removal of introns from premessenger-RNA transcripts. The consequences of alternative splicing on the phenotype are well documented, but the impact of the environment on alternative splicing is not yet clear. We studied variation in alternative splicing among four different temperatures, 13, 18, 23, and 29°, in two Drosophila melanogaster genotypes. We show plasticity of alternative splicing with up to 10% of the expressed genes being differentially spliced between the most extreme temperatures for a given genotype. Comparing the two genotypes at different temperatures, we found <1% of the genes being differentially spliced at 18°. At extreme temperatures, however, we detected substantial differences in alternative splicing-with almost 10% of the genes having differential splicing between the genotypes: a magnitude similar to between species differences. Genes with differential alternative splicing between genotypes frequently exhibit dominant inheritance. Remarkably, the pattern of surplus of differences in alternative splicing at extreme temperatures resembled the pattern seen for gene expression intensity. Since different sets of genes were involved for the two phenotypes, we propose that purifying selection results in the reduction of differences at benign temperatures. Relaxed purifying selection at temperature extremes, on the other hand, may cause the divergence in gene expression and alternative splicing between the two strains in rarely encountered environments.

Buffering of Genetic Regulatory Networks in Drosophila melanogaster

Regulatory variation in gene expression can be described by cis- and trans-genetic components. Here we used RNA-seq data from a population panel of Drosophila melanogaster test crosses to compare allelic imbalance (AI) in female head tissue between mated and virgin flies, an environmental change known to affect transcription. Indeed, 3048 exons (1610 genes) are differentially expressed in this study. A Bayesian model for AI, with an intersection test, controls type I error. There are ∼200 genes with AI exclusively in mated or virgin flies, indicating an environmental component of expression regulation. On average 34% of genes within a cross and 54% of all genes show evidence for genetic regulation of transcription. Nearly all differentially regulated genes are affected in cis, with an average of 63% of expression variation explained by the cis-effects. Trans-effects explain 8% of the variance in AI on average and the interaction between cis and trans explains an average of 11% of the total variance in AI. In both environments cis- and trans-effects are compensatory in their overall effect, with a negative association between cis- and trans-effects in 85% of the exons examined. We hypothesize that the gene expression level perturbed by cis-regulatory mutations is compensated through trans-regulatory mechanisms, e.g., trans and cis by trans-factors buffering cis-mutations. In addition, when AI is detected in both environments, cis-mated, cis-virgin, and trans-mated-trans-virgin estimates are highly concordant with 99% of all exons positively correlated with a median correlation of 0.83 for cis and 0.95 for trans We conclude that the gene regulatory networks (GRNs) are robust and that trans-buffering explains robustness.

Preadult life history variation determines adult transcriptome expression

Preadult determinants of adult fitness and behaviour have been documented in a variety of organisms with complex life cycles, but little is known about expression patterns of genes underlying these adult traits. We explored the effects of differences in egg-to-adult development time on adult transcriptome and cuticular hydrocarbon variation in order to understand the nature of the genetic correlation between preadult development time and premating isolation between populations of Drosophila mojavensis reared in different host cactus environments. Transcriptome variation was analysed separately in flies reared on each host and revealed that hundreds of genes in adults were differentially expressed (FDR P < 0.05) due to development time differences. For flies reared on pitaya agria cactus, longer preadult development times caused increased expression of genes in adults enriched for ribosome production, protein metabolism, chromatin remodelling and regulation of alternate splicing and transcription. Baja California flies reared on organ pipe cactus showed fewer differentially expressed genes in adults due to longer preadult development time, but these were enriched for ATP synthesis and the TCA cycle. Mainland flies reared on organ pipe cactus with shorter development times showed increased transcription of genes enriched for mitochondria and energy production, protein synthesis and glucose metabolism: adults with longer development times had increased expression of genes enriched for adult life span, cuticle proteins and ion binding, although most differentially expressed genes were unannotated. Differences due to population, sex, mating status and their interactions were also assessed. Adult cuticular hydrocarbon profiles also showed shifts due to egg-to-adult development time and were influenced by population and mating status. These results help to explain why preadult life history variation determines subsequent expression of the adult transcriptome along with traits involved with reproductive isolation and revealed previously undocumented connections between genetic and environmental influences over the entire life cycle in this desert insect.

Gene Expression Reaction Norms Unravel the Molecular and Cellular Processes Underpinning the Plastic Phenotypes of Alternanthera Philoxeroides in Contrasting Hydrological Conditions

Alternanthera philoxeroides is an amphibious invasive weed that can colonize both aquatic and terrestrial habitats. Individuals growing in different habitats exhibit extensive phenotypic variation but little genetic differentiation. Little is known about the molecular basis underlying environment-induced phenotypic changes. Variation in transcript abundance in A. philoxeroides was characterized throughout the time-courses of pond and upland treatments using RNA-Sequencing. Seven thousand eight hundred and five genes demonstrated variable expression in response to different treatments, forming 11 transcriptionally coordinated gene groups. Functional enrichment analysis of plastically expressed genes revealed pathway changes in hormone-mediated signaling, osmotic adjustment, cell wall remodeling, and programmed cell death, providing a mechanistic understanding of the biological processes underlying the phenotypic changes in A. philoxeroides. Both transcriptional modulation of environmentally sensitive loci and environmentally dependent control of regulatory loci influenced the plastic responses to the environment. Phenotypic responses and gene expression patterns to contrasting hydrological conditions were compared between A. philoxeroides and its alien congener Alternanthera pungens. The terricolous A. pungens displayed limited phenotypic plasticity to different treatments. It was postulated based on gene expression comparison that the interspecific variation in plasticity between A. philoxeroides and A. pungens was not due to environmentally-mediated changes in hormone levels but to variations in the type and relative abundance of different signal transducers and receptors expressed in the target tissue.

Differential expression of the chemosensory transcriptome in two populations of the stemborer Sesamia nonagrioides

Among the proposed mechanisms of local adaptation to different ecological environments, transcriptional changes may play an important role. In this study, we investigated whether such variability occurred within the chemosensory organs of a herbivorous insect, for which chemosensation guides most of its host preferences. A European and an African population of the noctuid Sesamia nonagrioides that display significant differences in their ecological preferences were collected on Zea mays and Typha domingensis, respectively. RNAseq were used between the two populations for digital expression profiling of chemosensory organs from larval antennae and palps. Preliminary data on adult female antennae and ovipositors were also collected. We found 6,550 differentially expressed transcripts in larval antennae and palps. Gene ontology enrichment analyses suggested that transcriptional activity was overrepresented in the French population and that virus and defense activities were overrepresented in the Kenyan population. In addition, we found differential expression of a variety of cytochrome P450s, which may be linked to the different host-plant diets. Looking at olfactory genes, we observed differential expression of numerous candidate odorant-binding proteins, chemosensory proteins, and one olfactory receptor, suggesting that differences in olfactory sensitivity participate in insect adaptation.

Transcriptome-Wide Differential Gene Expression in Bicyclus anynana Butterflies: Female Vision-Related Genes Are More Plastic

Vision is energetically costly to maintain. Consequently, over time many cave-adapted species downregulate the expression of vision genes or even lose their eyes and associated eye genes entirely. Alternatively, organisms that live in fluctuating environments, with different requirements for vision at different times, may evolve phenotypic plasticity for expression of vision genes. Here, we use a global transcriptomic and candidate gene approach to compare gene expression in the heads of a polyphenic butterfly. Bicyclus anynana have two seasonal forms that display sexual dimorphism and plasticity in eye morphology, and female-specific plasticity in opsin gene expression. Nonchoosy dry season females downregulate opsin expression, consistent with the high physiological cost of vision. To identify other genes associated with sexually dimorphic and seasonally plastic differences in vision, we analyzed RNA-sequencing data from whole head tissues. We identified two eye development genes (klarsicht and warts homologs) and an eye pigment biosynthesis gene (henna) differentially expressed between seasonal forms. By comparing sex-specific expression across seasonal forms, we found that klarsicht, warts, henna, and another eye development gene (domeless) were plastic in a female-specific manner. In a male-only analysis, white (w) was differentially expressed between seasonal forms. Reverse transcription polymerase chain reaction confirmed that warts and white are expressed in eyes only, whereas klarsicht, henna and domeless are expressed in both eyes and brain. We find that differential expression of eye development and eye pigment genes is associated with divergent eye phenotypes in B. anynana seasonal forms, and that there is a larger effect of season on female vision-related genes.

Genome-wide DNA methylation patterns in wild samples of two morphotypes of threespine stickleback (Gasterosteus aculeatus)

Epigenetic marks such as DNA methylation play important biological roles in gene expression regulation and cellular differentiation during development. To examine whether DNA methylation patterns are potentially associated with naturally occurring phenotypic differences, we examined genome-wide DNA methylation within Gasterosteus aculeatus, using reduced representation bisulfite sequencing. First, we identified highly methylated regions of the stickleback genome, finding such regions to be located predominantly within genes, and associated with genes functioning in metabolism and biosynthetic processes, cell adhesion, signaling pathways, and blood vessel development. Next, we identified putative differentially methylated regions (DMRs) of the genome between complete and low lateral plate morphs of G. aculeatus. We detected 77 DMRs that were mainly located in intergenic regions. Annotations of genes associated with these DMRs revealed potential functions in a number of known divergent adaptive phenotypes between G. aculeatus ecotypes, including cardiovascular development, growth, and neuromuscular development.

Transcriptional and epigenetic responses to mating and aging in Drosophila melanogaster

BACKGROUND

Phenotypic plasticity allows organisms to respond rapidly to changing environmental circumstances, and understanding its genomic basis can yield insights regarding the underlying genes and genetic networks affecting complex phenotypes. Female Drosophila melanogaster undergo dramatic physiological changes mediated by seminal fluid components transferred upon mating, including decreased longevity. Their physiological and behavioral effects have been well characterized, but little is known about resulting changes in regulation of gene expression or the extent to which mating-induced changes in gene expression are the same as those occurring during aging.

RESULTS

We assessed genome-wide mRNA, microRNA, and three common histone modifications implicated in gene activation for young and aged virgin and mated female D. melanogaster in a factorial design. We identified phenotypically plastic transcripts and epigenetic modifications associated with mating and aging. We used these data to derive phenotypically plastic regulatory networks associated with mating of young flies, and aging of virgin and mated flies. Many of the mRNAs, microRNAs and epigenetic modifications associated with mating of young flies also occur with age in virgin flies, which may reflect mating-induced accelerated aging. We functionally tested the plastic regulatory networks by overexpressing environmentally sensitive microRNAs. Overexpression resulted in altered expression of ~70% of candidate target genes, and in all cases affected oviposition.

CONCLUSIONS

Our results implicate microRNAs as mediators of phenotypic plasticity associated with mating and provide a comprehensive documentation of the genomic and epigenomic changes that accompany mating- and aging-induced physiological changes in female D. melanogaster.

Identification of differentially expressed genes in female Drosophila antonietae and Drosophila meridionalis in response to host cactus odor

Background

Studies of insect-plant interactions have provided critical insights into the ecology and evolution of adaptive processes within and among species. Cactophilic Drosophila species have received much attention because larval development occurs in the necrotic tissues of cacti, and both larvae and adults feed on these tissues. Such Drosophila-cactus interactions include effects of the host plant on the physiology and behavior of the flies, especially so their nutritional status, mating condition and reproduction. The aim of this work was to compare the transcriptional responses of two species, Drosophila antonietae and Drosophila meridionalis, and identify genes potentially related to responses to odors released by their host cactus, Cereus hildmannianus. The two fly species are sympatric in most of their populations and use this same host cactus in nature.

Results

We obtained 47 unique sequences (USs) for D. antonietae in a suppression subtractive hybridization screen, 30 of these USs had matches with genes predicted for other Drosophila species. For D. meridionalis we obtained 81 USs, 46 of which were orthologous with genes from other Drosophila species. Functional information (Gene Ontology) revealed that these differentially expressed genes are related to metabolic processes, detoxification mechanisms, signaling, response to stimuli, and reproduction. The expression of 13 genes from D. meridionalis and 12 from D. antonietae were further analyzed by quantitative real time-PCR, showing that four genes were significantly overexpressed in D. antonietae and six in D. meridionalis.

Conclusions

Our results revealed the differential expression of genes related to responses to odor stimuli by a cactus, in two associated fly species. Although the majority of activated genes were similar between the two species, we also observed that certain metabolic pathways were specifically activated, especially those related to signaling pathways and detoxification mechanisms. The activation of these genes may reflect different metabolic pathways used by these flies in their interaction with this host cactus. Our findings provide insight into how the use of C. hildmannianus may have arisen independently in the two fly species, through genetic differentiation in metabolic pathways to effectively explore this cactus as a host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0191-2) 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.

Sequence evolution and expression regulation of stress-responsive genes in natural populations of wild tomato

The wild tomato species Solanum chilense and S. peruvianum are a valuable non-model system for studying plant adaptation since they grow in diverse environments facing many abiotic constraints. Here we investigate the sequence evolution of regulatory regions of drought and cold responsive genes and their expression regulation. The coding regions of these genes were previously shown to exhibit signatures of positive selection. Expression profiles and sequence evolution of regulatory regions of members of the Asr (ABA/water stress/ripening induced) gene family and the dehydrin gene pLC30-15 were analyzed in wild tomato populations from contrasting environments. For S. chilense, we found that Asr4 and pLC30-15 appear to respond much faster to drought conditions in accessions from very dry environments than accessions from more mesic locations. Sequence analysis suggests that the promoter of Asr2 and the downstream region of pLC30-15 are under positive selection in some local populations of S. chilense. By investigating gene expression differences at the population level we provide further support of our previous conclusions that Asr2, Asr4, and pLC30-15 are promising candidates for functional studies of adaptation. Our analysis also demonstrates the power of the candidate gene approach in evolutionary biology research and highlights the importance of wild Solanum species as a genetic resource for their cultivated relatives.