Oviposition traits generate extrinsic postzygotic isolation between two pine sawfly species
© The Author(s). 2017
Received: 22 September 2016
Accepted: 5 January 2017
Published: 19 January 2017
Although empirical data indicate that ecological speciation is prevalent in nature, the relative importance of different forms of reproductive isolation and the traits generating reproductive isolation remain unclear. To address these questions, we examined a pair of ecologically divergent pine-sawfly species: while Neodiprion pinetum specializes on a thin-needled pine (Pinus strobus), N. lecontei utilizes thicker-needled pines. We hypothesized that extrinsic postzygotic isolation is generated by oviposition traits. To test this hypothesis, we assayed ovipositor morphology, oviposition behavior, and host-dependent oviposition success in both species and in F1 and backcross females.
Compared to N. lecontei, N. pinetum females preferred P. strobus more strongly, had smaller ovipositors, and laid fewer eggs per needle. Additionally, we observed host- and trait-dependent reductions in oviposition success in F1 and backcross females. Hybrid females that had pinetum-like host preference (P. strobus) and lecontei-like oviposition traits (morphology and egg pattern) fared especially poorly.
Together, these data indicate that maladaptive combinations of oviposition traits in hybrids contribute to extrinsic postzygotic isolation between N. lecontei and N. pinetum, suggesting that oviposition traits may be an important driver of divergence in phytophagous insects.
KeywordsEcological speciation Reproductive barriers Host adaptation Phytophagous insect Diprionidae Neodiprion
Evolutionary biologists have long recognized that natural selection plays an important role in the formation of new species [1–4]. However, it is only within the last two decades that ecological speciation—the process by which environmentally based divergent selection gives rise to reproductive isolation [5, 6]—has become the focus of sustained research effort. During this time, laboratory and field studies in a wide range of organisms have demonstrated unequivocally that ecological speciation occurs in nature [7–10]. Moreover, comparative data suggest that ecological divergence plays a fundamental and taxonomically general role in driving speciation . Nevertheless, while some aspects of ecological speciation are now fairly well understood, many major questions—including the relative importance of different forms of reproductive isolation (RI), and the types of traits that generate RI—remain unresolved .
Any form of RI can, in theory, contribute to ecological speciation so long as it arises as a consequence of divergent natural selection. However, one form of RI that may be especially important is extrinsic postzygotic isolation (hereafter, EPI), in which intermediacy or maladaptive combinations of traits in hybrids causes low fitness in both parental environments [7, 12–14]. EPI is thought to be a particularly common form of RI in ecological speciation because, so long as hybrids are intermediate and intermediate environments are lacking, it is a direct and automatic result of divergent selection . As such, EPI should be among the earliest barriers to arise during ecological speciation [8, 15]. Additionally, when there is gene flow between diverging populations, EPI will lead to direct selection for assortative mating via reinforcement [4, 16]. However, although EPI is one of only two forms of RI that are unique to ecological speciation (the other being immigrant inviability, ), it is understudied relative to other forms of RI. One possible reason for the dearth of EPI studies is that it is challenging to distinguish between extrinsic (ecologically dependent) and intrinsic (due to genetic incompatibilities) sources of reduced hybrid fitness .
To date, three techniques have been proposed to distinguish between extrinsic and intrinsic sources of postzygotic isolation. The simplest of these is to compare the fitness of F1 hybrids in the wild to their fitness in a benign environment, in which the source of ecologically based selection has presumably been removed. If reduced hybrid fitness disappears in the “benign” habitat, this implies that the reduction was environmentally dependent . The main limitation of this approach is that it does not control for stress-related expression of intrinsic hybrid incompatibilities . A second, more rigorous approach is to rear backcrosses of F1s to both parental forms in both parental environments. EPI predicts that each backcross type will perform best in the parental habitat to which it is most genetically similar . A final technique is to examine how specific hybrid traits impact fitness in parental habitats. This approach requires knowledge of the traits contributing to EPI and can be accomplished in one of two ways, the first of which is to experimentally manipulate parental individuals to resemble hybrids (modify-parental-phenotype) [7, 8]. Parents with the hybrid trait are expected to have reduced fitness in the parental environments if that trait is generating EPI. For many traits and organisms, however, these phenotypic modifications would be impractical, if not impossible. An alternative to direct modification is to take advantage of trait variation in F1 hybrids, F2, or backcross individuals and track how different trait values and combinations impact fitness in parental environments (e.g, [19, 20]).
One group of organisms that has featured prominently in empirical and theoretical studies of ecological speciation and EPI is plant-feeding insects. Several lines of evidence support the hypothesis that changes in host use are an important driver of ecological speciation in insects, including: (1) phylogenetic studies that show elevated rates of diversification among lineages of phytophagous insects compared to non-phytophagous insects [21–23], (2) comparative studies that demonstrate an association between changes in host use and speciation [11, 24, 25] (but see ), and (3) a growing list of empirical case studies that have confirmed key predictions of ecological speciation (reviewed in [8, 27–29]). However, while evidence supporting ecological speciation in insects is strong, the contribution of EPI remains unknown. In particular, although indirect evidence for host-related EPI exists for many taxa (reviewed in ), few direct tests exist (but see [20, 30–33]). Moreover, in most cases, the specific traits contributing to EPI have not been identified (but see ). Understanding the mechanistic basis of EPI is critical if we are to understand whether biases exist in the types of traits (e.g., morphological, physiological, behavioral) that contribute to reduced hybrid performance.
To investigate EPI—and the traits that produce it—we focus here on pine sawflies in the genus Neodiprion (Hymenoptera, Diprionidae), a Holarctic group of pine specialists that develop in intimate association with their host plants: adults mate on the host, females embed their eggs in the host tissue, and larvae complete their development on the host, spinning their cocoons on or beneath the host . Population genomic data from a single species, N. lecontei, indicate that divergence in host use contributes to population differentiation , and comparative data from multiple species indicate that host-associated population differentiation occasionally progresses to speciation . However, the mechanisms linking divergent host use to population differentiation and RI have not been identified.
To explore mechanistic links between host-use divergence and speciation in Neodiprion, we examined a pair of sister species that differ in host use, N. pinetum and N. lecontei [36, 37]. N. pinetum is a specialist that feeds on Pinus strobus, while N. lecontei feeds on a wider range of Pinus hosts, but generally avoids P. strobus. These species will mate under no-choice conditions in the lab and produce viable, fertile offspring (personal observation). Two lines of evidence indicate that they hybridize in the wild as well: (1) we have collected hybrids—which are identifiable via their intermediate larval coloration—at multiple field sites (personal observation), and (2) mitochondrial introgression has occurred between these two species . Nevertheless, despite widespread sympatry and occasional hybridization, these two species remain morphologically, behaviorally, and genetically distinct. These observations suggest that there are postzygotic barriers to gene flow. Given that lab-reared hybrids are viable and fertile (personal observation), we hypothesize that postzygotic barriers between N. lecontei and N. pinetum are largely extrinsic in nature, stemming from their specialization on different Pinus hosts.
There are several oviposition traits that vary among Neodiprion species that may be shaped by divergent natural selection stemming from host plant needle characteristics. First, there is interspecific variation in the number of eggs laid per needle: while some species lay a single egg per needle, others lay anywhere from 2 to 20 eggs per needle . Whereas laying many eggs per needle may help diffuse host resin defenses and ensure that at least some eggs survive , increasing the number of eggs per needle also increases the risk of needle desiccation and egg death . Second, Neodiprion species vary in their tendency to cut a “preslit” on egg-bearing needles. A preslit is a non-egg-bearing cut located near the base of the pine needle, proximal to the eggs. This cut disrupts the resin canal, allowing for resin to drain from the needle . Like increasing egg number, the cutting of a preslit is associated with both increased oviposition success on resinous hosts  and increased desiccation risk . Finally, ovipositor morphology is highly variable across the genus and ovipositor characters are the primary traits used in species identification . Although the fitness consequences of different ovipositor shapes and sizes have not been studied, we hypothesize that different needle characteristics favor different ovipositor morphologies. For example, host plants with thin needles are likely to favor a reduction in ovipositor size.
To test the hypothesis that divergence in oviposition traits produces EPI between N. pinetum and N. lecontei, we evaluated a series of predictions. First, we predicted that N. pinetum and N. lecontei would have behavioral and morphological traits that are suited to the needle characteristics of their respective hosts. Specifically, we predicted that N. pinetum females would have oviposition traits that reduce the probability of needle desiccation (few eggs per needle, no pre-slits, and small ovipositors), while N. lecontei females would have oviposition traits that reduce the risk of drowning in host resin (many eggs per needle, pre-slits, and large ovipositors). Second, we predicted that hybrids and backcrosses would have reduced oviposition success (i.e., egg hatching) compared to each species, and that this reduction in success would be host dependent. Finally, if EPI is generated by oviposition traits, we predicted that oviposition success of backcrosses would be dependent on their oviposition traits. Together, our results provide compelling evidence that maladaptive combinations of oviposition traits contribute to extrinsic postzygotic isolation in Neodiprion lecontei and Neodiprion pinetum.
Insect collection and rearing
We collected N. pinetum and N. lecontei larvae throughout the eastern United States (Additional file 1: Table S1). We brought larvae back to the lab, transferred them to plastic boxes (32.4 cm × 17.8 cm × 15.2 cm) with mesh lids, and fed them pine foliage from their natal host species ad libitum. We collected cocoons as they were spun and stored them in individual gelatin capsules until adult emergence. We maintained all larvae and cocoons at 22 °C, 70% relative humidity, and an 18–6 h light-dark cycle [40, 46]. Upon emergence, live adults were stored at 4 °C until needed for crosses, morphological measurements, or behavioral assays. To propagate additional generations, we placed adult females and males into a mesh cage (35.6 cm × 35.6 cm × 61 cm) with seedlings of the pine species they were collected on. We allowed the adults to mate and the females to oviposit. After the eggs hatched, we reared larvae as described above.
Host needle width
N. pinetum uses Pinus strobus (white pine) exclusively, while N. lecontei has eight primary pine hosts (P. banksiana, P. resinosa, P. echinata, P. palustris, P. elliottii, P. rigida, P. taeda, and P. virginiana) [34, 47, 48]. To characterize the oviposition environment, we measured the widths of needles collected from ten trees from each of six Pinus species, including Pinus strobus and five of N. lecontei’s primary hosts (P. banksiana, P. resinosa, P. echinata, P. virginiana, P. rigida). Host collection locations are indicated in Additional file 1: Table S2. For each pine tree, we measured the width of ten needles using digital calipers (Mitutoyo CD-6”PMX), then averaged these values to produce an average needle width per tree. Although we did not measure needles from three of N. lecontei’s primary hosts (P. palustris, P. elliottii, P. taeda), we note that all three species have needle widths greater than 1 mm , which exceeds P. strobus needle widths (see below).
We used P. strobus and P. banksiana (jack pine) seedlings as hosts for oviposition in all experiments (purchased from Itasca Greenhouse in Cohasset, MN and North Central Reforestation, Inc. in Evansville, MN). We used P. banksiana because N. lecontei larvae from different populations tend to perform well on this pine species (personal observation; ), and because P. banksiana seedlings are readily available from nurseries. To assess how seedling needles (experimental hosts) compare to needles from mature hosts (typical hosts in nature), we measured needle widths for P. strobus and P. banksiana seedlings using the same approach described above. To analyze the differences in host needle width among mature pine species, and between mature hosts and seedlings, we performed two ANOVAs with Tukeys post hoc tests.
Although we didn’t measure the resin content for the pine hosts we used in our experiments, variation in resin canal number among pines is well characterized in the literature. Pines in the subgenus Strobus, to which P. strobus belongs, have 2–3 resin canals per needle [38, 51]. In comparison, all of the hosts N. lecontei uses are in the subgenus Pinus, which has 2–12 resin canals per needle [38, 51]. Thus, on average, we expect P. strobus needles to contain less resin than the Pinus species that are typically utilized by N. lecontei.
Our hypothesis that divergent selection has shaped oviposition traits in N. pinetum and N. lecontei predicted that N. pinetum (thin-needle specialist) would have a stronger preference for P. strobus and would lay fewer eggs per needle than N. lecontei (thick-needle specialist). We also predicted that, compared to N. lecontei, N. pinetum would cut fewer “preslits,” . We evaluated these predictions via a choice experiment. We first placed females in a clear 3.25-ounce deli cup with a single male until mating occurred. Neodiprion, like most hymenopterans, have arrhenotokous haplodiploidy, in which unfertilized eggs develop into haploid males [46, 52]. Although both mated and unmated females will oviposit readily (personal observation), we used mated females in our N. lecontei/N. pinetum oviposition assays as a means of propagating these lines. We then placed each mated female in a mesh cage (35.6 cm × 35.6 cm × 61 cm) with two P. banksiana seedlings and two P. strobus seedlings. We checked the cage daily until the female either oviposited or died. In nature, adult Neodiprion females have a very short life span (3–4 days) that is completely dedicated to reproduction [42, 47, 53]; likewise, in our choice cages, oviposition (or death) reliably occurred within 1–4 days (personal observation). For each female, we scored whether or not oviposition occurred. When oviposition occurred, we recorded host choice. Because N. pinetum and N. lecontei tend to cluster all of their eggs on a single branch terminus [48, 53], host choice is best described as a categorical trait with two possible outcomes: P. banksiana or P. strobus. We excluded three females that laid eggs on both hosts (representing 3.15% of the total sample).
To describe oviposition pattern, we counted the number of eggs, the number of egg bearing needles (EBN), and the number of EBN with preslits. We then used these data to calculate, for each female, the average number of eggs per EBN (number of eggs/number of EBN) and the proportion of EBN with preslits (number of EBN with preslits/total number of EBN). We then placed the egg-bearing seedling into an individual mesh sleeve cage (25.4 cm × 50.8 cm), and watered as needed until egg hatching occurred. To more fully describe oviposition pattern differences between the species, we also examined egg spacing, which may or may not influence hatching success on different hosts. To quantify egg spacing, we removed all EBN after hatching and imaged ten randomly selected EBN per seedling with a Canon EOS Rebel t3i camera equipped with an Achromat S 1.0X FWD 63 mm lens. Using these images, we measured the space between eggs in ImageJ  and, for each female, averaged egg spacing data across each needle. Sample sizes for each oviposition trait we scored are given in Additional file 1: Table S1.
To determine whether the two species differ in willingness to oviposit, we analyzed the proportion of N. pinetum or N. lecontei females that oviposited on any host with a generalized linear mixed effects model using a logit link factor, species as a fixed effect, and population (where each collecting location/host species combination was considered a separate population) as a random effect nested within species. To determine whether the two species differ in host preference, we used the same generalized linear mixed effects model to analyze the proportion of ovipositing females that chose P. strobus. To determine whether the two species differ in the average number of eggs per needle we used an ANOVA with species, natal host, and population as fixed factors. To determine whether the two species differ in average spacing between the eggs, we used an ANOVA. To determine whether the two species differ in the proportion of needles with preslits, we arcsine transformed the data and performed an ANOVA on the transformed data with species and population as fixed factors.
In addition to behavior, we examined ovipositor morphology, with the prediction that, compared to N. lecontei, N. pinetum would have smaller ovipositors. To characterize ovipositor morphology, we used five females from each of five populations from each species (N = 25 females per species; Additional file 1: Table S1). We used females preserved at -80 °C from either the parental phenotyping experiments or that had been frozen upon emergence. To control for body size, we measured the length of the forewing from the anterior junction of the forewing with the body to the tip of the forewing. We then removed the ovipositor, and mounted a single lancet (inner saw) using an 80:20 permount:toluene solution. We photographed each mounted lancet at 5× magnification using a Zeiss DiscoveryV8 stereomicroscope with an Axiocam 105 color camera and ZEN lite 2012 software (Carl Zeiss Microscopy, LLC Thornwood, NY). Using this software, we measured the length from the top of the second annulus to the top of the penultimate annulus, and measured width at the second annulus. We then performed morphometric analysis, which allows us to test for shape differences while controlling for size of the ovipositor. For this analysis, we placed nine landmarks and 21 sliding landmarks on each ovipositor (see “Results”). We then examined ovipositor shape using Geomorph . We applied a general procrustes alignment by minimizing binding energy. To determine whether the two species differed in ovipositor shape, we performed a procrustes ANOVA with forewing (as an allometric measurement), species, and population as fixed factors. To visualize the differences in ovipositor shape, we performed a principle components analysis in Geomorph. To determine whether the two species differed in ovipositor length or width, we used ANOVAs that included species, population nested within species, and forewing length as fixed factors. We completed all measurements and landmark placements in ImageJ version 1.49 V .
Cross oviposition behavior and success
We placed individual females of each cross-type in a choice cage as described above (see “Oviposition behavior”) and recorded whether or not oviposition occurred and, when it did occur, the preferred host. The backcross females were unmated, while all other females were mated. As we are specifically interested in reduced fitness due to oviposition traits, we used oviposition success as our measure of female performance. A female was considered to have “successful” oviposition if at least one of her eggs hatched and “unsuccessful” oviposition if no eggs hatched within 4 weeks. We chose 4 weeks as a cut-off because this is well beyond the typical egg development time for both species under our rearing conditions (generally <16 days, personal observation).
Finally, we attempted to recover every female as soon as possible after death or oviposition occurred. Recovered females were preserved in 100% EtOH and stored at -20 °C for future use. Sample sizes are located in Additional file 1: Table S3.
To determine whether the direction of F1 cross (i.e., F1LP vs. F1PL) differed in oviposition willingness, preference, or success, we used Z-tests. Because we did not observe any significant differences (see “Results”), we combined both cross-types into a single F1 category for the remaining analyses. To determine whether female cross-type (L, P, F1, BCL, BCP) differed in willingness to oviposit or in host preference, we used GLMs with a logit link factor and cross-type as a fixed effect, followed by post hoc Z-tests.
When there is postzygotic isolation, hybrids have reduced fitness compared to parental forms. To determine whether hybrids had reduced oviposition success compared to the parental species, we analyzed our hatch success data with a GLM using a logit link factor with cross-type as a fixed factor, followed by post hoc Z-tests. Additionally, if postzygotic isolation is “extrinsic” (due to the host plant), then oviposition success should be host-dependent. More specifically, each backcross type is expected to have the highest fitness (oviposition success) in the environment corresponding to the parent to which it is most similar genetically (i.e., there should be a cross-type-by-host interaction, Additional file 1: Figure S1) . To test these predictions, we used the same GLM model as for total postzygotic isolation, but added the interaction between cross-type and chosen host to the model. Because both BCL and BCP females would have been reared on whatever host their F1 mother chose, it is possible that rearing host may have influenced oviposition success of backcross females. To control for a possible rearing host effect on oviposition success of backcross females, we added rearing host as a fixed effect to a GLM model that included backcross type (BCL or BCP), oviposition host, and their interaction. We used Z-tests for all post hoc tests.
Impact of oviposition traits on BCL oviposition success
If oviposition traits are under selection and contribute to reduced hybrid fitness, the oviposition success of hybrid females should be dependent on these traits. To test these predictions, we focused on the BCL females because they were the only cross-type for which we had an appreciable sample size for all relevant traits (ovipositors, oviposition pattern, and hatching success). Additionally, because there was very little variation in hatch success on P. banksiana (see “Results”), we focused our analyses on P. strobus, with the prediction that the BCL females with pinetum-like traits (small ovipositor and few eggs per needle) would have the highest oviposition success on P. strobus. For these analyses, we scored oviposition success as a binary trait (hatch or no hatch) as described above. The results of the parental oviposition behavior assay indicated that N. pinetum females have a highly consistent and diagnostic oviposition pattern of three or fewer widely spaced eggs per needle (see below). To describe oviposition pattern, we therefore assigned each female to one of two categories: “pinetum,” if she laid three or fewer widely spaced eggs per needle and “non-pinetum,” if she laid more than three eggs per needle and/or eggs were spaced close together. To describe ovipositor morphology, we dissected and mounted female ovipositors as described above, with the addition of a rehydration step for EtOH-preserved females. The rehydration step consisted of six 10-min incubations of the female abdomen (at room temperature) in decreasing EtOH concentrations (100, 95, 80, 65, 50, and 25% EtOH), followed by overnight incubation in water.
To determine whether having a pinetum-like oviposition pattern increased the proportion of BCL females whose eggs hatched on P. strobus, we used a GLM with a logit link factor and oviposition pattern as a fixed factor. Next, to determine whether “successful” females had more pinetum-like ovipositors than “unsuccessful” females, we used Geomorph (procrustes ANOVA accounting for forewing length) to compare ovipositor shape between females with and without egg hatching on P. strobus . To determine whether ovipositor size affected the hatching rate on P. strobus, we performed separate GLMs with length and width data, both with a logit link factor and size as a fixed factor. To determine if oviposition pattern was correlated with ovipositor morphology we performed separate ANOVAs for ovipositor length and width, and a procrustes ANOVA for ovipositor shape.
Finally, we also used the BCL data to determine whether there was any relationship between host choice and oviposition traits, which may occur if females exhibit behavioral plasticity (e.g., alter oviposition behavior based on chosen host or alter host preference based on having a particular ovipositor morphology). To determine whether host choice correlates with oviposition pattern we used a GLM with a logit link factor and chosen host as a fixed factor. To determine whether host choice correlates with ovipositor morphology, we performed a procrustes ANOVA. To determine if ovipositor length and width correlated with oviposition host we performed two ANOVAs. All statistical analyses were performed using R version 3.2.3 .
Host needle width
Mature P. strobus had significantly thinner needles than all of the N. lecontei hosts (F5,54 = 72.42, P < 0.001, Additional file 1: Table S4, Fig. 1a). Likewise, P. strobus seedlings had significantly thinner needles than P. banksiana seedlings (P < 0.001, Additional file 1: Figure S2 and Table S5). However, because needles from P. banksiana seedlings were thinner than mature foliage (P < 0.001) and needles from P. strobus seedlings did not differ significantly from mature foliage (P = 0.15), the differences between our experimental hosts (F3,36 = 188.47, P < 0.001) are likely to be less extreme than differences typically experienced by ovipositing females in nature. In the discussion, we consider possible implications for the difference between seedling needles (experimental hosts) and mature needles (typical hosts).
Cross oviposition behavior and success
The direction of the F1 hybrid cross (i.e., F1LP vs. F1PL) had no effect on the female’s willingness to oviposit (Z = 1.20, P = 0.23), preference (Z = 1.20, P = 0.23), or oviposition success (Z = 1.01, P = 0.31). Given these findings, we combined F1 cross directions in subsequent analyses.
Post hoc tests (Z-tests) for interspecific crosses
Success on P. strobus
Success on P. banksiana
N. pinetum vs. Bcp
N. pinetum vs. F1 Hybrids
<1 × 10 −04
1.88 × 10 −03
N. pinetum vs. BCL
<1 × 10 −04
3 × 10 −04
4.66 × 10 −03
N. pinetum vs. N, lecontei
<1 × 10 −04
BCP vs. F1 Hybrid
1.58 × 10 −03
BCP vs. BCL
<1 × 10 −04
<1 × 10 −04
BCP vs. N. lecontei
<1 × 10 −04
F1 Hybrid vs BCL
<1 × 10 −04
<1 × 10 −04
F1 Hybrid vs. N. lecontei
3 × 10 −04
4.8 × 10 −04
BCL vs. N. lecontei
3 × 10 −04
6.34 × 10 −03
The cross-types also differed in their oviposition success (χ2 4 = 13.03, P = 0.011, Fig. 6c), and the F1 females had significantly lower hatching success than any of the other cross-types (Table 1). When oviposition host and an interaction between host and cross-type were added, cross-type remained significant (χ2 4 = 14.92, P = 0.0049, Fig. 6d). Additionally, there was a significant effect of host on oviposition success (χ2 1 = 44.43, P <0.001): across all cross-types, females that chose P. strobus had lower hatching success than females that chose P. banksiana (Fig. 6d). Also, although none of the N. lecontei females involved in the cross chose P. strobus, four of the N. lecontei females from our multi-population preference experiment did chose P. strobus (Fig. 3b). Notably, all four of these females experienced complete hatching failure (Additional file 1: Figure S3). Although both cross-type and oviposition host significantly impacted hatching success, the interaction between them was not significant (χ2 3 = 2.37, P = 0.50). We also found that rearing host did not affect the oviposition success of backcross females (χ2 1 = 0.22, P = 0.64).
Impact of oviposition traits on BCL oviposition success
BCL females that had a pinetum-like oviposition pattern were significantly more likely to have eggs that hatched on P. strobus than if they deviated from this pattern (χ2 1 = 3.85, P =0.0498, Fig. 7a). Also, BCL females that successfully oviposited on P. strobus had significantly shorter ovipositors than unsuccessful females (χ2 1 = 9.50, P = 0.0021, Fig. 7b). In contrast, successful and unsuccessful females did not differ in ovipositor width (χ2 1 = 0.019, P = 0.89) or ovipositor shape (F1, 17 = 1.16, P = 0.24).
Empirical data from diverse taxa indicate that ecological speciation is common in nature [7–10], and that changes in host use frequently initiate ecological speciation in plant-feeding insects [27, 29]. However, the contributions of specific divergent traits to EPI are unknown in most systems. In this study, we evaluated evidence of oviposition traits generating extrinsic postzygotic isolation between a pair of Neodiprion sawfly species that specialize on different pines. We found compelling evidence of EPI stemming from maladaptive combinations of oviposition traits. Here, we discuss the limitations, as well as broader implications of our work for ecological specialization and speciation in plant-feeding insects and future research directions in this promising empirical system.
Although all sawflies in the genus Neodiprion feed on host plants in the family Pinaceae (mostly in the genus Pinus), different sawfly species tend to specialize on different pine hosts [34, 45]. Previous analyses at both the inter- and intraspecific levels indicate that changes in host use are associated with population differentiation and speciation in this genus [25, 35]. In this study, we investigated a potential causal relationship between adaptation to different hosts and reproductive isolation. In particular, we hypothesized that maladaptive combinations of divergent oviposition traits give rise to extrinsic postzygotic isolation between Neodiprion species. Consistent with this hypothesis, we found that sister species N. pinetum (a thin-needled specialist) and N. lecontei (occurs on thicker-needled hosts) differed in multiple behavioral and morphological traits related to oviposition (Figs. 3, 4 and 5). In terms of behavior, N. pinetum females preferred P. strobus (white pine), laid a small number of widely spaced eggs on each needle, and never cut resin-draining preslits. In contrast, N. lecontei females generally avoided P. strobus, laid many closely spaced eggs per needle, and almost always cut preslits. In terms of morphology, N. pinetum females had smaller, straighter ovipositors than N. lecontei females. Together, N. pinetum traits likely enable females to insert eggs into P. strobus without damaging the thin needles to the point that they dry out and the eggs die, while N. lecontei traits should better equip females to circumvent host defenses and prevent eggs from being overwhelmed by resin.
Although N. lecontei and N. pinetum appear to be specialized to oviposit on different hosts, they do hybridize in nature (personal observation, ), indicating that premating barriers are incomplete. Nevertheless, the strong genetic, behavioral, and morphological differentiation between these two sympatric species (; Figs. 3, 4 and 5) suggests that there are postzygotic barriers to gene exchange. Consistent with this prediction, we found that F1 females had reduced oviposition success relative to the two parental species (Fig. 6c). For these females, there were two potential sources of oviposition failure: botched oviposition (which would be host-dependent and therefore extrinsic in nature) and egg inviability (which could stem from intrinsic genetic incompatibilities or from extrinsic egg-host interactions). Our observation that hybrid females had reduced oviposition success only when they chose P. strobus suggests that postzygotic isolation between N. lecontei and N. pinetum is largely attributable to extrinsic, rather than intrinsic, factors. By contrast, oviposition success of hybrid females on the more “benign” P. banksiana seedlings was indistinguishable from oviposition success of pure N. lecontei and N. pinetum females. Although this finding is consistent with EPI, an alternative explanation for these results is that intrinsic genetic incompatibilities between the species are more pronounced in the P. strobus environment [7, 13, 18]. One way to control for intrinsic genetic incompatibilities is to compare the fitness of both backcross types in both parental environments . Using this method, we found that BCP females had high oviposition success on both hosts, while BCL females had high oviposition success on P. banksiana only (Fig. 6d). While patterns observed in BCL females are consistent with EPI, patterns observed in BCp are not.
While seemingly at odds with predictions under EPI, our observation that BCP females had high oviposition success on both hosts could be attributable to our experimental design. There are two main sources of experimental error that could have precluded us from detecting reduced hatch success on P. banksiana. First, by scoring oviposition success as a binary trait (hatch or no hatch), we lumped together females with a wide range of hatching success (from <10 to 100%). Failure to account for variation in non-zero hatch success would have reduced our power to detect all but the most extreme differences in oviposition success. In other words, while lecontei-like oviposition traits led to a complete failure on P. strobus consistently enough that we could detect it with our crude measure of success, we had little power to detect subtler reductions in hatching success.
The second potential source of experimental error in our assessment of EPI is that we used pine seedlings in lieu of larger trees, which we could not accommodate in our growth rooms. However, as our needle width data indicate (Additional file 1: Figure S2), the pine seedlings we used did not fully recapitulate differences in the host age classes that are typically selected by ovipositing N. lecontei and N. pinetum females in the wild [53, 57]. In particular, the needles of our P. banksiana seedlings were considerably thinner than needles from older trees. Moreover, resin content tends to increase as pine trees age . Thus, while the P. strobus seedlings we used replicated the challenge of laying eggs on a thin needled-host, the P. banksiana seedlings did not replicate the challenge of laying on a thick, resinous needle.
Despite these possible experimental artifacts, we do have an additional line of direct evidence supporting the existence of EPI due to oviposition traits: on P. strobus, BCL females with lecontei-like oviposition traits (ovipositor morphology and egg-laying behavior) had reduced oviposition success compared to BCL females with pinetum-like oviposition traits (Fig. 7). Because all BCL females share the same genetic makeup (i.e., same proportion of N. lecontei and N. pinetum alleles), these differences cannot be explained by intrinsic genetic incompatibilities. Taken together, our cross data indicate that maladaptive combinations of oviposition preference and oviposition traits in hybrids generate EPI between N. lecontei and N. pinetum. Intriguingly, maladaptive combinations of preference and performance traits have been reported in several other insect taxa [20, 29, 59, 60], suggesting that this might be a widespread cause of reduced hybrid fitness.
Our analysis of traits in BC females also demonstrates how examination of specific traits in hybrid individuals can be used as an alternative to the “modify-parental-phenotype” test of EPI that has been proposed, but never utilized . In our case, modifying parental phenotypes was not an option because our focal phenotypes were either behavioral (host preference, oviposition pattern) or involved a delicate morphological structure (ovipositor) that we could not alter readily—we suspect that the same is probably true of many organisms in which one might want to investigate EPI. However, as we have shown here, genetic crosses can serve a similar function as parental modification. In particular, by generating recombination among loci underlying ecologically relevant traits and assessing fitness in recombinant individuals, we could begin to tease apart how individual traits and interactions between them contribute to reduced fitness of hybrids in parental environments. To date, we know of only one other study that has taken advantage of trait variation in hybrids to make inferences regarding EPI in plant-feeding insects: McBride and Singer’s  study of EPI in Euphydryas butterflies (see also  for an example in Caribbean pupfishes). In their study, McBride and Singer reared F1 hybrids between allopatric, host-specialized populations on both parental hosts and, for four behavioral traits, found that trait intermediacy in the hybrids reduced their fitness on both hosts.
To date, numerous studies—many of which focused on plant-feeding insects—have reported evidence of EPI (see [8, 28, 29]). While only a handful of these have employed a more rigorous approach (e.g., reciprocal backcross or trait-focused studies) that controls for genetic incompatibilities [20, 30–33], the emerging picture from this body of work is that EPI frequently accompanies ecological speciation. However, in only a handful of cases have the traits underling EPI been identified [18, 20]. Importantly, although EPI is a direct consequence of adaptive divergence, adaptive divergence does not always produce EPI. For example, if intermediate trait values do not impact fitness in parental environments or if an intermediate environment is available in nature, hybrids will not experience ecologically based reductions in fitness . As more traits are explicitly tested for their role in EPI, we can begin to ask more specific questions about its mechanistic basis, such as: which traits (behavior, physiology, morphology) and which aspects of ecology (reproduction, food acquisition and processing, parasitism) are most likely to produce EPI?
Based on their findings in Euphydryas butterflies, McBride and Singer  proposed that behavioral traits—especially niche preferences—might be especially important drivers of EPI. In support of this argument they provided two additional examples. First, two European blackcap populations that migrate in opposite directions to their wintering grounds produce hybrids with a tendency to migrate in an intermediate and maladaptive direction . Second, hybrids between apple- and hawthorn host races of Rhagoletis pomonella have a tendency to avoid both parental hosts, making it difficult for them to locate suitable oviposition sites [62, 63]. By contrast, our hybrids did not exhibit a reduction in willingness to oviposit (in fact, for reasons that are currently unclear to us, BCL seemed more willing to oviposit than other cross types; Fig. 6), indicating that “host confusion” is not contributing to EPI in this system. Nevertheless, our data are consistent with the overall importance of behavioral traits (in our case, host preference and oviposition pattern) in driving EPI.
Additionally, similar to our finding that oviposition traits contribute to EPI in Neodiprion, three of the four traits implicated in reduced hybrid fitness in Euphydryas butterflies were related to oviposition. Experimental and natural history work in additional Neodiprion species suggest that this phenomenon might be widespread in pine sawflies as well. Specifically, for several different Neodiprion species, there are published observations of egg mortality caused by either needle desiccation or drowning in resin, and these outcomes are often associated with needle thickness, needle resin content, and female oviposition pattern (e.g., number of eggs per needle, presence of preslit, depth of egg slits and preslits) [41–44, 64–67]. Together, these observations suggest that oviposition traits are under strong selection both within and between Neodiprion species. Intriguingly, host preference, ovipositor morphology, and oviposition pattern are also among the most variable traits in the genus and are often useful in species identification [45, 68, 69]. If host-related selection has shaped this variation, inter- and intraspecific variation in host preference should correlate with variation in other oviposition traits; this prediction could be tested using a comparative approach.
Beyond Neodiprion, oviposition-related traits—which include traits related to finding and choosing a host, selecting a site within the host for egg deposition, depositing eggs in specific patterns on or within the host tissue, defusing host defenses, ovipositor morphology, and egg morphology—could profoundly impact the fitness of any egg-laying phytophagous insect female and are therefore likely to be frequent targets of natural selection . In support of this argument, numerous studies have reported host-associated differentiation in oviposition traits, including: clutch size in seed beetles , ovipositor morphology in yucca moths , ovipositor length in gall-inducing Asphodylia flies , ovipositor length in fig wasps , ovipositor size in Plateumaris leaf beetles , clutch size and oviposition site in butterflies , and multiple morphological and behavioral traits in pine sawflies (this study). However, in the context of traits driving ecological specialization and speciation in plant-feeding insects, research has focused almost exclusively on female host preference and larval performance (i.e., growth and survival rates when feeding on a particular host plant). To understand the role of host specialization in phytophagous insect speciation, it is critical that we examine additional host-related traits.
In this study, we have demonstrated that oviposition traits contribute to EPI between N. lecontei and N. pinetum. While these observations are consistent with ecological speciation, the evidence is not yet iron-clad and many important questions remain. First, while we focused here on oviposition traits, other traits—such as larval performance—could also contribute to EPI. In future work, we hope to quantify the impact of individual traits—and the interaction between them—on host-dependent reductions in hybrid fitness. Second, while we have focused here on EPI, there are other sources of reproductive isolation between these species (personal observation). To evaluate the contribution of EPI to total isolation, we must quantify the strength of EPI relative to other reproductive barriers [77, 78]. Finally, understanding how divergent traits give rise to RI requires that we identify the genetic mechanisms (i.e., linkage or pleiotropy) linking them [7, 8]. As we have demonstrated here, these species are interfertile in the lab; thus, a QTL mapping approach is feasible in this system. Additionally, identification of causal loci—which is required if we are to distinguish between pleiotropy and linkage—will be facilitated by the availability of annotated genome assemblies for N. lecontei  and N. pinetum (in progress).
While a long-term goal is to identify all host-related traits under selection, all reproductive barriers, and their underlying genes in N. lecontei and N. pinetum, these efforts will provide a single snapshot at one time point in speciation. Because these species have been diverging for up to several million years [25, 80], they have had time to accumulate many differences and barriers to reproduction, which will make it difficult to determine which reproductive barriers arose first. To get at this question, we can examine other Neodiprion species and populations at different stages along the “speciation continuum” . For example, there is evidence of host-associated differentiation in at least two Neodiprion species (Neodiprion lecontei; ; Neodiprion abietis, [42, 82]), and possibly other Neodiprion species as well. Although much work remains, extensive natural history data, experimental tractability, and growing genomic resources make Neodiprion an exceptionally rich system for addressing many long-standing questions regarding the evolution of host specialization and its role in generating the staggering diversity of phytophagous insects.
Egg bearing needle
Extrinsic postzygotic isolation
We would like to thank Ethan Adams, Julianne Horn, Melanie Hurst, and Anna Sosso for their assistance with data collection, and Robin K. Bagley for specimen collection. We thank members of the Linnen lab for insect rearing assistance. We also thank members of the Linnen lab, two anonymous reviewers, and Jeff Feder for helpful comments and discussion on previous versions of this manuscript.
This work was supported by the University of Kentucky and the National Science Foundation (DEB-1257739).
Availability of data and materials
All data (needle widths; oviposition behavior, ovipositor morphology, and hatch success for parents, F1s, and BCs) are archived at Dryad doi:10.5061/dryad.4g26c.
EEB, KLV, and CRL participated in the study design. EEB and KLV participated in the data collection. EEB and CRL analyzed the data. EEB and CRL drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Darwin C. On the origins of species by means of natural selection. London: John Murray; 1859.Google Scholar
- Mayr E. Ecological factors in speciation. Evolution. 1947;1:263–88.View ArticleGoogle Scholar
- Mayr E. Systematics and the origin of species, from the viewpoint of a zoologist. New York: Columbia University Press; 1942.Google Scholar
- Dobzhansky T. Genetics and the origin of species. 3rd ed. New York: Columbia University Press; 1951.Google Scholar
- Schluter D. Ecology and the origin of species. Trends Ecol Evol. 2001;16:372–80.View ArticlePubMedGoogle Scholar
- Schluter D. The ecology of adaptive radiation. Oxford: Oxford University Press; 2000.Google Scholar
- Rundle HD, Nosil P. Ecological speciation. Ecol Lett. 2005;8:336–52.View ArticleGoogle Scholar
- Nosil P. Ecological Speciation. Oxford: Oxford University Press; 2012.
- Schluter D. Evidence for Ecological Speciation and Its Alternative. Science. 2009;323:737–41.
- Van der Niet T, Peakall R, Johnson SD. Pollinator-driven ecological speciation in plants: new evidence and future perspectives. Ann Bot. 2014;113:199–211.View ArticlePubMedPubMed CentralGoogle Scholar
- Funk DJ, Nosil P, Etges WJ. Ecological divergence exhibits consistently positive associations with reproductive isolation across disparate taxa. Proc Natl Acad Sci U S A. 2006;103:3209–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Rice W, Hostert E. Laboratory experiments on speciation - what have We learned in 40 years. Evolution (N Y). 1993;47:1637–53.Google Scholar
- Rundle HD, Whitlock MC. A genetic interpretation of ecologically dependent isolation. Evolution. 2001;55:198–201.View ArticlePubMedGoogle Scholar
- Coyne J, Orr H. Speciation. Sunderland: Sinauer Associates, Inc; 2004.Google Scholar
- Seehausen O, Butlin RK, Keller I, Wagner CE, Boughman JW, Hohenlohe PA, et al. Genomics and the origin of species. Nat Rev Genet. 2014;15:176–92.View ArticlePubMedGoogle Scholar
- Dobzhansky T. Speciation as a stage in evolutionary divergence. Am Nat Science Press. 1940;74:312–21.Google Scholar
- Nosil P, Vines TH, Funk DJ. Perspective reproductive isolation caused by natural selection against immigrants from divergent habitats. Evolution. 2005;59:705–19.PubMedGoogle Scholar
- Hatfield T, Schluter D. Ecological speciation in sticklebacks : environment-dependent hybrid fitness. Evolution. 1999;54:866–73.View ArticleGoogle Scholar
- Martin CH, Wainwright PC. Multiple fitness peaks on the adaptive landscape drive adaptive radiation in the wild. Science. 2013;339:208–11.View ArticlePubMedGoogle Scholar
- Mcbride CS, Singer MC. Field studies reveal strong postmating isolation between ecologically divergent butterfly populations. PLoS Biol. 2010;8:e1000529.View ArticlePubMedPubMed CentralGoogle Scholar
- Mitter C, Farrell B, Wiegmann B. The phylogenetic study of adaptive zones: has phytophagy promoted insect diversification? Am Nat. 1988;132:107–28.View ArticleGoogle Scholar
- Farrell BD. “Inordinate fondness” explained: Why are there so many beetles? Sci. 1998;281:555–9.View ArticleGoogle Scholar
- Wiens J, Lapoint R, Whiteman N. Herbivory increases diversification across insect clades. Nat Commun. 2015;6:8370.
- Winkler IS, Mitter C. The Phylogenetic Dimension of Insect-Plant Interactions: A Review of Recent Evidence. In: Tilmon KJ. Specialization, Speciation, Radiation: The Evolutionary Biology of Herbivorous Insects. Berkeley: University of California Press; 2008. p. 240–63.
- Linnen CR, Farrell BD. A test of the sympatric host race formation hypothesis in Neodiprion (Hymenoptera: Diprionidae). Proc Biol Sci. 2010;277:3131–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Nyman T, Vikberg V, Smith D. How common is ecological speciation in plant-feeding insects? A‘Higher’Nematinae perspective. BMC. 2010;10:266.Google Scholar
- Berlocher SH, Feder JL. Sympatric speciation in phytophagous insects: moving beyond controversy ? Annu Rev Entomol. 2002;47:773–815.View ArticlePubMedGoogle Scholar
- Funk DJ, Filchak KE, Feder JL. Herbivorous insects: model systems for the comparative study of speciation ecology. Genetica. 2002;116:251–67.View ArticlePubMedGoogle Scholar
- Matsubayashi KW, Ohshima I, Nosil P. Ecological speciation in phytophagous insects. Entomol Exp Appl. 2010;134:1–27.View ArticleGoogle Scholar
- Rundle HD. A test of ecologically dependent postmating isolation between sympatric sticklebacks. Evolution. 2002;56:322–9.View ArticlePubMedGoogle Scholar
- Egan SP, Funk DJ. Ecologically dependent postmating isolation between sympatric host forms of Neochlamisus bebbianae leaf beetles. Proc Natl Acad Sci U S A. 2009;106:19426–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Kuwajima M, Kobayashi N. Detection of ecological hybrid inviability in a pair of sympatric phytophagous ladybird beetles (Henosepilachna spp.). Entomol Experimentalis et applicata. 2010;134:280–6.View ArticleGoogle Scholar
- Soudi S, Reinhold K, Engqvist L. Ecologically dependent and intrinsic genetic signatures of postzygotic isolation between sympatric host races of the leaf beetle Lochmaea capreae. Evolution. 2016;70:471–9.View ArticlePubMedGoogle Scholar
- Coppel H, Benjamin D. Bionomics of the- nearctic pine-feeding diprionids. Annu Rev Entomol. 1965;10:69–96.View ArticleGoogle Scholar
- Bagley RK, Sousa VC, Niemiller ML, Linnen CR. History, geography, and host use shape genome-wide patterns of genetic differentiation in the redheaded pine sawfly (Neodiprion lecontei). Mol Ecol. In Press.
- Linnen CR, Farrell BD. Mitonuclear discordance is caused by rampant mitochondrial introgression in Neodiprion (Hymenoptera: Diprionidae) sawflies. Evolution. 2007;61:1417–38.View ArticlePubMedGoogle Scholar
- Linnen CR, Farrell BD. Comparison of methods for species-tree inference in the sawfly genus Neodiprion (Hymenoptera: Diprionidae). Syst Biol. 2008;57:876–90.View ArticlePubMedGoogle Scholar
- Wu H, Hu Z. Comparative anatomy of resin ducts of the Pinaceae. Trees. 1997;11:135–43.View ArticleGoogle Scholar
- Wilkinson RC, Becker GC, Benjamin DM. The biology of Neodiprion rugifrons (Hymenoptera: Diprionidae), a sawfly infesting jack pine in Wisconsin. Ann Entomol Soc. 1966;59:786–92.View ArticleGoogle Scholar
- Knerer G. Morphological and physiological clines in Neodiprion pratti (Dyar) (Symphyta, Diprionidae) in eastern North America. J Appl Entomol. 1984;97:9–21.Google Scholar
- Wilkinson R. Slash-pine sawfly, Neodiprion merkeli. 1. oviposition pattern and descriptions of egg, female larva, pupa, and cocoon. Ann Entomol Soc. 1971;64:241–7.View ArticleGoogle Scholar
- Knerer G, Atwood CE. Diprionid sawflies: polymorphism and speciation. Science. 1973;179:1090–9.View ArticlePubMedGoogle Scholar
- McCullough D, Wagner MR. Defusing host defenses: ovipositional adaptations of sawflies to plant resin. Sawfly life Hist. Adapt. to woody plants. 1993. p. 151–71.Google Scholar
- Codella SG, Raffa KF. Desiccation of Pinus foliage induced by conifer sawfly oviposition: effect on egg viability. Ecol Entomol. 2002;27:618–21.View ArticleGoogle Scholar
- Ross H. The Taxonomy and Evolution of the Saw Genus Neodiprion. For Sci. 1955;1:196–209.
- Harper K, Bagley RK, Thompson K, Linnen C. Complementary sex determination, inbreeding depression and inbreeding avoidance in a gregarious sawfly. 2016.Google Scholar
- Benjamin DM. The Biology and Ecology of the Red-headed Pine Sawfly. Washington, D.C.: U.S. Department of Agriculture; 1955.
- Wilson LF, Wilkinson RC, Averill RC. Redheaded pine sawfly- its ecology and management. 1992.Google Scholar
- Stults DZ, Axsmith BJ, Liu Y-S. Evidence of white pine (Pinus subgenus Strobus) dominance from the Pliocene Northeastern Gulf of Mexico Coastal Plain. Palaeogeogr Palaeoclimatol Palaeoecol. 2010;287:95–100.View ArticleGoogle Scholar
- Knerer G. Diprionid sawflies: biological topics and rearing techniques (hymenoptera: symphyta). Bull Entomol Soc Am. 1984;30:53–7.Google Scholar
- Gernandt DS, Lopez GG, Garcia SO, Liston A. Phylogeny and classification of Pinus. Taxon. 2005;54:29–42.View ArticleGoogle Scholar
- Heimpel G, Boer J de. Sex determination in the Hymenoptera. Annu Rev Entomol. 2008;53:209–30.
- Rauf A, Benjamin DM. The biology of the white pine sawfly, Neodiprion pinetum (Hymenoptera: Diprionidae) in Wisconsin. Gt Lakes Entomol. 1980;13:219–24.Google Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved.View ArticlePubMedGoogle Scholar
- Adams DC, Otárola-Castillo E. Geomorph: an r package for the collection and analysis of geometric morphometric shape data. Methods Ecol Evol. 2013;4:393–9.View ArticleGoogle Scholar
- R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R foundation for statistical computing; 2015.
- Averill RD, Wilson LF, Fowler RF. Impact of the redheaded pine sawfly (hymenoptera: diprionidae) on young red pine plantations. Gt Lakes Entomol. 1982;15:65–96.Google Scholar
- Lin J, Sampson DA, Ceulemans R. The effect of crown position and tree age on resin-canal density in Scots pine ( Pinus sylvestris L.) needles. Can J Bot. 2001;79:1257–61. NRC Research Press Ottawa, Canada.Google Scholar
- Forister ML. Independent inheritance of preference and performance in hybrids between host races of Mitoura butterflies (Lepidoptera: Lycaenidae). Evolution. 2005;59:1149–55.View ArticlePubMedGoogle Scholar
- Ohshima I. Host race formation in the leaf-mining moth Acrocercops transecta (Lepidoptera: Gracillariidae). Biol J Linn Soc. 2008;93:135–45.
- Helbig A. Inheritance of migratory direction in a bird species: a cross-breeding experiment with SE- and SW-migrating blackcaps (Sylvia atricapilla). Behav Ecol Sociobiol. 1991;28:9–12. Springer-Verlag.View ArticleGoogle Scholar
- Linn C, Dambroski H, Feder J. Postzygotic isolating factor in sympatric speciation in Rhagoletis flies: reduced response of hybrids to parental host-fruit odors. Proc Natl Acad U S A. 2004;101:17753–8.View ArticleGoogle Scholar
- Forbes AA, Fisher J, Feder JL. Habitat avoidance: overlooking an important aspect of host-specific mating and sympatric speciation? Evolution. 2005;59:1552–9.View ArticlePubMedGoogle Scholar
- Warren LO, Coyne JF. The pine sawfly, Neodiprion taedae linearis Ross, in Arkansas. Univ Arkansas Agric Exp Stn Bull. 1958;602:1–23.
- Martineau R. On an infestation of the red-headed jack pine sawfly, Neodiprion virginianus complex in Quebec. Can Dep Agric For Biol Div Bimon Prog Rep. 1959;15:1.
- Kapler JE, Benjamin DM. The biology and ecology of the red-pine sawfly in Wisconsin. For Sci. 1960;6:253–68.Google Scholar
- Wilkinson RC. The biology and ecology of the Neodiprion virginianus complex in Wisconsin. Madison: University of Wisconsin-Madison; 1961.
- Ghent AW. Row type oviposition in Neodiprion sawflies as exemplified by the european pine sawfly, N. sertifer (Geoff). Can J Zool. 1959;37:267–81.View ArticleGoogle Scholar
- Linnen CR, Smith D. Recognition of Two Additional Pine-Feeding Neodiprion Species (Hymenoptera: Diprionidae) in the Eastern United States. Proc Entomol Soc Washingt. 2012;114:492–500.Google Scholar
- Ni J. Evolutionary ecology of oviposition strategies. In: Hilker M, Meiners T, editors. Chemoecology insect eggs egg depos. Berlin: Blackwell Publishing Ltd; 2003. p. 349–76.Google Scholar
- Messina FJ, Karren ME. Adaptation to a novel host modifies host discrimination by the seed beetle Callosobruchus maculatus. Anim Behav. 2003;65:501–7.View ArticleGoogle Scholar
- Groman JD, Pellmyr O. Rapid evolution and specialization following host colonization in a yucca moth. J Evol Biol. 2000;13:223–36.View ArticleGoogle Scholar
- Joy JB, Crespi BJ. Adaptive radiation of gall-inducing insects within a single host-plant species. Evolution. 2007;61:784–95.View ArticlePubMedGoogle Scholar
- Weiblen GD, Bush GL. Speciation in fig pollinators and parasites. Mol Ecol. 2002;11:1573–8.View ArticlePubMedGoogle Scholar
- Sota T, Hayashi M, Yagi T. Geographic variation in body and ovipositor sizes in the leaf beetle Plateumaris constricticollis (Coleoptera: Chrysomelidae) and its association with climatic conditions and host plants. Eur J Entomol. 2007;104:165–72.View ArticleGoogle Scholar
- Singer MC, McBride CS. Multitrait, host-associated divergence among sets of butterfly populations: implications for reproductive isolation and ecological speciation. Evolution. 2010;64:921–33.View ArticlePubMedGoogle Scholar
- Sobel JM, Chen GF. Unification of methods for estimating the strength of reproductive isolation. Evolution. 2014;68:1511–22.View ArticlePubMedGoogle Scholar
- Nosil P. Divergent host plant adaptation and reproductive isolation between ecotypes of Timema cristinae walking sticks. Am Nat. 2007;169:151–62.View ArticlePubMedGoogle Scholar
- Vertacnik KL, Linnen CR. Neodiprion lecontei genome assembly Nlec1.0. NCBI/GenBank.
- Linnen CR, Farrell BD. Phylogenetic analysis of nuclear and mitochondrial genes reveals evolutionary relationships and mitochondrial introgression in the sertifer species group of the genus Neodiprion (Hymenoptera: Diprionidae). Mol Phylogenet Evol. 2008;48:240–57.View ArticlePubMedGoogle Scholar
- Nosil P, Harmon LJ, Seehausen O. Ecological explanations for (incomplete) speciation. Trends Ecol Evol. 2009;24:145–56.View ArticlePubMedGoogle Scholar
- Knerer G, Atwood CE. Evolutionary trends in the subsocial sawflies belonging to the Neodiprion abietis complex (Hymeoptera: Tentherdinoidea). Amer Zool. 1972;12:407–18.View ArticleGoogle Scholar