Divergent host plant adaptation and reproductive isolation between ecotypes of Timema cristinae walking sticks

Theoretical and empirical studies have demonstrated that divergent natural selection can promote the evolution of reproductive isolation. Three unresolved questions concern the types of reproductive barriers involved, the role of geography, and the factors determining the extent of progress toward complete speciation. Here I synthesize studies of Timema cristinae host plant ecotypes to address these issues. The approach is to compare the magnitude of multiple reproductive barriers among different ecological and geographic scenarios, where pairs of populations within each scenario are the unit of replication. Application of this approach to T. cristinae revealed that divergent host adaptation can promote the evolution of diverse reproductive barriers, including those that are not inherently ecological. Gene flow in parapatry tended to constrain divergence, with the notable exception of the reinforcement of sexual isolation. Thus, geography affected progress toward speciation but did not influence all reproductive barriers in the same way. Studies of any single pair of taxa often capture only certain stages of the speciation process. For example, reproductive isolation between T. cristinae ecotypes is incomplete, and so only the stages before the completion of speciation have been examined. Studies of more divergent taxa within the genus are required to determine the factors that complete speciation.     

The evolution of host preference in allopatric vs. parapatric populations of Timema cristinae walking-sticks

Divergent habitat preferences can contribute to speciation, as has been observed for host-plant preferences in phytophagous insects. Geographic variation in host preference can provide insight into the causes of preference evolution. For example, selection against maladaptive host-switching occurs only when multiple hosts are available in the local environment and can result in greater divergence in regions with multiple vs. a single host. Conversely, costs of finding a suitable host can select for preference even in populations using a single host. Some populations of Timema cristinae occur in regions with only one host-plant species present (in allopatry, surrounded by unsuitable hosts) whereas others occur in regions with two host-plant species adjacent to one another (in parapatry). Here, we use host choice and reciprocal-rearing experiments to document genetic divergence in host preference among 33 populations of T. cristinae. Populations feeding on Ceanothus exhibited a stronger preference for Ceanothus than did populations feeding on Adenostoma. Both allopatric and parapatric pairs of populations using the different hosts exhibited divergent host preferences, but the degree of divergence tended to be greater between allopatric pairs. Thus, gene flow between parapatric populations apparently constrains divergence. Host preferences led to levels of premating isolation between populations using alternate hosts that were comparable in magnitude to previously documented premating isolation caused by natural and sexual selection against migrants between hosts. Our findings demonstrate how gene flow and different forms of selection interact to determine the magnitude of reproductive isolation observed in nature.     

Natural selection and divergence in mate preference during speciation

Sexual isolation can evolve due to natural selection against hybrids (reinforcement). However, many different forms of hybrid dysfunction, and selective processes that do not involve hybrids, can contribute to the evolution of sexual isolation. Here we review how different selective processes affect the evolution of sexual isolation, describe approaches for distinguishing among them, and assess how they contribute to variation in sexual isolation among populations of Timema cristinae stick-insects. Pairs of allopatric populations of T. cristinae living on different host-plant species exhibit greater sexual isolation than those on the same host, indicating that some sexual isolation has evolved due to host adaptation. Sexual isolation is strongest in regions where populations on different hosts are in geographic contact, a pattern of reproductive character displacement that is indicative of reinforcement. Ecological costs to hybridization do occur but traits under ecological selection (predation) do not co-vary strongly with the probability of between-population mating such that selection on ecological traits is not predicted to produce a strong correlated evolutionary response in mate preference. Moreover, F1 hybrid egg inviability is lacking and the factors contributing to reproductive character displacement require further study. Finally, we show that sexual isolation involves, at least in part, olfactory communication. Our results illustrate how understanding of the evolution of sexual isolation can be enhanced by isolating the roles of diverse ecological and evolutionary processes.     

Evolutionary history of aphid-plant associations and their role in aphid diversification

Aphids are intimately linked with their host plants that constitute their only food resource and habitat, and thus impose considerable selective pressure on their evolution. It is therefore commonly assumed that host plants have greatly influenced the diversification of aphids. Here, we review what is known about the role of host plant association on aphid speciation by examining both macroevolutionary and population-level studies. Phylogenetic studies conducted at different taxonomic levels show that, as in many phytophagous insect groups, the radiation of angiosperms has probably favoured the major Tertiary diversification of aphids. These studies also highlight many aphid lineages constrained to sets of related host plants, suggesting strong evolutionary commitment in aphids’ host plant choice, but they fail to document cospeciation events between aphid and host lineages. Instead, phylogenies of several aphid genera reveal that divergence events are often accompanied by host shifts, and suggest, without constituting a formal demonstration, that aphid speciation could be a consequence of adaptation to new hosts. Experimental and field studies below the species level support reproductive isolation between host races as partly due to divergent selection by their host plants. Selected traits are mainly feeding performances and life cycle adaptations to plant phenology. Combined with behavioural preference for favourable host species, these divergent adaptations can induce pre- and post-zygotic barriers between host-specialized aphid populations. However, the hypothesis of host-driven speciation is seldom tested formally and must be weighed against overlooked explanations involving geographic isolation and non-ecological reproductive barriers in the process of speciation.     

Host plant adaptation in Drosophila mettleri populations

The process of local adaptation creates diversity among allopatric populations, and may eventually lead to speciation. Plant-feeding insect populations that specialize on different host species provide an excellent opportunity to evaluate the causes of ecological specialization and the subsequent consequences for diversity. In this study, we used geographically separated Drosophila mettleri populations that specialize on different host cacti to examine oviposition preference for and larval performance on an array of natural and non-natural hosts (eight total). We found evidence of local adaptation in performance on saguaro cactus (Carnegiea gigantea) for populations that are typically associated with this host, and to chemically divergent prickly pear species (Opuntia spp.) in a genetically isolated population on Santa Catalina Island. Moreover, each population exhibited reduced performance on the alternative host. This finding is consistent with trade-offs associated with adaptation to these chemically divergent hosts, although we also discuss alternative explanations for this pattern. For oviposition preference, Santa Catalina Island flies were more likely to oviposit on some prickly pear species, but all populations readily laid eggs on saguaro. Experiments with non-natural hosts suggest that factors such as ecological opportunity may play a more important role than host plant chemistry in explaining the lack of natural associations with some hosts.     

Reproductive isolation driven by the combined effects of ecological adaptation and reinforcement

Recent years have seen a resurgence of interest in the process of speciation but few studies have elucidated the mechanisms either driving or constraining the evolution of reproductive isolation. In theory, the direct effects of reinforcing selection for increased mating discrimination where interbreeding produces hybrid offspring with low fitness and the indirect effects of adaptation to different environments can both promote speciation. Conversely, high levels of homogenizing gene flow can counteract the forces of selection. We demonstrate the opposing effects of reinforcing selection and gene flow in Timema cristinae walking-stick insects. The magnitude of female mating discrimination against males from other populations is greatest when migration rates between populations adapted to alternate host plants are high enough to allow the evolution of reinforcement, but low enough to prevent gene flow from eroding adaptive divergence in mate choice. Moreover, reproductive isolation is strongest under the combined effects of reinforcement and adaptation to alternate host plants. Our findings demonstrate the joint effects of reinforcement, ecological adaptation and gene flow on progress towards speciation in the wild.     

Use it or lose it: reproductive implications of ecological specialization in a haematophagous ectoparasite

Using experimentally induced disruptive selection, we tested two hypotheses regarding the evolution of specialization in parasites. The 'trade-off' hypothesis suggests that adaptation to a specific host may come at the expense of a reduced performance when exploiting another host. The alternative 'relaxed selection' hypothesis suggests that the ability to exploit a given host would deteriorate when becoming obsolete. Three replicate populations of a flea Xenopsylla ramesis were maintained on each of two rodent hosts, Meriones crassus and Dipodillus dasyurus, for nine generations. Fleas maintained on a specific host species for a few generations substantially decreased their reproductive performance when transferred to an alternative host species, whereas they generally did not increase their performance on their maintenance host. The results support the 'relaxed selection' hypothesis of the evolution of ecological specialization in haematophagous ectoparasites, while suggesting that trade-offs are unlikely drivers of specialization. Further work is needed to study the extent by which the observed specializations are based on epigenetic or genetic modifications.     

The genetic architecture of a niche: variation and covariation in host use traits in the Colorado potato beetle

The genetic basis of host plant use by phytophagous insects can provide insight into the evolution of ecological niches, especially phenomena such as specialization and phylogenetic conservatism. We carried out a quantitative genetic analysis of multiple host use traits, estimated on five species of host plants, in the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae). Mean values of all characters varied among host plants, providing evidence that adaptation to plants may require evolution of both behavioral (preference) and post-ingestive physiological (performance) characteristics. Significant additive genetic variation was detected for several characters on several hosts, but not in the capacity to use the two major hosts, a pattern that might be caused by directional selection. No negative genetic correlations across hosts were detected for any 'performance' traits, i.e. we found no evidence of trade-offs in fitness on different plants. Larval consumption was positively genetically correlated across host plants, suggesting that diet generalization might evolve as a distinct trait, rather than by independent evolution of feeding responses to each plant species, but several other traits did not show this pattern. We explored genetic correlations among traits expressed on a given plant species, in a first effort to shed light on the number of independent traits that may evolve in response to selection for host-plant utilization. Most traits were not correlated with each other, implying that adaptation to a novel potential host could be a complex, multidimensional 'character' that might constrain adaptation and contribute to the pronounced ecological specialization and the phylogenetic niche conservatism that characterize many clades of phytophagous insects.     

Evolution of host specificity drives reproductive isolation among RNA viruses

Ecological speciation hypotheses claim that assortative mating evolves as a consequence of divergent natural selection for ecologically important traits. Reproductive isolation is expected to be particularly likely to evolve by this mechanism in species such as phytophagous insects that mate in the habitats in which they eat. We tested this expectation by monitoring the evolution of reproductive isolation in laboratory populations of an RNA virus that undergoes genetic exchange only when multiple virus genotypes coinfect the same host. We subjected four populations of the RNA bacteriophage phi6 to 150 generations of natural selection on a novel host. Although there was no direct selection acting on host range in our experiment, three of the four populations lost the ability to infect one or more alternative hosts. In the most extreme case, one of the populations evolved a host range that does not contain any of the hosts infectible by the wild-type phi6. Whole genome sequencing confirmed that the resulting reproductive isolation was due to a single nucleotide change, highlighting the ease with which an emerging RNA virus can decouple its evolutionary fate from that of its ancestor. Our results uniquely demonstrate the evolution of reproductive isolation in allopatric experimental populations. Furthermore, our data confirm the biological credibility of simple "no-gene" mechanisms of assortative mating, in which this trait arises as a pleiotropic effect of genes responsible for ecological adaptation.     

DOES GENE FLOW CONSTRAIN ADAPTIVE DIVERGENCE OR VICE VERSA? A TEST USING ECOMORPHOLOGY AND SEXUAL ISOLATION IN TIMEMA CRISTINAE WALKING‐STICKS

Population differentiation often reflects a balance between divergent natural selection and the opportunity for homogenizing gene flow to erode the effects of selection. However, during ecological speciation, trait divergence results in reproductive isolation and becomes a cause, rather than a consequence, of reductions in gene flow. To assess both the causes and the reproductive consequences of morphological differentiation, we examined morphological divergence and sexual isolation among 17 populations of Timema cristinae walking-sticks. Individuals from populations adapted to using Adenostoma as a host plant tended to exhibit smaller overall body size, wide heads, and short legs relative to individuals using Ceonothus as a host. However, there was also significant variation in morphology among populations within host-plant species. Mean trait values for each single population could be reliably predicted based upon host-plant used and the potential for homogenizing gene flow, inferred from the size of the neighboring population using the alternate host and mitochondrial DNA estimates of gene flow. Morphology did not influence the probability of copulation in between-population mating trials. Thus, morphological divergence is facilitated by reductions in gene flow, but does not cause reductions in gene flow via the evolution of sexual isolation. Combined with rearing data indicating that size and shape have a partial genetic basis, evidence for parallel origins of the host-associated forms, and inferences from functional morphology, these results indicate that morphological divergence in T. cristinae reflects a balance between the effects of host-specific natural selection and gene flow. Our findings illustrate how data on mating preferences can help determine the causal associations between trait divergence and levels of gene flow.     

Oviposition preference and larval performance within a diverging lineage of lycaenid butterflies

Abstract. 1. The butterfly genus Mitoura in Northern California includes three nominal species associated with four host plants having parapatric or interdigitated ranges. Genetic analyses have shown the taxa to be very closely related, and adults from all host backgrounds will mate and produce viable offspring in the laboratory. Oviposition preference and larval performance were investigated with the aim of testing the hypothesis that variation in these traits can exist in a system in which non‐ecological barriers to gene flow (i.e. geographic barriers and genetic incompatibilities) appear to be minimal. 2. Females were sampled from 12 locations throughout Northern California, including sympatric and parapatric populations associated with the four different host‐plant species. Oviposition preference was assayed by confining wild‐caught females with branches of all four host species and counting the number of eggs laid on each. Offspring were reared on the same host species and two measures of larval success were taken: per cent survival and pupal weight. 3. For populations associated with one of the hosts, incense cedar, the preference–performance relationship is simple: the host that females chose is the plant which results in the highest pupal weights for offspring. The preference–performance relationship for populations associated with the other hosts is more complex and may reflect different levels of local adaptation. The variation in preference and performance reported here suggests that these traits can evolve when non‐ecological barriers to gene flow are low, and that differences in these traits may be important for the evolution of reproductive isolation within Mitoura.     

Mate preference across the speciation continuum in a clade of mimetic butterflies

Premating behavioral isolation is increasingly recognized as an important part of ecological speciation, where divergent natural selection causes the evolution of reproductive barriers. A number of studies have now demonstrated that traits under divergent natural selection also affect mate preferences. However, studies of single species pairs only capture a snapshot of the speciation process, making it difficult to assess the role of mate preferences throughout the entire process. Heliconius butterflies are well known for their brightly colored mimetic warning patterns, and previous studies have shown that these patterns are also used as mate recognition cues. Here, we present mate preference data for four pairs of sister taxa, representing different stages of divergence, which together allow us to compare diverging mate preferences across the continuum of Heliconius speciation. Using a novel Bayesian approach, our results support a model of ecological speciation in which strong premating isolation arises early, but continues to increase throughout the continuum from polymorphic populations through to "good," sympatric ecologically divergent species.     

Experimental evidence of host race formation in Mitoura butterflies (Lepidoptera: Lycaenidae)

A population of herbivorous insects that shifts to a novel host can experience selection pressures that result in adaptation to the new resource. Host race formation, considered an early stage of the speciation process, may result. The current study investigates host shifts and variation in traits potentially involved in the evolution of reproductive isolation among populations of the juniper hairstreak butterfly, Mitoura gryneus. Mitoura are closely associated with their host trees (Cupressaceae) and exhibit host plant fidelity: in addition to larval development and oviposition, host trees support male leks and mating. Female oviposition preference for the natal host, and differential fitness of larvae when reared on natal versus alternate hosts, are indications that specialization and local adaptation to the natal host plant are occurring. Populations with single host plant associations (Juniperus ashei, J. pinchotii and J. virginiana) as well as populations with multiple hosts (both J. ashei and J. pinchotii) were examined. Concordance between female preference and larval performance was found for J. ashei‐associated populations. Population‐level variation in the patterns of female preference and larval performance, both within and among host associations, may reflect differences in the timing and direction of colonization of hosts. For a single nominal species that otherwise exhibits no morphological or phenological differences, the experimental assessment of specialization and host fidelity in M. gryneus provides strong support for the hypothesis of ongoing host race formation in these butterflies.     

A test of adaptive divergence in a newly discovered host association of the soapberry bug Jadera haematoloma on Mexican buckeye,Ungnadia speciosa

Host races represent an important step in the speciation process of phytophagous insects as they reflect the maintenance of genetically divergent host-associated populations in the face of appreciable gene flow. The red-shouldered soapberry bug, Jadera haematoloma (Herrich-Schäffer) (Hemiptera: Rhopalidae), is an oligophagous seed predator with a history of host race evolution on plant associations in the (soapberry) family Sapindaceae. Soapberry bugs are a model group for understanding rapid ecological adaptation to their hosts, and hence good candidates for investigating evolutionary divergence in host associations over short timescales. Here, we describe the recent discovery of Mexican buckeye, Ungnadia speciosa Endl., as a host of J. haematoloma in a region of the Chihuahuan desert including west Texas and southeastern New Mexico, USA. This host differs from J. haematoloma’s previously recorded hosts in the Sapindaceae in seed chemistry, ecology, and phylogeny. The tendency toward rapid, host-associated adaptations by populations of J. haematoloma and the unique biology of the newly discovered Ungnadia host create the opportunity for potential host race formation, as it overlaps geographically with two previously recorded host plants in this region – the native western soapberry tree, Sapindus saponaria var. drummondii (Hook & Arn.), and the non-native goldenrain tree, Koelreuteria paniculata Laxm. We explore the possibility of host race formation on Ungnadia-associated insects by testing for host-associated differentiation in morphology and feeding behaviors. We find evidence of differentiation in the length of the mouth parts, which is an ecologically relevant feeding trait between host plant species with larger or smaller seed capsules. This divergence is maintained in the face of potential gene flow by reproductive isolation in the form of habitat isolation, which we detect in host preference trials. Together, our results demonstrate that soapberry bugs associated with this newly discovered host exhibit morphological and behavioral traits consistent with host race formation, but additional work is required to confirm its state along the speciation continuum.     

Host‐associated divergence in sympatric host races of the leaf beetle Lochmaea capreae: implications for local adaptation and reproductive isolation

Ecological specialization is widely recognized as a major determinant of the emergence and maintenance of biodiversity. We studied two critical facets of specialization – local adaptation and habitat choice – in the host races of the leaf beetle Lochmaea capreae on willow and birch. Our results revealed that there is asymmetric disruptive selection for host use traits, and host races achieved different adaptive sets of life history traits through association with their host plant. Beetles from each host race exhibited food and oviposition preference for their own host plant. Reciprocal transplant displayed significant variation in host acceptance and performance: all families from the willow race rejected the alternative host plant before initiation of feeding and all died on this host plant. By contrast, all families from the birch race accepted willow for feeding, but they consumed less and performed less well. Intriguingly, families that performed well on birch also performed well on willow, suggesting positive genetic correlation rather than genetic trade‐offs. Our results suggest that the major proximal determinant of host specialization in the willow race is the behavioural acceptance of a plant rather than the toxicity of the food resource. However, in the birch race a combination of behavioural host acceptance and performance may play a role in specialization. Our study sheds light on the mechanisms by which divergent host adaptation might influence the evolution of reproductive isolation between herbivorous populations.     

Lack of evidence for reproductive isolation among ecologically specialised lycaenid butterflies

Abstract 1. The evolution of reproductive isolation between recently diverged or incipient species is a critical component of speciation and a major focus of speciation models. In phytophagous insects, host plant fidelity (the habit of mating and ovipositing on a single host species) can contribute to assortative mating and reproductive isolation between populations adapting to alternative hosts. The potential role of host plant fidelity in the evolution of reproductive isolation was examined in a pair of North American blue butterfly species, Lycaeides idas and L. melissa .
2. These species are morphologically distinct and populations of each species utilise different host plants; however they share 410 bp haplotypes of the mitochondrial cytochrome oxidase subunit I (COI) gene, indicating recent divergence.
3. Some populations using native hosts exhibited strong fidelity for their natal host plant over the hosts used by nearby populations. Because these butterflies mate on or near the host plant, the development of strong host fidelity may create reproductive isolation among populations on different hosts and restrict gene flow.
4. Tests of population differentiation using allozyme allele frequency data did not provide convincing evidence of restricted gene flow among populations. Based on morphological differences, observed ecological specialisation, and the sharing of genetic markers, these butterflies appear to be undergoing adaptive radiation driven at least partially by host shifts. Neutral genetic markers may fail to detect the effects of very recent host shifts in these populations due to gene flow and/or the recency of divergence and shared ancestral polymorphism.     

Sexual and local divergence in host exploitation in the marine herbivore Idotea baltica (Isopoda)

We studied sexual and habitat-specific variation in herbivores' host exploitation patterns and tested for a hypothesis on sex-specific adaptation in replicated landscapes. The hypothesis has not to our knowledge been previously tested either in a marine environment or for the host exploitation traits of an herbivore. The hypothesis states, first, that populations may show different adaptations due to spatially differing selective environments, and second, that males and females that differ ecologically (e.g. in reproductive behavior, in host use patterns) will respond to selective environments in distinct ways. We investigated possible differences in the host exploitation patterns of the marine generalist crustacean grazer Idotea baltica between the sexes or among populations originating from different habitats: plant assemblages dominated by either a brown alga or an angiosperm species. We determined the preferences for both the structure and nutritive quality of five common host species and compared these to the actual performance on the hosts. We found that performance on the hosts differed between the sexes and differently so in the two habitat types, supporting the hypothesis of sex-specific adaptation between distinct selective environments. Habitat-specific preference for the structural host further supported ecological divergence between habitat types. The different approaches for testing host exploitation resulted in a very different rank order of the hosts, indicating the complexity of the factors involved. The herbivore differentiates for the chemical quality of the hosts, but this alone does not explain host exploitation; rather it supports the theory of enemy-free space as a determinant of host choice. We also discuss the chemical characteristics of the host associated with the exploitation patterns found.     

Diversification across a heterogeneous landscape

Adaptation to contrasting environments across a heterogeneous landscape favors the formation of ecotypes by promoting ecological divergence. Patterns of fitness variation in the field can show whether natural selection drives local adaptation and ecotype formation. However, to demonstrate a link between ecological divergence and speciation, local adaptation must have consequences for reproductive isolation. Using contrasting ecotypes of an Australian wildflower, Senecio lautus in common garden experiments, hybridization experiments, and reciprocal transplants, we assessed how the environment shapes patterns of adaptation and the consequences of adaptive divergence for reproductive isolation. Local adaptation was strong between ecotypes, but weaker between populations of the same ecotype. F1 hybrids exhibited heterosis, but crosses involving one native parent performed better than those with two foreign parents. In a common garden experiment, F2 hybrids exhibited reduced fitness compared to parentals and F1 hybrids, suggesting that few genetic incompatibilities have accumulated between populations adapted to contrasting environments. Our results show how ecological differences across the landscape have created complex patterns of local adaptation and reproductive isolation, suggesting that divergent natural selection has played a fundamental role in the early stages of species diversification.     

Phylogenetic evidence for multiple sympatric ecological diversification in a marine snail

Parallel speciation can occur when traits determining reproductive isolation evolve independently in different populations that experience a similar range of environments. However, a common problem in studies of parallel evolution is to distinguish this hypothesis from an alternative one in which different ecotypes arose only once in allopatry and now share a sympatric scenario with substantial gene flow between them. Here we show that the combination of a phylogenetic approach with life-history data is able to disentangle both hypotheses in the case of the intertidal marine snail Littorina saxatilis on the rocky shores of Galicia in northwestern Spain. In this system, numerous phenotypic and genetic differences have evolved between two sympatric ecotypes spanning a sharp ecological gradient, and as aside effect of the former have produced partial reproductive isolation. A mitochondrial phylogeny of these populations strongly suggests that the two sympatric ecotypes have originated independently several times. Building upon earlier work demonstrating size-based assortative mating as the main contributor to reproductive isolation among ecotypes, our analysis provides strong evidence that divergent selection across a sharp ecological gradient promoted the parallel divergence of body size and shape between two sympatric ecotypes. Thus, divergent selection occurring independently in different populations has produced the marine equivalent of host races, which may represent the first step in speciation.     

WIDESPREAD HOST‐DEPENDENT HYBRID UNFITNESS IN THE PEA APHID SPECIES COMPLEX

Linking adaptive divergence to hybrid unfitness is necessary to understand the ecological factors contributing to reproductive isolation and speciation. To date, this link has been demonstrated in few model systems, most of which encompass ecotypes that occupy relatively early stages in the speciation process. Here we extend these studies by assessing how host‐plant adaptation conditions hybrid fitness in the pea aphid, Acyrthosiphon pisum. We made crosses between and within five pea aphid biotypes adapted to different host plants and representing various stages of divergence within the complex. Performance of F1 hybrids and nonhybrids was assessed on a “universal” host that is favorable to all pea aphid biotypes in laboratory conditions. Although hybrids performed equally well as nonhybrids on the universal host, their performance was much lower than nonhybrids on the natural hosts of their parental populations. Hence, hybrids, rather than being intrinsically deficient, are maladapted to their parents’ hosts. Interestingly, the impact of this maladaptation was stronger in certain hybrids from crosses involving the most divergent biotype, suggesting that host‐dependent postzygotic isolation has continued to evolve late in divergence. Even though host‐independent deficiencies are not excluded, hybrid maladaptation to parental hosts supports the hypothesis of ecological speciation in this complex.