Spatial turnover in host-plant availability drives host-associated divergence in a South African leafhopper (Cephalelus uncinatus)
© The Author(s). 2017
Received: 5 October 2016
Accepted: 17 February 2017
Published: 9 March 2017
The evolution of reproductive isolation between herbivorous insect populations is often initiated by shifts to novel host-plants, a process that underlies some of the best examples of ecological speciation. However, it is not well understood why host-shifts occur. Arguably the most common hypothesis is that host-shifts occur in response to competition, while a less frequently invoked hypothesis is that herbivores adapt locally to geographic differences in potential host-plant communities. Here we investigate whether geographic variation in host-plant availability is likely to have driven host-shifts in restio leafhoppers. We studied local adaptation of a camouflaged restio leafhopper species, Cephalelus uncinatus, to host-plants in the Restionaceae (restios); a family of plants with exceptional diversity in the anomalously species-rich Cape Floristic Region (CFR). To determine whether C. uncinatus experiences heterogeneous host communities across its range, we first quantified the degree of geographic overlap between C. uncinatus and each of its associated host-plant species. Then we quantified trait divergence (host preference, body size and colour) for three pairs of C. uncinatus populations found on different host-plant species differing in their degree of spatial overlap. Spectral reflectance was modelled in bird visual space to investigate whether body colour divergence in C. uncinatus corresponds to leaf sheath colour differences between restio species as perceived by potential predators.
We demonstrate that C. uncinatus is forced to use different restio species in different regions because of turnover in available host species across its range. Comparisons between geographically separated populations were consistent with local adaptation: restio leafhoppers had preferences for local host-plants over alternative host-plants and matched local plants better in terms of size and colour.
Spatial turnover in host-plant availability has likely facilitated host-shifts in C. uncinatus. Spatial turnover in host-plant availability may be an important driver of insect diversification in the CFR and globally.
KeywordsBeta diversity Ecological divergence Geographic mosaic Host-shift Interaction turnover Local adaptation Magic traits Phytophagous insects Specialization Speciation
A simpler, but surprisingly, less commonly invoked hypothesis for host shifts is that spatially separated populations of insects experience different host-plants. Local adaptation occurs in allopatry because the different hosts have different chemical or morphological characteristics. This may occur if herbivores expand their ranges into areas with novel host assemblages or if host community composition changes in parts of the existing range of herbivorous insects. Once a population establishes on the novel host-plant, local adaptation can lead to a shift in host preference  and lower survival on the ancestral host relative to the novel host . Consequently, specialization on different host-plant species can lead to phenotypically divergent host-ecotypes which become reproductively isolated through assortative host-use and immigrant inviability [21, 22]. We refer to this process as the allopatric host-shift model (Fig. 1b). The pattern produced by the sympatric host-shift model may not be easily distinguished from the allopatric host-shift model if range expansion occurs after a sympatric host-shift (compare Fig. 1a, b). However, the allopatric host-shift model is dependent on high spatial heterogeneity in host-plant availability, whereas the sympatric model is not.
Spatial environmental heterogeneity (for example, in soils, topography, pollinators) is thought to underlie plant diversification in the Cape Floristic Region (CFR) , a biodiversity hotspot for plants . In addition, plant diversity has largely accrued through spectacular radiations of a limited number of plant clades . This mode of allopatric plant speciation has resulted in exceptional species turnover between communities, evidenced by the unusually high beta diversity that characterises the CFR . However, these spatially separated plant communities, differing strikingly in species composition, may often comprise phylogenetically closely related species from radiating clades [27–29]. Climatic stability in the region is also thought to have facilitated plant diversification by limiting population extinction and providing ample time for divergence to take place  across stable environmental gradients. Together these characteristics of the CFR flora suggest that an allopatric host-shift model of host-ecotype formation is a likely scenario for diversification of the associated insect herbivore lineages. Strong and temporally stable geographic mosaics of host-plant availability are likely to result in divergent selection across spatially separated populations of insect species. Furthermore, the presence of related plant species across spatially separated communities may facilitate the survival of insects dispersing outside the range of their ancestral host-plants. Therefore, spatial heterogeneity is likely to be important in causing allopatric host-shifts in herbivorous insects in the CFR, but it remains to be investigated.
In contrast, fire, a key driver of CFR ecology, may reduce the likelihood of host-shifts in sympatry. Fire disturbance occurs frequently in this region, and the main vegetation type, Fynbos, consists of plants showing a variety of adaptations to survive fires (reviewed in ). While some insect species can survive Fynbos fires (reviewed in ), many species must recolonise burned vegetation (probably through dispersal from neighbouring unburned vegetation patches), which can take 3 years . Additionally, local plant communities differ stochastically across fire cycles . Divergence between different host-associated insect populations in sympatry may be impeded by these consequences of fire, i.e. regular resetting of population size and insect community composition reducing the likelihood of strong competition, frequent extinction of locally adapting insect populations, and temporally variable selection resulting from both the fire-driven local dispersal/re-colonisation cycle and fire-induced changes in local plant composition .
One of the diagnostic plant families within the fynbos is the African Restionaceae (restios hereafter). Restios form a highly diverse, monophyletic plant clade, which comprises about 350 species. They are thought to have originated about 65 million years ago , making them one of the oldest clades in the CFR with only 10 species occurring outside the CFR (reviewed in ). Restios are generally reed-like in appearance. Their photosynthetic stems (culms) have regular nodes with dried-out leaf sheaths that persist in most species, but drop off in others. The leaf sheaths of restios appear to be mimicked by the morphology of locally host-specific herbivores called restio leafhoppers (Cicadellidae: Cephalelini) , the dominant insects on restios . Restio leafhoppers are characterised by small, slender bodies and elongated crowns (with the exception of Duospina capensis) resembling the bracts and dried out leaf sheaths of restios . Currently there are 21 described restio leafhopper species from two genera, namely Cephalelus and Duospina . Evidence suggests that restio leafhoppers did not co-diversify with restios. Instead they diversified much more recently (1–6 MYA) than restios , which are approximately 15 times more diverse than restio leafhoppers in terms of the number of described species. Thus, many restio species are not exploited by restio leafhoppers, leaving many unfilled niches onto which restio leafhoppers can potentially radiate.
Our study focuses on Cephalelus uncinatus, a broadly distributed restio leafhopper species which completes its entire life cycle on restios. Oviposition occurs on the host-plant and since the eggs have no protective coverings (see e.g. ) they are unable to survive fires (WJA, personal observation). C. uncinatus uses several genera of restios from the Willdenowieae sub tribe, as well as several species in the genus Elegia which belongs to the Restioneae sub tribe . Augustyn et al.  demonstrated experimentally that C. uncinatus from a single site actively chooses its predominant field host and also survives better on it than on unused restio species. As experiments were performed in the absence of predators, the authors suggest that preference is linked to performance through plant chemistry. However, C. uncinatus may also gain protection from predators by choosing restios that serve as good camouflage backgrounds. Previous studies have shown that restio leafhoppers maintain consistent host preference in the presence and absence of interspecific competitors , and provide no evidence that intraspecific competition broadens host preference in restio leafhoppers . These studies suggest that competition has not been important in the diversification of restio leafhoppers, leaving geographic mosaics of phenotypically different hosts as a likely explanation for host-shifts.
Across its distribution range, C. uncinatus uses different host species, and Prendini  suggested that it consists of several ecotypes. He reported consistent, but slight differences in genitalia between populations using Willdenowia incurvata, Mastersiella digitata and Elegia nuda as host-plants (species that differ in culm thickness and leaf sheath colour). Since competition appears to be an unlikely driver of host-shifts in restio leafhoppers [41, 42], we determine whether spatial heterogeneity in plant distributions could have initiated host-shifts in C. uncinatus. If so, we expect to find: (1) that host-plant distributions are smaller than the distribution of C. uncinatus, thus generating a geographic mosaic in the availability of different hosts across the insect species range, and (2) that host preference, colour and shape of restio leafhoppers are locally adapted to their host-plants across this geographic mosaic. 3) Lastly, only finding evidence for divergence between populations using host-plants with spatially non-overlapping distributions would support an allopatric over a sympatric host-shift model (see Fig. 1).
Distribution overlap between C. uncinatus and its hosts
We assessed the spatial heterogeneity of host-plant availability with plant and insect distribution records. We sourced 19 records of C. uncinatus from a MSc thesis  and 14 from an honours thesis . Through our own opportunistic collecting and standardised sampling we obtained 11 other distribution records . All C. uncinatus individuals were either identified by WJA, Davies  or Prendini . Restio distribution records were sourced from the leading expert on the Restionaceae, HP Linder (unpublished data). This data set mostly includes records from the Bolus herbarium and fieldwork conducted by HP Linder.
For the abovementioned restio species, and C. uncinatus, we estimated ranges by calculating convex hulls (polygons based on the outlines of point distributions) around all distribution records using the gConvexhull function in the R package rgeos . We then used the gIntersection function in rgeos to measure the proportion of the C. uncinatus range which overlapped with each restio species. Projections and data conversions were performed using the R packages sp  and rgdal . By means of a bootstrapping procedure, overlap data were analysed by obtaining a mean and upper 95% CI using the one.boot function in the R package simpleboot . An upper 95% CI lower than 1 indicates that host-plants have ranges that don’t cover the total distribution range of C. uncinatus. We reasoned that if none of the host species’ ranges overlap fully with the range of C. uncinatus, then by necessity, populations of C. uncinatus will vary in their host-use.
Local adaptation in preference
The host-plant arrangement at the site level also differed between the three different ecotype pair comparisons. Restio leafhoppers using hosts with low range overlap were collected from either M. digitata at Pringle Bay or W. incurvata at Rondeberg. These sites are 113 km away from each other, and gene flow should be minimal (Fig. 3). C. uncinatus were first collected from W. incurvata at Rondeberg on the 28th of October 2013. On the following day C. uncinatus was collected from M. digitata in Pringle Bay. When the experiment was repeated later, insects were first collected from M. digitata on the 19th of November 2013 at Pringle bay and then from W. incurvata the following day.
For local adaptation of herbivore populations using host species with moderate overlap, C. uncinatus were all collected at De Hoop on stands of either E. nuda or M. spathulata 800 m away from each other (Fig. 3). It is therefore likely that there is some level of gene flow between populations using different host-plants. Collections were made on the 6th and the 8th of January 2014.
To investigate local adaptation in populations with high host species overlap, insects were collected at Pringle Bay from either M. digitata or H. aristatus plants that were in both monospecific and mixed stands. It is therefore likely that, at the site level, there is frequent gene flow between populations using either M. digitata or H. aristatus. Insects were collected on the 6th and the 13th of November 2013.
For all experiments, insects were collected by vacuuming them off plants using a modified leaf blower/shredder. Captured insects were placed singly into clean Eppendorf vials, which were then placed in a cooling box. Later, insects were transferred to a fridge at 10 °C until experiments started the next day (approximately 12 h later) or two days later for the low overlap comparison (36 h). Restio culms of the relevant hosts were collected at the insect collecting sites and placed into distilled water to keep them fresh. Once in the laboratory, culms were cut to the length of 135 mm and kept fresh in a fridge at 10 °C.
We analysed host preference data using separate binomial generalised estimating equations (GEE) with log link functions and exchangeable correlation structures for each of the three experiments (correlation structure choice had little effect on models). The response variable was always a binary choice for Mastersiella (either M. digitata in the low range overlap and high range overlap experiments or M. spathulata in the moderate range overlap experiment). Each model tested for the role of the variables time (after start of experiment), sex and population, in determining host choice (whilst accounting for non-independence of repeat measurements from an individual). For each GEE we back transformed and plotted estimates of means and 95% CIs on a scale ranging from 0 to 1. When 95% CIs did not overlap with a preference of 0.5, a preference for one of the compared restio species was inferred. We also calculated standardised effect sizes (Hedges’ d) for the population variable using GEE estimated means and standard deviations . All GEEs were implemented in the geepack package  in R version 3.2.2.
Body colour is another trait likely to be under divergent selection if the colours of sheaths differ on different restio hosts. More specifically we expect insect body colour to match the colour of host-plant leaf sheaths. This may occur if selection from predation favours restio leafhoppers that colour-match the leaf sheaths of restios on which they occur. Using an Ocean Optics USB4000 spectrometer we recorded spectral colour data of C. uncinatus caught on the five different host-plant species used in the reciprocal restio presentation experiments. Before taking colour measurements, the spectrometer was allowed to heat up for approximately 45 min. Thereafter, light and dark calibrations were performed every 10 min. All insects were measured once, while two repeat measures of each restio sheath were taken and averaged. C. uncinatus sample sizes are shown in Fig. 5. Restio sheath sample sizes were as follows: 18 individuals of W. incurvata, 15 of M. digitata, 12 of H. aristatus, 10 of M. spathulata and 11 of E. nuda.
We modelled the spectral data in tetrahedral colour space of potential visual predators (i.e. birds) of restio leafhoppers using the R package PAVO. Specifically, we modelled colour under the average avian VS model , which assumes that the birds hunting restio leafhoppers are not very UV (ultraviolet) sensitive. Modelling was conducted assuming D65 standard daylight. Using the colour space model, we tested whether each population colour matches the leaf sheaths of their host better than non-hosts. Specifically, we made the same pairwise comparisons as above (low host overlap, moderate host overlap and high host overlap). For each individual insect, we determined the ΔS in just noticeable differences (JNDs) between itself and the closest restio sheath point within the cluster of sheath measurements of each restio species. A lower ΔS indicates a better colour match between insect and sheath. We used two-way repeated measures ANOVAs including the variables, population origin (i.e. which restio species it is from) and restio species identity. As all individual insects were represented by two measurements (i.e. ΔS relative to host and to non-host), we included individual identity as a random variable. Female and male insects were analysed separately. A significant interaction between population origin and plant species identity would be indicative of local adaptation to sheath colour .
Distribution overlap between C. uncinatus and its hosts
Restios are not homogenously distributed across the distribution range of C. uncinatus and there is no restio species that occurs across the entire range of C. uncinatus (C. u in Fig. 2). On average, each host species was only present in 33% of the C. uncinatus range (48% upper 95% CI). Two broadly distributed species, W. teres (Fig. 2, W. i) and E. filacea (Fig. 2, E. fi), occur across 74 and 76% of the distribution of C. uncinatus, but these are both very rarely used hosts (Additional file 1: Table S1). The majority of species (7 out of 10) occurred across less than 50% of the C. uncinatus distribution range. Many of these species only overlapped with a small fraction of the C. uncinatus distribution range (as low as 2%). Some species that occurred in only a small part of the C. uncinatus distribution range were frequently used by C. uncinatus. For example, M. digitata (Fig. 2, M. d) occurs in only 16% of C. uncinatus’ distribution, but is always used when encountered (4/4 times in distribution records). It is also interesting that host species tended to have either North to Southwest or East to Southwest distributions (Fig. 2), so that many restio species overlapped more towards the Southwest. C. uncinatus populations towards the Northern or Eastern extremes of its distribution range therefore encounter very different suites of host-plants.
Divergent host preferences were found in C. uncinatus populations using hosts with low and moderate range overlap but not in populations using hosts with high geographic overlap. Individuals captured on M. digitata at Pringle Bay and W. incurvata at Rondeberg (from hosts with low range overlap, Fig. 3) had strongly divergent host preferences (Wald = 26.51, P < 0.001, Fig. 4a). However, there was still a general preference for W. incurvata. Females and males from W. incurvata showed a strong preference for this species (95% CI lower than 0.5). Females from M. digitata showed a weaker preference for W. incurvata than females from W. incurvata and males had no preference (95% CI overlapped with 0.5). Generally, males had a stronger preference for M. digitata than females (Wald = 15.90, P < 0.001, Fig. 4). Time since the start of the experiment did not influence host preference in this model (Wald = 1.70, P = 0.192).
Individuals captured on either M. spathulata or E. nuda at De Hoop (from hosts with moderate range overlap, Fig. 3) also had strongly divergent host preferences (Wald = 67.23, P < 0.001, Fig. 3). Both sexes consistently preferred the host-plants on which they were found. Preference differences between the sexes for the moderate host overlap comparison were weak (Wald = 4.11, P = 0.043, Fig. 3). Time since the start of the experiment did not have a strong effect on host preference (Wald = 2.87, P = 0.090).
Individuals captured on M. digitata and H. aristatus at Pringle Bay (hosts with high range overlap, Fig. 3) did not have significantly divergent host preferences (Wald = 3.20, P = 0.074, Fig. 4). Time since the start of the experiment influenced host preference (Wald = 5.02, P = 0.025). Specifically, there was a tendency for insects collected from both plants to prefer H. aristatus over time.
Divergence in the body dimensions of C. uncinatus was found in populations using hosts with both low and moderate range overlap, but not in populations using hosts with high range overlap. Insects from the thinner culmed M. digitata at Pringle bay were thinner than those caught from thick culmed W. incurvata at Rondeberg (F = 273.29, df = 1, P < 0.001), and females tended to be thicker than males (F = 55.31, df = 1, P < 0.001) (comparison with populations using hosts with low range overlap, Fig. 5). In this comparison insects from M. digitata also tended to be less stocky (body width divided by elytron length) than those on W. incurvata (F = 42.235, df = 1, P < 0.001), and females tended to be stockier than males (F = 92.185, df = 1, P < 0.001) (Additional file 3: Figure S2).
At De Hoop where populations use plants with moderately overlapping distribution ranges (Fig. 5), insects caught from the thicker culmed M. spathulata were thicker than those caught from the thin culmed E. nuda (F = 52.092, df = 1, P < 0.001). Females were also thicker than males (F = 113.163, df = 1, P < 0.001). Stockiness, however, did not differ between populations (F = 0.612, df = 1, P = 0.436), but females tended to be stockier than males (F = 192.438, P < 0.001) (Additional file 3: Figure S2). In other words, unlike for the low overlap comparison, body width divergence in this pair reflected body size divergence alone.
At Pringle bay (where populations used hosts with high range overlap, Fig. 5), there were no notable differences in any body dimensions between insects using different host-plants. Insects caught from M. digitata and H. aristatus had the same body widths (F = 1.860, df = 1, P = 0.176), and like for all the other comparisons, females were thicker than males (F = 65.160, df = 1, P < 0.001). Similarly, insects using different host-plants did not differ in stockiness (F = 1.254, df = 1, P = 0.265), but females tended to be stockier than males (F = 95.259, df = 1, P < 0.001) (Additional file 3: Figure S2). Complete ANOVA results of body width and stockiness are shown in in Additional file 4: Tables S2 and Additional file 5: Table S3.
Restio sheath colour matching
Unlike most taxa, where allopatric over sympatric divergence is typically invoked as the most important pathway to speciation (reviewed in ); a great deal of emphasis has been placed on sympatric speciation in herbivorous insects [17, 18]. This manuscript clearly demonstrates that the geographic heterogeneity required for allopatric divergence is ubiquitous in this system. Furthermore, it also demonstrates instances of insect local adaptation to different allopatric host-plants. We suggest that like most other taxa, allopatry may be the dominant pathway to ecological speciation in herbivorous insects and that the role of sympatric divergence in herbivorous insects may have been overplayed in the literature. Below, we evaluate the evidence in support of geographic host-plant mosaics and ecotype formation in allopatry as assessed through host preference and morphological divergence.
Here we provide rare empirical evidence suggesting that host-plant mosaics underlie ecologically driven divergence of insects in the CFR. Seven out of the ten known host-plant species, potentially acting as different selective environments (sensu ), have distribution ranges that overlap with less than 50% of the distribution range of C. uncinatus. A similar finding was reported by Kemp et al.  who found that spatially separated herbivorous insect populations (and species) on restios often use a different suite of host-plants because of spatial turnover of host-plants. This pattern is mirrored by herbivorous insects in the tropics , suggesting that it may be common for spatially separated herbivorous insect populations of the same species to experience divergent selection as a result of plant species turnover. While there is also evidence that host-shifts can occur in sympatry [1, 56], host-shifts resulting from geographic turnover of host-plants may be commonplace. Allopatric shift models are also thought to be important in initiating ecological speciation in plants [57–59] and most other taxa . We suggest that divergence of herbivorous insects may not be different from other taxa in terms of the strong role played by allopatry.
Divergence in host preference
Divergence in host preference plays an important role in insect speciation because it can directly result in assortative mating . For this reason, host preference can be viewed as a “magic trait” promoting reproductive isolation (sensu ). However, no studies to our knowledge have explicitly shown that female and male insects prefer the same host-plants (in the context of divergence), which is essential for divergent host preference to result in assortative mating (implied by ). In our study, females and males differed slightly in host preference. Nonetheless, different sexes of the same ecotype always had more similar host preferences than different sexes from different ecotypes (see Fig. 4). Thus, considering that C. uncinatus mates on its host-plants (WJA, personal observation), divergent host preference can result in assortative mating. Host preference divergence was strongest in the comparison where populations used hosts exhibiting moderate range overlap, and absent in the comparison where populations were using host-plants with high range overlap (Fig. 4). This partly suggests that, like in other systems, divergent host preference can be maintained in the presence of moderate gene flow , or that assortative mating via divergent host preference is preventing maladaptive gene flow . Therefore, even if host preference is developmentally plastic, preferences for different host-plants should decrease gene flow, thereby facilitating adaptive divergence in heritable traits.
In this study, we measured traits that are likely under selection from predation. Similar to the preference experiments we cannot rule out the possibility of developmental plasticity or that trait values resulted from non-adaptive processes. Nonetheless, consistent with the findings from preference experiments, we only detected divergent body width in population pairs using host-plants with low and moderate range overlap, suggesting that divergence is adaptive and not due to developmental plasticity. The population pair that used host-plants exhibiting low range overlap also differed in body width relative to body size, while no such scaling difference was found between the moderate overlap population pair. Additionally, larger individuals were always found on plants with thicker culms, suggesting that natural selection has driven body size divergence. Divergent body size could have resulted from predators selecting for smaller (or thinner) insects on thinner culms, allowing them to hide behind culms. Alternatively, large bodied restio leafhoppers may not be able to hold on to thin culms in strong winds and vice versa. Besides camouflage there might, therefore, be a fitness cost to poor mechanical fit (reviewed in ), and therefore divergent selection on body size and width.
Regardless of what drives body size, our findings suggest that there is divergence in body size, which is often regarded as a classic magic trait because of its direct involvement in mate choice . Body size in Bahamas mosquitofish, for example, is involved in mate selection whilst simultaneously under divergent selection from predation . This is thought to play an important role in early stage reproductive isolation between populations of this species . While we are not aware of any reported cases of divergent body size causing sexual isolation between populations of herbivorous insect species, male insects often choose mates based on their body size . Thus, body size in C. uncinatus might be an important trait involved in assortative mating and a promoter of rapid reproductive isolation. Our findings, therefore, warrant further investigation into body size divergence and sexual isolation in C. uncinatus.
Local adaptation in body colour partly reflected preference experiments and body width/size divergence. In bird vision, only C. uncinatus from populations using host-plants with moderate range overlap showed evidence for local adaptation in leaf sheath colour matching. No such evidence was found for population pairs using host-plants with low and high distribution overlap (Fig. 6 left and right panels). Nonetheless, both populations using host-plants with moderate range overlap (Fig. 3 middle panel) matched M. spathulata leaf sheaths best (despite a significant host x population interaction effect which is indicative of local adaptation, Fig. 6 middle panels). In addition, E. nuda females showed more between individual colour variation than females from M. digitata (some E. nuda females are dark red while others are browner) (Fig. 6 middle panels). One explanation for this apparent maladaptation  is that gene flow is occurring from the M. spathulata to the E. nuda ecotype. Nevertheless, some degree of divergence is maintained despite the possibility of gene flow between the M. spathulata and E. nuda populations. Although reduced fitness in alternative environments is not directly related to assortative mating, it can indirectly increase reproductive isolation through immigrant inviability . This process only requires that maladapted immigrants have reduced survivorship prior to mating, which as a by-product results in assortative mating [22, 62]. Considering that predatory birds should be better at detecting immigrants in the moderate distribution range overlap comparison (especially females), immigrant inviability may cause reproductive isolation.
Multiple axes of phenotypic divergence (divergent habitat preference, body size and body colour) driven by host-shifts might facilitate reproductive isolation between C. uncinatus ecotypes. Another divergent selection pressure that might cause immigrant inviability is plant physiology. We previously showed that C. uncinatus adults using H. aristatus have reduced survival on Elegia filacea, a lower ranking host species . It is therefore conceivable that reproductive isolation is, in addition to e.g. predation, maintained by immigrant inviability through physiological trade-offs. Closely related restio leafhopper species tend to use closely related restio species  suggesting that shifts occur between physiologically and morphologically similar species. Therefore, populations may be able to establish on novel host-plants without escape from competition playing an important role in facilitating the shift.
The presence of restio leafhopper ecotypes in allopatry and a geographic mosaic of host-plants makes it very likely that host-shifts occur in allopatry when herbivores encounter novel host communities. While this manuscript does not provide strong evidence to refute the alternative sympatric pathway of host-shifting; previous manuscripts [41, 42], find no evidence for interspecific  and intraspecific  competition in restio leafhoppers, a key element of sympatric host shifting [12–14]. Consequently, we suggest that host-shifts in this system are likely driven by restio leafhoppers adapting allopatrically to novel host-plants. This may occur if restio leafhoppers expand their range into sites with different suites of restios (Fig. 1b). One additional line of evidence supports this pathway; divergent host preference and morphology is evident in population pairs exploiting hosts with little range overlap (i.e. in allopatry), and absent in the high range overlap (i.e. in sympatry) comparison. This line of evidence would be much stronger if experiments were replicated across multiple allopatric and sympatric host-plant pairs , and we suggest that it could be a fertile area for future study, in this and other insect herbivore systems.
Cape Floristic Region
Generalised estimating equation
Just noticeable differences
Million years ago
We are grateful to the NRF South African Biosystematics Initiative and Stellenbosch University for providing the funding that made this study possible. The Western Cape Nature Conservation Board provided permits (0035-AAA007-00082 and 0028-AAA007-00005), and we thank Callum Beattie for permission to work in De Hoop Nature Reserve. We also thank Marion Carrier for help in the field and Corneile Minnaar for advice and conversations on data analyses. Jaatje Goussard provided advice on graphics. Lastly, we thank two anonymous reviewers for helpful comments on the manuscript.
Funding was provided by the NRF South African Biosystematics Initiative (Grant No. 74458 to AGE) and Stellenbosch University (to AGE).
Availability of data and materials
Data are available from the corresponding authors on request.
WJA, BA and AGE conceived the ideas; WJA. conducted all fieldwork, preference trials and statistical analyses. JFV measured morphological traits, and constructed visual models. WJA led the writing with substantial contributions from BA and AGE. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All experiments comply with the current laws of South Africa (ethical clearance is not required for work on insects).
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.
- Bush GL. Sympatric host race formation and speciation in frugivorous flies of the genus Rhagoletis (Diptera, Tephritidae). Evolution (N Y). 1969;23:237–51.Google Scholar
- Bush GL. Mating behavior, host specificity, and the ecological significance of sibling species in frugivorous flies of the genus Rhagoletis (Diptera-Tephritidae). Am Nat. 1969;130:669–72.View ArticleGoogle Scholar
- Ehrlich P, Raven P. Butterflies and plants: a study in coevolution. Evolution (N Y). 1964;18:586–608.Google Scholar
- Imada Y, Kawakita A, Kato M. Allopatric distribution and diversification without niche shift in a bryophyte-feeding basal moth lineage (Lepidoptera: Micropterigidae). Proc R Soc B. 2011;278:3026–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Després L, Jaeger N. Evolution of oviposition strategies and speciation in the globeflower flies Chiastocheta spp. (Anthomyiidae). J Evol Biol. 1999;12:822–31.View ArticleGoogle Scholar
- Nyman T, Vikberg V, Smith DR, Boevé J. How common is ecological speciation in plant-feeding insects? A “higher” Nematinae perspective. BMC Evol Biol. 2010;10:266.View ArticlePubMedPubMed CentralGoogle Scholar
- Gripenberg S, Mayhew PJ, Parnell M, Roslin T. A meta-analysis of preference–performance relationships in phytophagous insects. Ecol Lett. 2010;13:383–93.View ArticlePubMedGoogle Scholar
- Nosil P, Crespi BJ. Experimental evidence that predation promotes divergence in adaptive radiation. Proc Natl Acad Sci U S A. 2006;103:9090–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Brower L. Bird predation and foodplant specificity in closely related procryptic insects. Am Nat. 1958;92:183–7.View ArticleGoogle Scholar
- Winkler IS, Mitter C. The phylogenetic dimension of insect-plant interactions: a review of recent evidence. In: Tilmon KJ, editor. Spec. speciation, Radiat. Evol. Biol. Herbiv. insects. Berkeley: University of California Press; 2008. p. 240–63.
- Wilson DS, Turelli M. Stable underdominance and the evolutionary invasion of empty niches. Am Nat. 1986;127:835.View ArticleGoogle Scholar
- Rosenzweig ML. Competitive speciation. Biol J Linn Soc. 1978;10:275–89.View ArticleGoogle Scholar
- Feder J, Reynolds K, Go W, Wang E. Intra-and interspecific competition and host race formation in the apple maggot fly, Rhagoletis pomonella (Diptera: Tephritidae). Oecologia. 1995;101:416–25.View ArticleGoogle Scholar
- Rausher MD. The evolution of habitat preference in subdivided populations. Evolution (N Y). 1984;38:596–608.Google Scholar
- Rice WR. Disruptive selection on habitat preference and the evolution of reproductive isolation: an exploratory experiment. Evolution (N Y). 1984;38:1251–60.Google Scholar
- Thoday JM, Gibson JB. Isolation by Disruptive Selection. Nature. 1962;193:1164–6.View ArticlePubMedGoogle Scholar
- Drès M, Mallet J. Host races in plant–feeding insects and their importance in sympatric speciation. Phil Trans Royal Soc Lon B. 2002;357:471–92.View ArticleGoogle Scholar
- Berlocher SH, Feder JL. Sympatric speciation in phytophagous insects: moving beyond controversy? Annu Rev Entomol. 2002;47:773–815.View ArticlePubMedGoogle Scholar
- Gurevitch J, Morrow LL, Wallace A, Walsh JS. A meta-analysis of competition in field experiments. Am Nat. 1992;140:539.View ArticleGoogle Scholar
- Nosil P, Sandoval C, Crespi B. The evolution of host preference in allopatric vs. parapatric populations of Timema cristinae walking sticks. J Evol Biol. 2006;19:929–42.View ArticlePubMedGoogle Scholar
- Feder J, Opp S, Wlazlo B, Reynolds K, Go W, Spisak S. Host fidelity is an effective premating barrier between sympatric races of the apple maggot fly. Proc Natl Acad Sci U S A. 1994;91:7990–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Nosil P, Vines TH, Funk DJ. Reproductive isolation caused by natural selection against immigrants from divergent habitats. Evolution (N Y). 2005;59:705–19.Google Scholar
- Ellis AG, Verboom A, van der Niet T, Johnson SD, Linder HP. Speciation and extinction in the greater Cape Floristic Region. In: Allsopp N, Colville JF, Verboom GA, editors. Fynbos Ecol. Evol. Conserv. A megadiverse Reg. New York: Oxfort University Press; 2014. p. 119–41.View ArticleGoogle Scholar
- Mittermeier RA, Myers N, Thomsen JB, da Fonseca GAB, Olivieri S. Biodiversity hotspots and major tropical wilderness areas: approaches to setting conservation priorities. Conserv Biol. 1998;12:516–20.View ArticleGoogle Scholar
- Linder HP. The radiation of the Cape flora, southern Africa. Biol Rev Camb Philos Soc. 2003;78:597–638.View ArticlePubMedGoogle Scholar
- Cowling R, Rundel P. Plant diversity in Mediterranean-climate regions. Trends Ecol Evol. 1996;11:362–6.View ArticlePubMedGoogle Scholar
- Slingsby J, Verboom G. Phylogenetic relatedness limits co-occurrence at fine spatial scales: evidence from the Schoenoid sedges (Cyperaceae: Schoeneae) of the Cape Floristic Region, South. Am Nat. 2006;168:14–27.PubMedGoogle Scholar
- Schnitzler J, Barraclough TG, Boatwright JS, Goldblatt P, Manning JC, Powell MP, et al. Causes of plant diversification in the Cape biodiversity hotspot of South Africa. Syst Biol. 2011;60:343–57.View ArticlePubMedGoogle Scholar
- van der Niet T, Johnson SD. Patterns of plant speciation in the Cape floristic region. Mol Phylogenet Evol. 2009;51:85–93.View ArticlePubMedGoogle Scholar
- Kruger FJ, Bigalke RC. Fire in Fynbos. In: de Booysen PV, Tainton NM, editors. Ecol. Eff. Fire South African ecosyst. Berlin, Heidelberg: Springer Berlin Heidelberg; 1984. p. 67–114.View ArticleGoogle Scholar
- Pryke JS, Samways MJ. Differential resilience of invertebrates to fire. Austral Ecol. 2012;37:460–9.View ArticleGoogle Scholar
- Thuiller W, Slingsby JA, Privett SDJ, Cowling RM. Stochastic species turnover and stable coexistence in a species-rich, fire-prone plant community. PLoS One. 2007;2:e938.View ArticlePubMedPubMed CentralGoogle Scholar
- Kassen R. The experimental evolution of specialists, generalists, and the maintenance of diversity. J Evol Biol. 2002;15:173–90.View ArticleGoogle Scholar
- Linder HP, Eldenäs P, Briggs BG. Contrasting patterns of radiation in African and Australian Restionaceae. Evolution (N Y). 2003;57:2688–702.Google Scholar
- Augustyn WJ, Anderson B, Stiller M, Ellis AG. Specialised host-use and phenophase tracking in restio leafhoppers (Cicadellidae: Cephalelini) in the Cape Floristic Region. J Insect Conserv. 2013;17:1267–74.View ArticleGoogle Scholar
- Kemp JE, Ellis AG. Significant local-scale plant-insect species richness relationship independent of abiotic effects in the temperate Cape Floristic Region biodiversity hotspot. PLoS One. 2017;12:e0168033.View ArticlePubMedPubMed CentralGoogle Scholar
- Osborn H. An interesting case of mimicry. Psyche (Stuttg). 1903;9:327.Google Scholar
- Prendini L. Two new host-restricted restioid leafhoppers of the genus Cephalelus Percheron (Cicadellidae: Cephalelelini), with descriptions of the females of C. brevipilus Davies, C. daviesi Davies and C. rawsonia Davies. African Entomol. 1997;5:273–81.Google Scholar
- Wiese TJ. The evolution of fynbos-endemic Cephalelini leafhoppers specialising on Restionaceae. MSc Thesis University of Stellenbosch (Stellenbosch, South Africa), Botany and Zoology department; 2014.
- Sandoval C. Persistence of a walking-stick population (Phasmatoptera: Timematodea) after a wildfire. Southwest Nat. 2000;45:123–7.View ArticleGoogle Scholar
- Augustyn WJ, Anderson B, Ellis AG. Experimental evidence for fundamental, and not realised, niche partitioning in a plant-herbivore community interaction network. J Anim Ecol. 2016;85:994–1003.View ArticlePubMedGoogle Scholar
- Augustyn WJ. Evolutionary drivers of temporal and spatial host use patterns in restio leafhoppers Cephalelini (Cicadellidae). PhD Thesis University of Stellenbosch (Stellenbosch, South Africa), Botany and Zoology department; 2015.
- Davies DM. A taxonomic description of the leafhopper fauna (Homoptera: Cicadellidae) in association with Restionaceae. MSc Thesis University of Stellenbosch (Stellenbosch, South Africa) Department of Entomology; 1986.
- Prendini L. Parallel cladogenesis, coevolution, and insect-plant associations: evidence from leafhoppers (Cicadellidae, Cephalelini) and Restionaceae. BSc Hons. Thesis University of Cape Town (Cape Town, South Africa), Zoology department; 1995.
- Bivand R, Rundel C, Pebesma E, Stuetz R, Hufthammer KO. rgeos: interface to geometry engine - open source (GEOS). 2016. Available from: https://cran.r-project.org/package=rgeos.
- Pebesma E, Bivand R, Rowlingson B, Gomez-Rubio V, Hijmans R, Sumner M, et al. sp: classes and methods for spatial data. 2016.Google Scholar
- Bivant R, Keitt T, Rowlingson B. rgdal: bindings for the geospatial data abstraction library. 2016. Available from: https://cran.r-project.org/package=rgdal.Google Scholar
- Peng RD. simpleboot: simple bootstrap routines. 2008. Available from: https://cran.r-project.org/package=simpleboot.Google Scholar
- Hedges LV. Distribution theory for glass’s estimator of effect size and related estimators. J Educ Behav Stat. 1981;6:107–28.View ArticleGoogle Scholar
- Højsgaard S. geepack: generalized estimating equation package. 2016. Available from: https://cran.r-project.org/package=geepack.Google Scholar
- Endler JA, Mielke PW. Comparing entire colour patterns as birds see them. Biol J Linn Soc. 2005;86:405–31.View ArticleGoogle Scholar
- Kawecki TJ, Ebert D. Conceptual issues in local adaptation. Ecol Lett. 2004;7:1225–41.View ArticleGoogle Scholar
- Bolnick DI, Fitzpatrick BM. Sympatric speciation: models and empirical evidence. Annu Rev Ecol Evol Syst. 2007;38:459–87.View ArticleGoogle Scholar
- Kemp JE, Evans DM, Augustyn WJ, Ellis AG. Invariant antagonistic network structure despite high spatial and temporal turnover of species and their interactions in a biodiversity hotspot. Ecography (Cop). 2017;40:001–10.View ArticleGoogle Scholar
- Novotny V. Beta diversity of plant-insect food webs in tropical forests: a conceptual framework. Insect Conserv Divers. 2009;2:5–9.View ArticleGoogle Scholar
- Feder JL, Chilcote A, Bush GL. Genetic differentiation between sympatric host races of the apple maggot fly Rhagoletis pomonella. Nature. 1988;336:61–4.View ArticleGoogle Scholar
- Ellis AG, Weis A, Gaut B. Evolutionary radiation of “stone plants” in the genus Argyroderma (Aizoaceae): unraveling the effects of landscape, habitat, and flowering time. Evolution (N Y). 2006;60:39–55.
- Grant V, Grant KA. The process of differentiation. Flower pollinat. Phlox Fam. 1st ed. New York, London: Columbia University Press; 1965. p. 162–5.Google Scholar
- Johnson SD. Pollinator-driven speciation in plants. Ecol. Evol. flowers. 2006. p. 295–230.
- Rundle HD, Nosil P. Ecological speciation. Ecol Lett. 2005;8:336–52.View ArticleGoogle Scholar
- Via S. Reproductive isolation between sympatric races of Pea aphids. I. Gene flow restriction and habitat choice. Evolution (N Y). 1999;53:1446–57.Google Scholar
- Servedio MR, Van Doorn GS, Kopp M, Frame AM, Nosil P. Magic traits in speciation: “magic” but not rare? Trends Ecol Evol. 2011;26:389–97.View ArticlePubMedGoogle Scholar
- Rice WR. Disruptive selection on habitat preference and the evolution of reproductive isolation: a simulation study. Evolution (N Y). 1984;38:1251–60.Google Scholar
- Bolnick DI, Snowberg LK, Patenia C, Stutz WE, Ingram T, Lau OL. Phenotype-dependent native habitat preference facilitates divergence between parapatric lake and stream stickleback. Evolution (N Y). 2009;63:2004–16.Google Scholar
- Matsubayashi KW, Ohshima I, Nosil P. Ecological speciation in phytophagous insects. Entomol Exp Appl. 2010;134:1–27.View ArticleGoogle Scholar
- Langerhans RB, Gifford ME, Joseph EO. Ecological speciation in Gabusia fishes. Evolution (N Y). 2007;61:2056–74.Google Scholar
- Bonduriansky R. The evolution of male mate choice in insects: a synthesis of ideas and evidence. Biol Rev. 2001;76:305–39.View ArticlePubMedGoogle Scholar
- Crespi BJ. The evolution of maladaptation. Heredity (Edinb). 2000;84:623–9.View ArticleGoogle Scholar
- Janz N, Nylin S. The oscillation hypothesis of host-plant range and speciation. In: Tilmon KJ, editor. Spec. Speciation, Radiat. Evol. Biol. Herbiv. Insects. Berkeley: University of California Press; 2008. p. 203–15.Google Scholar