Natural variation in chemosensation: lessons from an island nematode

All organisms must interact with their environment, responding in behavioral, chemical, and other ways to various stimuli throughout their life cycles. Characterizing traits that directly represent an organism's ability to sense and react to their environment provides useful insight into the evolution of life-history strategies. One such trait for the nematode Pristionchus pacificus, chemosensation, is involved in navigation to beetle hosts. Essential for the survival of the nematode, chemosensory behavior may be subject to variation as nematodes discriminate among chemical cues to complete their life cycle. We examine this hypothesis using natural isolates of P. pacificus from La Réunion Island. We select strains from a variety of La Réunion beetle hosts and geographic locations and examine their chemoattraction response toward organic compounds, beetle washes, and live beetles. We find that nematodes show significant differences in their response to various chemicals and are able to chemotax to live beetles in a novel assay. Further, strains can discriminate among different cues, showing more similar responses toward beetle washes than to organic compounds in cluster analyses. However, we find that variance in chemoattraction response is not significantly associated with temperature, location, or beetle host. Rather, strains show a more concerted response toward compounds they most likely directly encounter in the wild. We suggest that divergence in odor-guided behavior in P. pacificus may therefore have an important ecological component.


Introduction
Understanding the ways in which organisms interact with their surroundings requires identification and characterization of the phenotypic traits involved in environmental response. Ecologically relevant traits, such as those involved in environmental sensing and monitoring, are particularly useful, because they represent an organism's ability to identify and cope with environmental change. In nematodes, an important component of the sensing mechanisms mediating interactions with external stimuli is olfactory-based (Bargmann et al. 1993;Jovelin et al. 2003;Hallem et al. 2011). Through odor detection, nematodes have been shown to navigate toward to or away from host excretory products (Grewal et al. 1993), bacterial symbionts and changes in pH (Pye and Burman 1981), temperature (Burman and Pye 1980), electrical field (Shapiro-Ilan et al. 2012), and various plant volatile compounds (Laznik and Trdan 2013).
The hermaphroditic nematode Pristionchus pacificus has an additional environmental interaction that utilizes its olfactory apparatusa necromenic association with host beetles (Herrmann et al. 2006(Herrmann et al. , 2010Hong and Sommer 2006). This association begins as P. pacificus selectively follows remote volatile cues from a potential host. Upon contact and recognition, the nematode passively inhabits the host in an arrested "dauer" stage until the beetles' demise (Weller et al. 2010). Thereafter, P. pacificus detects a suite of chemical substances that cause it to resume development, feeding on the microbes that flourish on the beetle carcass (Herrmann et al. 2006(Herrmann et al. , 2007Mayer and Sommer 2011).
Pristionchus pacificus has been found in association with multiple scarab beetle species throughout a cosmopolitan distribution that encompasses Africa, Asia, Europe, America, and the Mascareigne Islands of the Indian Ocean (Herrmann et al. 2006(Herrmann et al. , 2010. In the latter location, La R eunion Island (2-3 Ma; 2512 km 2 ; Fig. 1) appears to be an oasis for P. pacificus (e.g., Herrmann et al. 2010;Morgan et al. 2012;McGaughran et al. 2013). Considered to be a major biodiversity hotspot (Myers et al. 2000;Thebaud et al. 2009), La R eunion is the largest and most complex island in the Mascareigne group, both topographically and ecologically. With a mix of inland steep relief and coastal lowlands, the island is home to a complex of habitat types or "ecozones", within which species diversity and endemism broadly differ (Strasberg et al. 2005).
In recent years (Herrmann et al. 2010;Morgan et al. 2012), nematode samples have been collected across the R eunion ecozone tableau from several scarab beetle species, both recently introduced (<2000 ya; Vercambre et al. 1991;Cheke and Hume 2008) and endemic. Thus, ecozone (i.e., geological, climatic) and host differences among strains together provide a unique template from which examination of ecological trait differentiation in the context of host choice may be possible.
Chemosensory behavior has been extensively studied in the model nematode Caenorhabditis elegans, which has been shown to follow both water-soluble and volatile cues (e.g., Bargmann et al. 1993). In Pristionchus nematodes, closely related pairs of species show significantly diverged attraction to insect-and plant-related organic compounds in a species-specific manner (Hong and Sommer 2006;Rae et al. 2008), suggesting that the genus has evolved to intercept the chemical communication system of its members' insect hosts (Hong et al. 2008;Sokolowski and Fitzpatrick 2008). Intraspecific chemoattraction differentiation has also been demonstrated among globally distributed P. pacificus strains in their response to beetle pheromones (Hong and Sommer 2006;Hong et al. 2008), but has currently not been studied systematically utilizing the unique La Réunion system.
Here, we use chemoattraction assays to examine whether chemosensory behavior in P. pacificus may be subject to high degrees of variation as nematodes discriminate among chemical cues to complete their life cycle. To achieve this, we examine chemosensation in the context of both host beetle and habitat differences among strains. We demonstrate that La Réunion P. pacificus strains show extreme variation in their odor-guided behavior and that this variance has an ecological component that is driven by concerted responses among strains toward compounds most likely directly encountered in the wild.

Nematode strains
In total, 62 nematode strains were selected for transfer from laboratory culture to nematode growth medium (NGM) agar plates to start isogenic lines for various chemosensation assays. Of these, 61 strains were originally isolated from several locations/scarab beetle hosts across La R eunion Island (Fig. 1, Table 1). One further strain (RS2333) originally isolated from soil in Pasadena (California, USA) serves as a reference strain for our laboratory and was used for quality control in the current study. During the various experimental periods, all strains were maintained, following standard protocols, at 20°C on NGM-agar media and fed with Escherichia coli OP50 strain (Brenner 1974).

Standard assays
Standard chemosensation assays largely followed the procedure outlined in Hong and Sommer (2006), which represents the standard procedure for free-living nematode chemoattraction studies (see Hong and Sommer 2006;and references therein). Briefly, a total of five organic compounds, four beetle wash extracts, and one pure pheromone extract were surveyed for their attractiveness to each of 21 P. pacificus strains (Table 1). These strains were selected because they are closely related genetically (all lineage "C" strains; see Morgan et al. 2012), thus largely removing genetic covariance effects from subsequent analyses. They were also chosen to cover a range of beetle hosts/locations (  were purchased from Sigma-Aldrich (>97% purity; MO) and diluted in 100% ethanol (v/v) to varying concentrations (0.1%, 1%, 10% 100%; Table 2). These compounds/ dilutions were selected based upon their ecological relevance (e.g., as insect allomone/pheromones orbecause host finding can also involve indirect stimuli such as response to plant signals that are released upon attack by beetle larvaeas plant volatiles) and for comparison to the literature protocols (e.g., Hong and Sommer 2006). Further, the dilution range applied reflects the fact that true concentrations in nature remain unknown. For beetle wash extracts, adult beetles of four species collected from La R eunion in January 2010 and 2011 were soaked in dichloromethane for 24 h at room temperature or 23°C. The liquid was then transferred to a new tube, which was either air-or vacuum-dried at 30°C and resuspended in 120 lL of pure ethanol. A solution of pure beetle pheromone extract was purchased from ChemTica Internacional, S.A. This extract was isolated from the Oryctes genus (Scarabaeidae, Coleoptera) and is commonly referred to as "Oryctalure." Its formal designation and chemical formula are ethyl-(RS)-4-methyloctanoate and C 11 H 22 O 2 , respectively. It is commonly used in its pure (undiluted) form in agricultural applications as an insect attractant.
As per Hong and Sommer (2006), each standard chemosensation assay was performed on a plate of mixed-stage nematodes taken from 5-day-old cultures that were started with 10 hermaphrodites grown at 20°C. Before addition to the agar plate, nematodes were washed twice in PBS (+MgCl 2 ) buffer and centrifuged briefly in 15-mL polypropylene tubes. A final concentrated suspension of nematodes was created by removing excess PBS buffer until approximately 100 nematodes were present per 100 lL buffer (as determined over several trials by counting the number of nematodes in 100-lL aliquots following sample preparation). On each plate, 1.5 lL of attractant/counter-attractant and control (ethanol 100% or, in the beetle wash assays, dichloromethane) was spotted at opposite axes in addition to 1.5 lL of 1 mol/L sodium azide, a mild anesthetic to prevent early attracted nematodes from adapting to the odor and leaving the area. Approximately 100 nematodes (see above) were added to the bottom of the plate at a point equidistant from the odor/control sources ( Fig. 2A). Following a 15-h assay period (see Hong and Sommer 2006), plates were first placed briefly into a cold room (8°C) to immobilize nematodes, which were then scored for their location within a 0.5-cm radius circle of each odor source.    A chemotaxis index (CI) was calculated as a ratio for each assay: CI = [No. of animals at attractantÀno. of animals at control]/[Total no. of nematodes at attractant and control] (see Hong and Sommer 2006). Final values for each assay were taken at the same time point (15 h), and the mean CI of all replicates was calculated. Three to five replicates of each assay containing at least ten scorable animals were carried out over the course of the experimental period for every condition tested. In all assays, we define compounds with a mean CI of ≥0.5 as attractive, 0.3-0.5 as weakly attractive, 0.29 to À0.29 as neutral, and ≤0.3 as repulsive (see Hong and Sommer 2006).

Temperature-based chemosensation assays
Given the variety of habitats across La R eunion Island, temperature is likely to be an important factor in environmental response among P. pacificus strains. To examine its potential impact on chemosensation variation, we performed a series of temperature-based chemosensation assays on a subset of six strains (Table 1). These strains were selected to encompass three beetle host/geographic location pairs and again were all lineage "C" strains (Morgan et al. 2012) to remove genetic covariance effects (Table 1). We tested these strains against three compounds: ethyl myristate (1%), Oryctes borbonicus beetle wash, and the pheromone extract, oryctalure ( Table 2). All assay conditions replicated those outlined above, except, following addition of nematodes, assay plates were incubated at four different assay temperatures (15°, 20°, 25°, 30°C).

Multiple-choice assays
Standard chemoattraction assays (above) are capable of demonstrating a preference among strains for the potential attractant relevant to an ethanol control. To test whether strains could discriminate among one attractant relative to another, we performed a series of multiplechoice-based assays using beetle washes (Table 1). Here, the assay conditions replicated those of the standard assay, but five attractants were spotted evenly around the assay plate and nematodes were spotted in the center of the plate equidistant from these odor sources. The CI was then calculated as proportions among attractants relative to each other. A total of fifty assays were performed on strains that were selected to cover ten location/host beetle combinations (n = 5 per combination; Table 1), with each assay replicated three times, to determine whether nematodes would show relative chemoattraction preferences to beetle washes and whether any chemosensation result would relate to the original beetle species of nematode origin.

Live-beetle assays
Odor-guided behavior in the wild most likely occurs in response to blends of chemicals. We and others (e.g., Hong et al. 2008) have used beetle washes to overcome the limitation of chemosensation assays that are based on single compounds. In addition, a limited number of beetles collected from La R eunion Island in January 2011 were able to be returned alive to the Max Planck Institute for Developmental Biology (T€ ubingen, Germany) and maintained at 4°C. This provided an opportunity to develop a novel assay, using live beetles as the attractant source. Thus, several individuals of the beetles Hoplochelus retusa and Hoplia marginalis and a subset of nematode strains (Table 1) were selected in order to determine whether nematodes would chemotax to live beetles. Due to the limited number of live beetles available (n = 36 and 33 for H. retusa and H. marginalis, respectively), and the requirement to use three beetles per assay in H. retusa assays (since this beetle is comparatively small in size), this assay was only available for a select number of strains (four and 11 for H. retusa and H. marginalis assays, respectively). Hence, strains were chosen to cover a variety of beetle hosts and geographic locations (Table 1). These assays followed the standard protocol, with three replicates per assay, except beetles were held in a plastic vestibule that was inserted through the lid of an agar plate so that the tip of the vestibule hovered directly above one axis of the plate. This vestibule was open at the tip so that biological compounds from the live beetles could diffuse down onto the plate, which was sealed at all points of entry. The control in these assays consisted of an empty vestibule adhered to the opposite axis of the agar plate. All other conditions followed those employed for standard assays (above).

Data analysis
All statistical analyses were performed in the programs Microsoft Office Excel ver. 2003 (Microsoft Corp., Redmond, WA) and Minitab ver. 14 (Minitab Inc., State College, PA).
The standard chemoattraction assays were analyzed using summary statistics (e.g., mean, standard deviation) to examine natural variation in chemotaxis index among strains across the various compounds and concentrations tested. An analysis of variance (ANOVA; general linear model, GLM) was used to detect whether response to compounds varied significantly among strains and to examine two potential factors (geographic location and host beetle species) for their effect on chemoattraction variation. Finally, cluster analysis was used to determine how similar the response of strains was toward the various organic compounds versus the beetle washes. Analysis of variance (GLM) was also used to determine whether temperature accounted for variance in chemoattraction. In this analysis strain, temperature and their interaction were tested as factors.
Multiple-choice assays were analyzed using one-way ANOVA (with Tukey's ad hoc pairwise comparisons) to see if there were significant differences among strains in their response when presented with a choice of beetle washes. To examine whether any beetle wash preference in these assays was related to the original beetle host, strain location, or genetic lineage, ANOVA (GLM) was used.
Live-beetle assays were evaluated using summary statistics (mean, standard deviation) and one-way ANOVA to examine whether strains would chemotax to live beetles, and if so, whether the response among strains varied significantly.

Standard assays
Strong levels of natural variation underlie chemosensation in P. pacificus, with the chemoattraction response among strains differing significantly (P < 0.001) with tested compound (Table 3, Fig. 2).
Finally, across compounds, variation was detected (Fig. 2). Generally, most strains were repulsed by benzaldehyde and attracted to ethyl myristate, but strains showed a mix of attraction and repulsion responses to methyl myristate, toluene, and trans-caryophyllene. When presented with beetle washes, strains generally showed repulsion toward Adoretus sp. washes, and a medley of attractive and repulsive strain responses characterized the other washes and the oryctalure pheromone (Fig. 2F).
No dominant signal for preferential attraction among strains toward the beetle wash of the species they were originally isolated from was detected (Fig. 2F). For example, strains originally collected from Adoretus sp. generally showed repulsion or neutral responses to an Adoretus sp. beetle wash (À0.33 AE 0.12, À0.41 AE 0.01, 0.02 AE 0.05, for RS5410, RS5411, and RS5409, respectively). Conversely, strains collected from H. retusa were attracted to H. retusa beetle washes (0.44 AE 0.02 and 0.37 AE 0.08 for RS5405 and RS5350, respectively), while neutral responses characterized H. marginalis-collected strains (0.17 AE 0.06, À0.09 AE 0.03 for RS5427 and RS5429, respectively); responses toward O. borbonicus-collected strains ranged from repulsive (À0.52 AE 0.08 for RS5378) to attractive (0.55 AE 0.03 for RS5399) for the O. borbonicus beetle wash and from À0.30 AE 0.01 (RS5385) to 0.57 AE 0.10 (RS5397) for oryctalure (Fig. 2F). Analysis of variance (GLM) was subsequently used to formally test whether the original beetle host, as well as geographic location, significantly partitioned chemoattrac-tion variance. In this analysis, neither factor significantly predicted CI variance (Table 3). This result was also seen in the cluster analysis, which showed that strain responses toward all compounds were quite dissimilar ( Fig. 3A; Table 4). However, cluster analysis also showed that the response among strains was more similar in assays using beetle washes compared to organic compounds ( Fig. 3B and C). For example, classifying the data into two clusters on the basis of strain CI values resulted in similarity/correlation coefficient distance levels among clusters of 84.037/ 0.319 for beetle wash assays, but 66.591/0.668 for organic compounds (Table 4). When we classified the data into 10 clusters, these values corresponded to 96.186/0.076 and 77.121/0.458 for beetle washes and organic compounds, respectively (Table 4); thus, similarity among strain responses was some 20% higher toward beetle washes than toward organic compounds (Fig. 3).

Temperature-based chemosensation assays
The finding of natural variation at constant temperature (above) was reiterated in our temperature-based assays. In Figure 4, this is demonstrated by the large degree of scatter among sample responses generally, which varied within compounds, across temperatures, and across compounds.  Fig. 4). Response of a single strain across temperatures also varied. For example, strain RS5397 (TB/ O. borbonicus) responses were 0.41 AE 0.24, À0.73 AE 0.36, 0.16 AE 0.30, and À0.64 AE 0.54 for 15°, 20°, 25°, and 30°C, respectively (Fig. 4). However, ANOVA (GLM) analysis found that differences in temperature among strains did not account for differences in chemoattraction response in a global analysis (Table 5). This analysis also found that the interaction between strain and temperature did not significantly partition CI variance (Table 5).

Multiple-choice assays
When presented with a choice of five different beetle washes (Table 1), P. pacificus strains varied significantly in their chemoattraction response (F 4,245 = 4.43; P = 0.002; Table 6). Tukey's pairwise analysis identified that these differences were due to mean strain responses to H. retusa beetle washes from TB and San Souci (SS), which were significantly different from both each other, and from 20% (the null proportion since five choices were presented; Fig. 5). Further analysis of the chemoattraction variance in a GLM framework showed that none of the analyzed factors (beetle host, genetic lineage, and location) significantly partitioned variance in CI for multiple-choice assays. Table 4. Results of cluster analysis examining the similarity among 21 strains of Pristionchus pacificus in their chemoattraction responses toward: all compounds, organic compounds only, and beetle washes only. Step No. of clusters   Figure 4. Results of the temperature-based chemoattraction assays for six Pristionchus pacificus strains. Results are presented for mean AE standard deviation chemoattraction to the compounds ethyl myristate (1%), Oryctes borbonicus beetle wash, and "Oryctalure" pheromone, at four measurement temperatures. Nematode origin information (original beetle host association and geographic location) is given to the right of the strain key for each isolate (Ory, Oryctes borbonicus; Ado, Adoretus sp.; Mal, Maladera affinis; TB, Trois Bassin; GE, Grand Etang; SB, Saint Benoit). All strains are of genetic lineage, "C" (Morgan et al. 2012). See Tables 1, 5 for further information.

Live-beetle assays
In live-beetle assays, strains showed they were able to chemotax toward the signal produced by live beetles.  (Fig. 6).

Discussion
An understanding of organism-environment interactions requires, as a first step, characterization of the potential traits involved in an organisms' sensation of its surroundings. Here, we show that one trait that is particularly important for nematode response to environmental stimuli, chemosensation, exhibits strong levels of variation among strains. This variation was present in P. pacificus across all of our assays, with strains varying in response in a given assay over a range covering up to 85% of the possible chemosensation index. As well as strains showing different responses to a given compound, variation was present across concentrations of the same compound and across compounds themselves.
Cluster analysis of the standard assay data showed that responses were indeed dissimilar among strains, and further statistical analysis indicated that temperature, beetle host, and sampling location were not responsible for chemoattraction variance. However, cluster analysis also showed that responses among strains were more similar when the test compound was a beetle wash than when it was an organic compound. Given that these nematodes show a more concerted response toward the compounds they most likely directly encounter in the wild, this suggests that variation in the odor-guided response of P. pacificus may have an important ecological component.
Our mix of standard, temperature-based, multiplechoice and live-beetle assays showed that P. pacificus is able to discriminate between test compounds. For example, in standard assays, strains showed a general pattern of repulsion toward benzaldehyde and attraction to ethyl myristate. Further differences existed in comparisons of response of organic versus beetle wash assays, and our live-beetle assays demonstrated that strains are capable of responding to volatile cues released by a live beetle. Our multiple-choice assays also indicated that different nematode strains could show relative chemoattraction preferences to beetle washes. The ability of P. pacificus to discriminate among various chemical cues in the current study further supports the ecological importance of this trait. A related nematode, Pristionchus maupasi, can use such discrimination to specifically recognize older feeding beetles, thereby reducing its dauer period to the shorter-lived adult host stage (Hong et al. 2008). Whether P. pacificus is also capable of discrimination across beetle developmental stages is an interesting consideration for future work.
Discrimination across beetle species in the form of preferential host choices could have several important effects  Rae et al. 2008), then, on certain beetle host species, a greater number of generations could occur before nematodes must arrest development, thus possibly influencing strain and/or population survival. Given these potential host choice effects, the lack of a dominant signal for preferential attraction among strains to their original beetle host in our multiple-choice assays is somewhat surprising. However, cluster analysis did show a concerted response among strains toward beetle washes compared to organic compounds, and it may be that our dataset lacked the power (n = 5 strains per host/ location combination and n = 3 replicates per assay) to  Figure 5. Results of multiple-choice assays: (A) Graphic to show the set-up of the assays, where five beetle wash attractants were spotted evenly around the assay plate at points A-E and nematodes were spotted in the center of the plate equidistant from these odor sources. The chemoattraction index (CI) was calculated as proportions among attractants relative to each other. A total of fifty assays using ten location/host beetle combinations (Table 1) were performed, with each assay replicated three times, to determine whether nematodes would show relative chemoattraction preferences to beetle washes; (B) Relative percentage of nematodes that chemotaxed to each point (A-E) on the assay plate for each assay for each beetle wash. Beetle wash is indicated by the key at the bottom of the figure (SB_Mal, Saint Benoit Maladera affinis beetle wash; TB_Hop, Trois Bassin Hoplia retusa beetle wash; SS_Hop, San Souci H. retusa beetle wash; TK_Ado, Takamaka Adoretus sp. beetle wash; TB_Ory, Trois Bassin Oryctes borbonicus beetle wash). Differences in attraction in response to the TB_Hop* and SS_Hop* washes were significant, as indicated by the "*" (Table 6). Nematode origin information (original beetle host association and geographic location) is given below the x-axis for each strain (Ado, Adoretus sp.; Ory, Oryctes borbonicus; Hop, Hoplia retusa; Mal, Maladera affinis; BV, Basse Vall ee; GE, Grand Etang; PA, Palm des Palmiste; PC, Plan des Cafr es; SB, Saint Benoit; SS, San Souci; TB, Trois Bassin; TK, Takamaka). All strains are of genetic lineage, "A", "C", and "D" (Morgan et al. 2012), as indicated in Table 1 The variability of chemoattraction in P. pacificus suggests that response to different volatile cues in this freeliving nematode, aside from showing similarity among beetle washes, may be more of a strain-specific characteristic than one influenced by location or temperature. Strain-specific responses independent of temperature have also been found previously in both entomopathogenic and parasitic nematodes (Ali et al. 2010(Ali et al. , 2011, suggesting that plasticity both within and among populations and/or hosts may be an important feature of the chemoattractive response for some nematodes. Conversely, chemosensation in other species has been shown to be an isolating force among populations. For example, chemosensation has been shown to lead to premating isolation among populations and species of moth (Smadja and Butlin 2009; and references therein) and to behavioral isolation in flies of the genus Drosophila (Etges and Ahrens 2001;Ferveur 2005;Smadja and Butlin 2009). In C. elegans, phenotypic variation in chemosensation has been suggested to be closely related to fitness (Schulenburg and M€ uller 2004;Barri ere and F elix 2005) and therefore potentially under selection.
In P. pacificus, a large degree of variability within both individuals and populations characterizes chemosensation, as shown here and in Hong et al. (2008), but also a number of other traits (e.g., Bento et al. 2010;Mayer and Sommer 2011). Across species, much experimental data suggest that individual differences in behavioral strategies employed at the intrapopulation level are important in optimizing fitness (e.g. , Ehlinger 1990;Paffenh€ ofer 1994;Utne et al. 1997;Sneddon 2003;Van Oers et al. 2005). Thus, the interindividual variability identified here may suggest that chemosensation variation plays a role in population persistence for P. pacificus by allowing strains to better match their phenotype to their environment (Viney and Diaz 2012).   TB  TBG  SB  TB  TB  TBG  SB  TB  TBG  TB  TB  TB  TB  TB  TBG   Live beetles  Control Plate with agar Nematodes (A)

(B)
Hoplochelus beetle wash Hoplia beetle wash Figure 6. Results of live-beetle assays for Pristionchus pacificus; (A) Schematic of the experimental set-up; (B) Variation in odor-response profile of selected P. pacificus strains in live-beetle assays using the beetles Hoplia retusa (left four strains) and Hoplochelus marginalis (right 11 strains) from La R eunion Island. Results are presented as follows: mean chemotaxis index AE standard deviation. Nematode origin information (original beetle host association and geographic location) is given below the x-axis for each strain (Hop, Hoplia retusa; Hch, Hoplochelus marginalis; Mal, Maladera affinis; Ory, Oryctes borbonicus; TB, Trois Bassin; TBG, Trois Bassin Garden; SB, Saint Benoit). All strains are of genetic lineage, "A" or "C" (Morgan et al. 2012). See Table 1 for further information. We have demonstrated a substantial degree of natural variation in chemosensation for P. pacificus and highlighted the ability of this nematode to discriminate among a variety of chemical cues. We suggest that divergence in odor-guided behavior in P. pacificus may have an important ecological component, whereby strains show a concerted response toward the compounds they most likely directly encounter in the wild. Our results provide a comprehensive basis for the evolution and ecology of P. pacificus chemosensation. It would be useful to investigate whether other ecological traits of the nematodes of La R eunion are diverging, and future work should focus on elucidation of the genetic basis and evolutionary mechanisms underlying such trait divergence (e.g., Hong et al. 2008;Bento et al. 2010;Mayer and Sommer 2011). With such studies, the future of this system promises success in bridging the boundaries of ecological context, phenotypic expression, gene function, and molecular variation.