Genetic linkage of ecological specialization and reproductive isolation in pea aphids

Adaptive genetic divergence along narrow environmental gradients in four stream insects

A central question linking ecology with evolutionary biology is how environmental heterogeneity can drive adaptive genetic divergence among populations. We examined adaptive divergence of four stream insects from six adjacent catchments in Japan by combining field measures of habitat and resource components with genome scans of non-neutral Amplified Fragment Length Polymorphism (AFLP) loci. Neutral genetic variation was used to measure gene flow and non-neutral genetic variation was used to test for adaptive divergence. We identified the environmental characteristics contributing to divergence by comparing genetic distances at non-neutral loci between sites with Euclidean distances for each of 15 environmental variables. Comparisons were made using partial Mantel tests to control for geographic distance. In all four species, we found strong evidence for non-neutral divergence along environmental gradients at between 6 and 21 loci per species. The relative contribution of these environmental variables to each species' ecological niche was quantified as the specialization index, S, based on ecological data. In each species, the variable most significantly correlated with genetic distance at non-neutral loci was the same variable along which each species was most narrowly distributed (i.e., highest S). These were gradients of elevation (two species), chlorophyll-a, and ammonia-nitrogen. This adaptive divergence occurred in the face of ongoing gene flow (Fst = 0.01-0.04), indicating that selection was strong enough to overcome homogenization at the landscape scale. Our results suggest that adaptive divergence is pronounced, occurs along different environmental gradients for different species, and may consistently occur along the narrowest components of species' niche.

Divergent warning patterns contribute to assortative mating between incipient Heliconius species

Theoretical models suggest that traits under divergent ecological selection, which also contribute to assortative mating, will facilitate speciation with gene flow. Evidence for these so-called "magic traits" now exists across a range of taxa. However, their importance during speciation will depend on the extent to which they contribute to reproductive isolation. Addressing this requires experiments to determine the exact cues involved as well as estimates of assortative mating in the wild. Heliconius butterflies are well known for their diversity of bright warning color patterns, and their amenability to experimental manipulation has provided an excellent opportunity to test their role in reproductive isolation. Here, we reveal that divergent color patterns contribute to mate recognition between the incipient species Heliconius himera and H. erato, a taxon pair for which assortative mating by color pattern has been demonstrated among wild individuals: First, we demonstrate that males are more likely to attempt to mate conspecific females; second, we show that males are more likely to approach pinned females that share their own warning pattern. These data are valuable as these taxa likely represent the early stages of speciation, but unusually also allow comparisons with rates of interbreeding between divergent ecologically relevant phenotypes measured in the wild.

Differences in tolerance to host cactus alkaloids in Drosophila koepferae and D. buzzatii

The evolution of cactophily in the genus Drosophila was a major ecological transition involving over a hundred species in the Americas that acquired the capacity to cope with a variety of toxic metabolites evolved as feeding deterrents in Cactaceae. D. buzzatii and D. koepferae are sibling cactophilic species in the D. repleta group. The former is mainly associated with the relatively toxic-free habitat offered by prickly pears (Opuntia sulphurea) and the latter has evolved the ability to use columnar cacti of the genera Trichocereus and Cereus that contain an array of alkaloid secondary compounds. We assessed the effects of cactus alkaloids on fitness-related traits and evaluated the ability of D. buzzatii and D. koepferae to exploit an artificial novel toxic host. Larvae of both species were raised in laboratory culture media to which we added increasing doses of an alkaloid fraction extracted from the columnar cactus T. terschekii. In addition, we evaluated performance on an artificial novel host by rearing larvae in a seminatural medium that combined the nutritional quality of O. sulphurea plus amounts of alkaloids found in fresh T. terschekii. Performance scores in each rearing treatment were calculated using an index that took into account viability, developmental time, and adult body size. Only D. buzzatii suffered the effects of increasing doses of alkaloids and the artificial host impaired viability in D. koepferae, but did not affect performance in D. buzzatii. These results provide the first direct evidence that alkaloids are key determinants of host plant use in these species. However, the results regarding the artificial novel host suggest that the effects of alkaloids on performance are not straightforward as D. koepferae was heavily affected. We discuss these results in the light of patterns of host plan evolution in the Drosophila repleta group.

Pleiotropic effects of a single gene on skeletal development and sensory system patterning in sticklebacks

Background

Adaptation to a new environment can be facilitated by co-inheritance of a suite of phenotypes that are all advantageous in the new habitat. Although experimental evidence demonstrates that multiple phenotypes often map to the same genomic regions, it is challenging to determine whether phenotypes are associated due to pleiotropic effects of a single gene or to multiple tightly linked genes. In the threespine stickleback fish (Gasterosteus aculeatus), multiple phenotypes are associated with a genomic region around the gene Ectodysplasin (Eda), but only the presence of bony lateral plates has been conclusively shown to be caused by Eda.

Results

Here, we ask whether pleiotropy or linkage is responsible for the association between lateral plates and the distribution of the neuromasts of the lateral line. We first identify a strong correlation between plate appearance and changes in the spatial distribution of neuromasts through development. We then use an Eda transgene to induce the formation of ectopic plates in low plated fish, which also results in alterations to neuromast distribution. Our results also show that other loci may modify the effects of Eda on plate formation and neuromast distribution.

Conclusions

Together, these results demonstrate that Eda has pleiotropic effects on at least two phenotypes, highlighting the role of pleiotropy in the genetic basis of adaptation.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.

Genetic mapping of two components of reproductive isolation between two sibling species of moths, Ostrinia nubilalis and O. scapulalis

We report the quantitative trait loci (QTL) mapping of reproductive isolation traits between Ostrinia nubilalis (the European corn borer) and its sibling species O. scapulalis (the Adzuki bean borer), focusing on two traits: mating isolation (mi) and pheromone production (Pher). Four genetic maps were generated from two backcross families, with two maps (one chromosomal map and one linkage map) per backcross. We located 165-323 AFLP markers on these four maps, resulting in the identification of 27-31 linkage groups, depending on the map considered. No-choice mating experiments with the offspring of each backcross led to the detection of at least two QTLs for mi in different linkage groups. QTLs underlying Pher were located in a third linkage group. The Z heterochromosome was identified by a specific marker (Tpi) and did not carry any of these QTLs. Finally, we considered the global divergence between the two sibling species, distortions of segregation throughout the genome, and the location and effect of mi and Pher QTLs in light of the known candidate genes for reproductive isolation within the genus Ostrinia and, more broadly, in phytophagous insects.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

Is ecological speciation a major trend in aphids? Insights from a molecular phylogeny of the conifer-feeding genus Cinara

Introduction

In the past decade ecological speciation has been recognized as having an important role in the diversification of plant-feeding insects. Aphids are host-specialised phytophagous insects that mate on their host plants and, as such, they are prone to experience reproductive isolation linked with host plant association that could ultimately lead to species formation. The generality of such a scenario remains to be tested through macroevolutionary studies. To explore the prevalence of host-driven speciation in the diversification of the aphid genus Cinara and to investigate alternative modes of speciation, we reconstructed a phylogeny of this genus based on mitochondrial, nuclear and Buchnera aphidicola DNA sequence fragments and applied a DNA-based method of species delimitation. Using a recent software (PhyloType), we explored evolutionary transitions in host-plant genera, feeding sites and geographic distributions in the diversification of Cinara and investigated how transitions in these characters have accompanied speciation events.

Results

The diversification of Cinara has been constrained by host fidelity to conifer genera sometimes followed by sequential colonization onto different host species and by feeding-site specialisation. Nevertheless, our analyses suggest that, at the most, only half of the speciation events were accompanied by ecological niche shifts. The contribution of geographical isolation in the speciation process is clearly apparent in the occurrence of species from two continents in the same clades in relatively terminal positions in our phylogeny. Furthermore, in agreement with predictions from scenarios in which geographic isolation accounts for speciation events, geographic overlap between species increased significantly with time elapsed since their separation.

Conclusions

The history of Cinara offers a different perspective on the mode of speciation of aphids than that provided by classic models such as the pea aphid. In this genus of aphids, the role of climate and landscape history has probably been as important as host-plant specialisation in having shaped present-day diversity.

A locus in Drosophila sechellia affecting tolerance of a host plant toxin

Many insects feed on only one or a few types of host. These host specialists often evolve a preference for chemical cues emanating from their host and develop mechanisms for circumventing their host's defenses. Adaptations like these are central to evolutionary biology, yet our understanding of their genetics remains incomplete. Drosophila sechellia, an emerging model for the genetics of host specialization, is an island endemic that has adapted to chemical toxins present in the fruit of its host plant, Morinda citrifolia. Its sibling species, D. simulans, and many other Drosophila species do not tolerate these toxins and avoid the fruit. Earlier work found a region with a strong effect on tolerance to the major toxin, octanoic acid, on chromosome arm 3R. Using a novel assay, we narrowed this region to a small span near the centromere containing 18 genes, including three odorant binding proteins. It has been hypothesized that the evolution of host specialization is facilitated by genetic linkage between alleles contributing to host preference and alleles contributing to host usage, such as tolerance to secondary compounds. We tested this hypothesis by measuring the effect of this tolerance locus on host preference behavior. Our data were inconsistent with the linkage hypothesis, as flies bearing this tolerance region showed no increase in preference for media containing M. citrifolia toxins, which D. sechellia prefers. Thus, in contrast to some models for host preference, preference and tolerance are not tightly linked at this locus nor is increased tolerance per se sufficient to change preference. Our data are consistent with the previously proposed model that the evolution of D. sechellia as a M. citrifolia specialist occurred through a stepwise loss of aversion and gain of tolerance to M. citrifolia's toxins.

Evolution, selection and isolation: a genomic view of speciation in fungal plant pathogens

895 I. 895 II. 896 III. 898 IV. 900 V. 902 VI. 904 VII. 905 905 References 905 SUMMARY: Speciation of fungal plant pathogens has been associated with host jumps, host domestication, clonal divergence, and hybridization. Although we have substantial insight into the speciation histories of several important plant pathogens, we still know very little about the underlying genetics of reproductive isolation. Studies in Saccharomyces cerevisiae, Neurospora crassa, and nonfungal model systems illustrate that reproductive barriers can evolve by different mechanisms, including genetic incompatibilities between neutral and adaptive substitutions, reinforcement selection, and chromosomal rearrangements. Advances in genome sequencing and sequence analyses provide a new framework to identify those traits that have driven the divergence of populations or caused reproductive isolation between species of fungal plant pathogens. These traits can be recognized based on signatures of strong divergent selection between species or through the association of allelic combination conferring hybrid inferiority. Comparative genome analyses also provide information about the contribution of genome rearrangements to speciation. This is particularly relevant for species of fungal pathogens with extreme levels of genomic rearrangements and within-species genome plasticity.

Ecological and genetic factors influencing the transition between host-use strategies in sympatric Heliconius butterflies

Shifts in host-plant use by phytophagous insects have played a central role in their diversification. Evolving host-use strategies will reflect a trade-off between selection pressures. The ecological niche of herbivorous insects is partitioned along several dimensions, and if populations remain in contact, recombination will break down associations between relevant loci. As such, genetic architecture can profoundly affect the coordinated divergence of traits and subsequently the ability to exploit novel habitats. The closely related species Heliconius cydno and H. melpomene differ in mimetic colour pattern, habitat and host-plant use. We investigate the selection pressures and genetic basis underlying host-use differences in these two species. Host-plant surveys reveal that H. melpomene specializes on a single species of Passiflora. This is also true for the majority of other Heliconius species in secondary growth forest at our study site, as expected under a model of interspecific competition. In contrast, H. cydno, which uses closed-forest habitats where both Heliconius and Passiflora are less common, appears not to be restricted by competition and uses a broad selection of the available Passiflora. However, other selection pressures are likely involved, and field experiments reveal that early larval survival of both butterfly species is highest on Passiflora menispermifolia, but most markedly so for H. melpomene, the specialist on that host. Finally, we demonstrate an association between host-plant acceptance and colour pattern amongst interspecific hybrids, suggesting that major loci underlying these important ecological traits are physically linked in the genome. Together, our results reveal ecological and genetic associations between shifts in habitat, host use and mimetic colour pattern that have likely facilitated both speciation and coexistence.

Inferring the population structure of Myzus persicae in diverse agroecosystems using microsatellite markers

Diverse agroecosystems offer phytophagous insects a wide choice of host plants. Myzus persicae is a polyphagous aphid common in moderate climates. During its life cycle it alternates between primary and secondary hosts. A spatial genetic population structure may arise due to environmental factors and reproduction modes. The aim of this work was to determine the spatial and temporal genetic population structure of M. persicae in relation to host plants and climatic conditions. For this, 923 individuals of M. persicae collected from six plant families between 2005 and 2008 in south-eastern Spain were genotyped for eight microsatellite loci. The population structure was inferred by neighbour-joining, analysis of molecular variance (AMOVA) and Bayesian analyses. Moderate polymorphism was observed for the eight loci in almost all the samples. No differences in the number of alleles were observed between primary and secondary hosts or between geographical areas. The proportion of unique genotypes found in the primary host was similar in the north (0.961 ± 0.036) and the south (0.987 ± 0.013), while in the secondary host it was higher in the north (0.801 ± 0.159) than in the south (0.318 ± 0.063). Heterozygosity excess and linkage disequilibrium suggest a high representation of obligate parthenogens in areas with warmer climate and in the secondary hosts. The F ST-values pointed to no genetic differentiation of M. persicae on the different plant families. F ST-values, AMOVA and Bayesian model-based cluster analyses pointed to a significant population structure that was related to primary and secondary hosts. Differences between primary and secondary hosts could be due to the overrepresentation of parthenogens on herbaceous plants.

Evolution of mate choice and the so-called magic traits in ecological speciation

Non-random mating provides multiple evolutionary benefits and can result in speciation. Biological organisms are characterised by a myriad of different traits, many of which can serve as mating cues. We consider multiple mechanisms of non-random mating simultaneously within a unified modelling framework in an attempt to understand better which are more likely to evolve in natural populations going through the process of local adaptation and ecological speciation. We show that certain traits that are under direct natural selection are more likely to be co-opted as mating cues, leading to the appearance of magic traits (i.e. phenotypic traits involved in both local adaptation and mating decisions). Multiple mechanisms of non-random mating can interact so that trait co-evolution enables the evolution of non-random mating mechanisms that would not evolve alone. The presence of magic traits may suggest that ecological selection was acting during the origin of new species.

Environmental dependence of genetic constraint

The epistatic interactions that underlie evolutionary constraint have mainly been studied for constant external conditions. However, environmental changes may modulate epistasis and hence affect genetic constraints. Here we investigate genetic constraints in the adaptive evolution of a novel regulatory function in variable environments, using the lac repressor, LacI, as a model system. We have systematically reconstructed mutational trajectories from wild type LacI to three different variants that each exhibit an inverse response to the inducing ligand IPTG, and analyzed the higher-order interactions between genetic and environmental changes. We find epistasis to depend strongly on the environment. As a result, mutational steps essential to inversion but inaccessible by positive selection in one environment, become accessible in another. We present a graphical method to analyze the observed complex higher-order interactions between multiple mutations and environmental change, and show how the interactions can be explained by a combination of mutational effects on allostery and thermodynamic stability. This dependency of genetic constraint on the environment should fundamentally affect evolutionary dynamics and affects the interpretation of phylogenetic data.

Neo-sex chromosomes and adaptive potential in tortricid pests

Changes in genome architecture often have a significant effect on ecological specialization and speciation. This effect may be further enhanced by involvement of sex chromosomes playing a disproportionate role in reproductive isolation. We have physically mapped the Z chromosome of the major pome fruit pest, the codling moth, Cydia pomonella (Tortricidae), and show that it arose by fusion between an ancestral Z chromosome and an autosome corresponding to chromosome 15 in the Bombyx mori reference genome. We further show that the fusion originated in a common ancestor of the main tortricid subfamilies, Olethreutinae and Tortricinae, comprising almost 700 pest species worldwide. The Z-autosome fusion brought two major genes conferring insecticide resistance and clusters of genes involved in detoxification of plant secondary metabolites under sex-linked inheritance. We suggest that this fusion significantly increased the adaptive potential of tortricid moths and thus contributed to their radiation and subsequent speciation.

Cospeciation vs host-shift speciation: methods for testing, evidence from natural associations and relation to coevolution

Hosts and their symbionts are involved in intimate physiological and ecological interactions. The impact of these interactions on the evolution of each partner depends on the time-scale considered. Short-term dynamics - 'coevolution' in the narrow sense - has been reviewed elsewhere. We focus here on the long-term evolutionary dynamics of cospeciation and speciation following host shifts. Whether hosts and their symbionts speciate in parallel, by cospeciation, or through host shifts, is a key issue in host-symbiont evolution. In this review, we first outline approaches to compare divergence between pairwise associated groups of species, their advantages and pitfalls. We then consider recent insights into the long-term evolution of host-parasite and host-mutualist associations by critically reviewing the literature. We show that convincing cases of cospeciation are rare (7%) and that cophylogenetic methods overestimate the occurrence of such events. Finally, we examine the relationships between short-term coevolutionary dynamics and long-term patterns of diversification in host-symbiont associations. We review theoretical and experimental studies showing that short-term dynamics can foster parasite specialization, but that these events can occur following host shifts and do not necessarily involve cospeciation. Overall, there is now substantial evidence to suggest that coevolutionary dynamics of hosts and parasites do not favor long-term cospeciation.

Genetic Diversity and Population Structure of Busseola segeta Bowden (Lepidoptera; Noctuidae): A Case Study of Host Use Diversification in Guineo-Congolian Rainforest Relic Area, Kenya

Habitat modification and fragmentation are considered as some of the factors that drive organism distribution and host use diversification. Indigenous African stem borer pests are thought to have diversified their host ranges to include maize [Zea mays L.] and sorghum [Sorghum bicolor (L.) Moench] in response to their increased availability through extensive cultivation. However, management efforts have been geared towards reducing pest populations in the cultivated fields with few attempts to understand possible evolution of "new" pest species. Recovery and growing persistence of Busseola segeta Bowden on maize (Zea mays L.) in Kakamega called for studies on the role of wild host plants on the invasion of crops by wild borer species. A two-year survey was carried out in a small agricultural landscape along the edge of Kakamega forest (Kenya) to assess host range and population genetic structure of B. segeta. The larvae of B. segeta were found on nine different plant species with the majority occurring on maize and sorghum. Of forty cytochrome b haplotypes identified, twenty-three occurred in both wild and cultivated habitats. The moths appear to fly long distances across the habitats with genetic analyses revealing weak differentiation between hosts in different habitats (F_ST = 0.016; p = 0.015). However, there was strong evidence of variation in genetic composition between growing seasons in the wild habitat (F_ST = 0.060; p < 0.001) with emergence or disappearance of haplotypes between habitats. Busseola segeta is an example of a phytophagous insect that utilizes plants with a human induced distribution range, maize, but does not show evidence of host race formation or reduction of gene flow among populations using different hosts. However, B. segeta is capable of becoming an important pest in the area and the current low densities may be attributed to the general low infestation levels and presence of a wide range of alternative hosts in the area.

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome-wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host-adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.

Quantitative genetic analysis of subspecific differences in body shape in the snail-feeding carabid beetle Damaster blaptoides

A dimorphic pattern of macrocephalic (wide, short) and stenocephalic (narrow, long) body shapes is observed in snail-feeding carabid beetles globally. The former exhibits high performance in crushing snail shells with powerful jaws, whereas the latter specializes in eating snails' soft body directly by inserting the head into the shell. In the snail-feeding species Damaster blaptoides, the subspecies D. b. capito has a wide, short forebody, and D. b. fortunei has a narrow, long forebody. They exhibit distinct morphologies despite their geographic and phylogenetic proximity. To examine the genetic basis of the morphological differences between these two subspecies, we conducted quantitative genetic analyses by crossing these subspecies and producing F(1) and backcross hybrids. The hybrids had body shapes intermediate between the parental subspecies. The variation between wide, short and narrow, long forebodies was based on negative genetic correlations between width and length of the head and thorax. Between one and eight genetic factors were involved in the morphological differences between subspecies. We suggest that the morphological integration of forebody parts in a small number of loci has facilitated the marked morphological diversification between subspecies of D. blaptoides.

Localizing F(ST) outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation-with-gene-flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence-with-gene-flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F(ST) outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between-population recombination and gene exchange are reduced over large 'divergence hitchhiking' (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F(ST) outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation-with-gene-flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.

Specificity, rank preference, and the colonization of a non-native host plant by the Melissa blue butterfly

Animals often express behavioral preferences for different types of food or other resources, and these preferences can evolve or shift following association with novel food types. Shifts in preference can involve at least two phenomena: a change in rank preference or a change in specificity. The former corresponds to a change in the order in which hosts are preferred, while a shift in specificity can be an increase in the tendency to utilize multiple hosts. These possibilities have been examined in relatively few systems that include extensive population-level replication. The Melissa blue butterfly, Lycaeides melissa, has colonized exotic alfalfa, Medicago sativa, throughout western North America. We assayed the host preferences of 229 females from ten populations associated with novel and native hosts. In four out of five native-associated populations, a native host was preferred over the exotic host, while preference for a native host characterized only two out of five of the alfalfa-associated populations. Across all individuals from alfalfa-associated populations, there appears to have been a decrease in specificity: females from these populations lay fewer eggs on the native host and more eggs on the exotic relative to females from native-host populations. However, females from alfalfa-associated populations did not lay more eggs on a third plant species, which suggests that preferences for specific hosts in this system can potentially be gained and lost independently. Geographic variation in oviposition preference in L. melissa highlights the value of surveying a large number of populations when studying the evolution of a complex behavioral trait.

The genetic architecture of a complex ecological trait: host plant use in the specialist moth, Heliothis subflexa

We used genetic mapping to examine the genetic architecture of differences in host plant use between two species of noctuid moths, Heliothis subflexa, a specialist on Physalis spp., and its close relative, the broad generalist H. virescens. We introgressed H. subflexa chromosomes into the H. virescens background and analyzed 1462 backcross insects. The effects of H. subflexa-origin chromosomes were small when measured as the percent variation explained in backcross populations (0.2-5%), but were larger when considered in relation to the interspecific difference explained (1.5-165%). Most significant chromosomes had effects on more than one trait, and their effects varied between years, sexes, and genetic backgrounds. Different chromosomes could produce similar phenotypes, suggesting that the same trait might be controlled by different chromosomes in different backcross populations. It appears that many loci of small effect contribute to the use of Physalis by H. subflexa. We hypothesize that behavioral changes may have paved the way for physiological adaptation to Physalis by the generalist ancestor of H. subflexa and H. virescens.

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris

BACKGROUND

Hybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years.

RESULTS

QTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups.

CONCLUSIONS

Two different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).

Admixture mapping of male nuptial colour and body shape in a recently formed hybrid population of threespine stickleback

Despite recent progress, we still know relatively little about the genetic architecture that underlies adaptation to divergent environments. Determining whether the genetic architecture of phenotypic adaptation follows any predictable patterns requires data from a wide variety of species. However, in many organisms, genetic studies are hindered by the inability to perform genetic crosses in the laboratory or by long generation times. Admixture mapping is an approach that circumvents these issues by taking advantage of hybridization that occurs between populations or species in the wild. Here, we demonstrate the utility of admixture mapping in a naturally occurring hybrid population of threespine sticklebacks (Gasterosteus aculeatus) from Enos Lake, British Columbia. Until recently, this lake contained two species of sticklebacks adapted to divergent habitats within the lake. This benthic-limnetic species pair diverged in a number of phenotypes, including male nuptial coloration and body shape, which were previously shown to contribute to reproductive isolation between them. However, recent ecological disturbance has contributed to extensive hybridization between the species, and there is now a single, admixed population within Enos Lake. We collected over 500 males from Enos Lake and found that most had intermediate nuptial colour and body shape. By genotyping males with nuptial colour at the two extremes of the phenotypic distribution, we identified seven genomic regions on three chromosomes associated with divergence in male nuptial colour. These genomic regions are also associated with variation in body shape, suggesting that tight linkage and/or pleiotropy facilitated adaptation to divergent environments in benthic-limnetic species pairs.

The genetic basis of a complex functional system

The relationship between form and function can have profound effects on evolutionary dynamics and such effects may differ for simple versus complex systems. In particular, functions produced by multiple structural configurations (many-to-one mapping, MTOM) may dampen constituent trade-offs and promote diversification. Unfortunately, we lack information about the genetic architecture of MTOM functional systems. The skulls of teleost fishes contain both simple (lower jaw levers) as well as more complex (jaws modeled as 4-bar linkages) functional systems within the same craniofacial unit. We examined the mapping of form to function and the genetic basis of these systems by identifying quantitative trait loci (QTL) in hybrids of two Lake Malawi cichlid species. Hybrid individuals exhibited novelty (transgressive segregation) in morphological components and function of the simple and complex jaw systems. Functional novelty was proportional to the prevalence of extreme morphologies in the simple levers; by contrast, recombination of parental morphologies produced transgression in the MTOM 4-bar linkage. We found multiple loci of moderate effect and epistasis controlling jaw phenotypes in both the simple and complex systems, with less phenotypic variance explained by QTL for the 4-bar. Genetic linkage between components of the simple and complex systems partly explains phenotypic correlations and may constrain functional evolution.

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.

Island phytophagy: explaining the remarkable diversity of plant-feeding insects

Plant-feeding insects have undergone unparalleled diversification among different plant taxa, yet explanations for variation in their diversity lack a quantitative, predictive framework. Island biogeographic theory has been applied to spatially discrete habitats but not to habitats, such as host plants, separated by genetic distance. We show that relationships between the diversity of gall-inducing flies and their host plants meet several fundamental predictions from island biogeographic theory. First, plant-taxon genetic distinctiveness, an integrator for long-term evolutionary history of plant lineages, is a significant predictor of variance in the diversity of gall-inducing flies among host-plant taxa. Second, range size and structural complexity also explain significant proportions of the variance in diversity of gall-inducing flies among different host-plant taxa. Third, as with other island systems, plant-lineage age does not predict species diversity. Island biogeographic theory, applied to habitats defined by genetic distance, provides a novel, comprehensive framework for analysing and explaining the diversity of plant-feeding insects and other host-specific taxa.

Unifying concepts and mechanisms in the specificity of plant-enemy interactions

Host ranges are commonly quantified to classify herbivores and plant pathogens as either generalists or specialists. Here, we summarize patterns and mechanisms in the interactions of plants with these enemies along different axes of specificity. We highlight the many dimensions within which plant enemies can specify and consider the underlying ecological, evolutionary and molecular mechanisms. Host resistance traits and enemy effectors emerge as central players determining host utilization and thus host range. Finally, we review approaches to studying the causes and consequences of variation in the specificity of plant-enemy interactions. Knowledge of the molecular mechanisms that determine host range is required to understand host shifts, and evolutionary transitions among specialist and generalist strategies, and to predict potential host ranges of pathogens and herbivores.

Speciation by symbiosis

In the Origin of Species, Darwin struggled with how continuous changes within a species lead to the emergence of discrete species. Molecular analyses have since identified nuclear genes and organelles that underpin speciation. In this review, we explore the microbiota as a third genetic component that spurs species formation. We first recall Ivan Wallin's original conception from the early 20th century on the role that bacteria play in speciation. We then describe three fundamental observations that justify a prominent role for microbes in eukaryotic speciation, consolidate exemplar studies of microbe-assisted speciation and incorporate the microbiota into classic models of speciation.

The genomics of speciation-with-gene-flow

The emerging field of speciation genomics is advancing our understanding of the evolution of reproductive isolation from the individual gene to a whole-genome perspective. In this new view it is important to understand the conditions under which 'divergence hitchhiking' associated with the physical linkage of gene regions, versus 'genome hitchhiking' associated with reductions in genome-wide rates of gene flow caused by selection, can enhance speciation-with-gene-flow. We describe here a theory predicting four phases of speciation, defined by changes in the relative effectiveness of divergence and genome hitchhiking, and review empirical data in light of the theory. We outline future directions, emphasizing the need to couple next-generation sequencing with selection, transplant, functional genomics, and mapping studies. This will permit a natural history of speciation genomics that will help to elucidate the factors responsible for population divergence and the roles that genome structure and different forms of hitchhiking play in facilitating the genesis of new biodiversity.

Large-scale candidate gene scan reveals the role of chemoreceptor genes in host plant specialization and speciation in the pea aphid

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

Genetic trade-offs and conditional neutrality contribute to local adaptation

Divergent natural selection promotes local adaptation and can lead to reproductive isolation of populations in contrasting environments; however, the genetic basis of local adaptation remains largely unresolved in natural populations. Local adaptation might result from antagonistic pleiotropy, where alternate alleles are favoured in distinct habitats, and polymorphism is maintained by selection. Alternatively, under conditional neutrality some alleles may be favoured in one environment but neutral at other locations. Antagonistic pleiotropy maintains genetic variation across the landscape; however, there is a systematic bias against discovery of antagonistic pleiotropy because the fitness benefits of local alleles need to be significant in at least two environments. Here, we develop a generally applicable method to investigate polygenic local adaptation and identify loci that are the targets of selection. This approach evaluates allele frequency changes after selection at loci across the genome to distinguish antagonistic pleiotropy from conditional neutrality and deleterious variation. We investigate local adaptation at the qualitative trait loci (QTL) level in field experiments, in which we expose 177 F(6) recombinant inbred lines and parental lines of Boechera stricta (Brassicaceae) to their parental environments over two seasons. We demonstrate polygenic selection for native alleles in both environments, with 2.8% of the genome exhibiting antagonistic pleiotropy and 8% displaying conditional neutrality. Our study strongly supports antagonistic pleiotropy at one large-effect flowering phenology QTL (nFT): native homozygotes had significantly greater probabilities of flowering than foreign homozygotes in both parental environments. Such large-scale field studies are essential to elucidate the genetic basis of adaptation in natural populations.

Habitat Choice and Speciation

The role of habitat choice in reproductive isolation and ecological speciation has often been overlooked, despite acknowledgement of its ability to facilitate local adaptation. It can form part of the speciation process through various evolutionary mechanisms, yet where habitat choice has been included in models of ecological speciation little thought has been given to these underlying mechanisms. Here, we propose and describe three independent criteria underlying ten different evolutionary scenarios in which habitat choice may promote or maintain local adaptation. The scenarios are the result of all possible combinations of the independent criteria, providing a conceptual framework in which to discuss examples which illustrate each scenario. These examples show that the different roles of habitat choice in ecological speciation have rarely been effectively distinguished. Making such distinctions is an important challenge for the future, allowing better experimental design, stronger inferences and more meaningful comparisons among systems. We show some of the practical difficulties involved by reviewing the current evidence for the role of habitat choice in local adaptation and reproductive isolation in the intertidal gastropodLittorina saxatilis, a model system for the study of ecological speciation, assessing whether any of the proposed scenarios can be reliably distinguished, given current research.

Fatty-acid preference changes during development in Drosophila melanogaster

Fatty-acids (FAs) are required in the diet of many animals throughout their life. However, the mechanisms involved in the perception of and preferences for dietary saturated and unsaturated FAs (SFAs and UFAs, respectively) remain poorly explored, especially in insects. Using the model species Drosophila melanogaster, we measured the responses of wild-type larvae and adults to pure SFAs (14, 16, and 18 carbons) and UFAs (C18 with 1, 2, or 3 double-bonds). Individual and group behavioral tests revealed different preferences in larvae and adults. Larvae preferred UFAs whereas SFAs tended to induce both a strong aversion and a persistent aggregation behavior. Adults generally preferred SFAs, and laid more eggs and had a longer life span when ingesting these substances as compared to UFAs. Our data suggest that insects can discriminate long-chain dietary FAs. The developmental change in preference shown by this species might reflect functional variation in use of FAs or stage-specific nutritional requirements, and may be fundamental for insect use of these major dietary components.

Hybrid breakdown weakens under thermal stress in population crosses of the copepod Tigriopus californicus

The outcome of hybridization can be impacted by environmental conditions, which themselves can contribute to reproductive isolation between taxa. In crosses of genetically divergent populations, hybridization can have both negative and positive impacts on fitness, the balance between which might be tipped by changes in the environment. Genetically divergent populations of the intertidal copepod Tigriopus californicus have been shown to differ in thermal tolerance at high temperatures along a latitudinal gradient. In this study, a series of crosses were made between pairs of genetically divergent populations of T. californicus, and the thermal tolerance of these hybrids was tested. In most cases, the first-generation hybrids had relatively high thermal tolerance and the second-generation hybrids were not generally reduced below the less-tolerant parental population for high temperature tolerance. This pattern contrasts with previous studies from crosses of genetically divergent populations of this copepod, which often shows hybrid breakdown in these second-generation hybrids for other measures of fitness. These results suggest that high temperature stress could either increase the positive impacts of hybridization or decrease the negative impacts of hybridization resulting in lowered hybrid breakdown in these population crosses.