Genetic linkage of ecological specialization and reproductive isolation in pea aphids

Ecological divergence despite common mating sites: Genotypes and symbiotypes shed light on cryptic diversity in the black bean aphid species complex

Different host plants represent ecologically dissimilar environments for phytophagous insects. The resulting divergent selection can promote the evolution of specialized host races, provided that gene flow is reduced between populations feeding on different plants. In black bean aphids belonging to the Aphis fabae complex, several morphologically cryptic taxa have been described based on their distinct host plant preferences. However, host choice and mate choice are largely decoupled in these insects: they are host-alternating and migrate between specific summer host plants and shared winter hosts, with mating occurring on the shared hosts. This provides a yearly opportunity for gene flow among aphids using different summer hosts, and raises the question if and to what extent the ecologically defined taxa are reproductively isolated. Here, we analyzed a geographically and temporally structured dataset of microsatellite genotypes from A. fabae that were mostly collected from their main winter host Euonymus europaeus, and additionally from another winter host and fourteen summer hosts. The data reveals multiple, strongly differentiated genetic clusters, which differ in their association with different summer and winter hosts. The clusters also differ in the frequency of infection with two heritable, facultative endosymbionts, separately hinting at reproductive isolation and divergent ecological selection. Furthermore, we found evidence for occasional hybridization among genetic clusters, with putative hybrids collected more frequently in spring than in autumn. This suggests that similar to host races in other phytophagous insects, both prezygotic and postzygotic barriers including selection against hybrids maintain genetic differentiation among A. fabae taxa, despite a common mating habitat.

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

The effectiveness of pseudomagic traits in promoting premating isolation

Upon the secondary contact of populations, speciation with gene flow is greatly facilitated when the same pleiotropic loci are both subject to divergent ecological selection and induce non-random mating, leading to loci with this fortuitous combination of functions being referred to as 'magic trait' loci. We use a population genetics model to examine whether 'pseudomagic trait' complexes, composed of physically linked loci fulfilling these two functions, are as efficient in promoting premating isolation as magic traits. We specifically measure the evolution of choosiness, which controls the strength of assortative mating. We show that, surprisingly, pseudomagic trait complexes, and to a lesser extent also physically unlinked loci, can lead to the evolution of considerably stronger assortative mating preferences than do magic traits, provided polymorphism at the involved loci is maintained. This is because assortative mating preferences are generally favoured when there is a risk of producing maladapted recombinants, as occurs with non-magic trait complexes but not with magic traits (since pleiotropy precludes recombination). Contrary to current belief, magic traits may not be the most effective genetic architecture for promoting strong premating isolation. Therefore, distinguishing between magic traits and pseudomagic trait complexes is important when inferring their role in premating isolation. This calls for further fine-scale genomic research on speciation genes.

Molecular Mechanisms Underlying Host Plant Specificity in Aphids

Aphids are serious pests of agricultural and ornamental plants and important model systems for hemipteran-plant interactions. The long evolutionary history of aphids with their host plants has resulted in a variety of systems that provide insight into the different adaptation strategies of aphids to plants and vice versa. In the past, various plant-aphid interactions have been documented, but lack of functional tools has limited molecular studies on the mechanisms of plant-aphid interactions. Recent technological advances have begun to reveal plant-aphid interactions at the molecular level and to increase our knowledge of the mechanisms of aphid adaptation or specialization to different host plants. In this article, we compile and analyze available information on plant-aphid interactions, discuss the limitations of current knowledge, and argue for new research directions. We advocate for more work that takes advantage of natural systems and recently established molecular techniques to obtain a comprehensive view of plant-aphid interaction mechanisms.

Evolution of physical linkage between loci controlling ecological traits and mating preferences

Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.

Population Genetics of the Blueberry Gall Midge, Dasineura oxycoccana (Diptera: Cecidomyiidae), on Blueberry and Cranberry and Testing Invasion Scenarios

We compared the population genetic structure between populations of the blueberry gall midge-Dasineura oxycoccana (Johnson) (Diptera: Cecidomyiidae)-from blueberry and cranberry and determined the genetic relationships among geographical subgroups by genotyping 632 individuals from 31 different populations from their native USA regions (New Jersey, Michigan, and Georgia) and from invaded Korean regions using 12 microsatellite loci. Our population genetic analyses showed a clear separation between the two host-associated D. oxycoccana populations from blueberry and cranberry. Using data from only the blueberry-associated D. oxycoccana populations, we identified five genetically isolated subgroups. An analysis of the approximate Bayesian computation suggests that the invasive D. oxycoccana population from Korea appears to have been introduced from an unsampled source population rather than directly from its native range. Our findings will allow for an easier identification of the source of D. oxycoccana into newly invaded regions, as well as to determine their association with blueberry and cranberry, which based on our results can be considered as two distinct species.

Additive genetic effects in interacting species jointly determine the outcome of caterpillar herbivory

Plant-insect interactions are common and important in basic and applied biology. Trait and genetic variation can affect the outcome and evolution of these interactions, but the relative contributions of plant and insect genetic variation and how these interact remain unclear and are rarely subject to assessment in the same experimental context. Here, we address this knowledge gap using a recent host-range expansion onto alfalfa by the Melissa blue butterfly. Common garden rearing experiments and genomic data show that caterpillar performance depends on plant and insect genetic variation, with insect genetics contributing to performance earlier in development and plant genetics later. Our models of performance based on caterpillar genetics retained predictive power when applied to a second common garden. Much of the plant genetic effect could be explained by heritable variation in plant phytochemicals, especially saponins, peptides, and phosphatidyl cholines, providing a possible mechanistic understanding of variation in the species interaction. We find evidence of polygenic, mostly additive effects within and between species, with consistent effects of plant genotype on growth and development across multiple butterfly species. Our results inform theories of plant-insect coevolution and the evolution of diet breadth in herbivorous insects and other host-specific parasites.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Evolution of parasitoid host preference and performance in response to an invasive host acting as evolutionary trap

The invasion of a novel host species can create a mismatch in host choice and offspring survival (performance) when native parasitoids attempt to exploit the invasive host without being able to circumvent its resistance mechanisms. Invasive hosts can therefore act as evolutionary trap reducing parasitoids' fitness and this may eventually lead to their extinction. Yet, escape from the trap can occur when parasitoids evolve behavioral avoidance or a physiological strategy compatible with the trap host, resulting in either host‐range expansion or a complete host‐shift. We developed an individual based model to investigate which conditions promote parasitoids to evolve behavioral preference that matches their performance, including host‐trap avoidance, and which conditions lead to adaptations to the unsuitable hosts. The model was inspired by solitary endo‐parasitoids attacking larval host stages. One important aspect of these conditions was reduced host survival during incompatible interaction, where a failed parasitization attempt by a parasitoid resulted not only in death of her offspring but also in host killing. This non‐reproductive host mortality had a strong influence on the likelihood of establishment of novel host–parasitoid relationship, in some cases constraining adaptation to the trap host species. Moreover, our model revealed that host‐search efficiency and genetic variation in host‐preference play a key role in the likelihood that parasitoids will include the suboptimal host in their host range, or will evolve behavioral avoidance resulting in specialization and host‐range conservation, respectively. Hence, invasive species might change the evolutionary trajectory of native parasitoid species, which is important for predicting biocontrol ability of native parasitoids towards novel hosts.

Interactions With Plant Defences Isolate Sympatric Populations of an Herbivorous Mite

Host plant specialisation can promote evolutionary divergence between herbivore populations associated with different plant species. While the mechanisms by which specialist species exploit their hosts have been studied widely across taxa, less is known about the mechanisms that allow intraspecific variants to arise and to be maintained across spatial and temporal scales. To understand whether adaptations to plant defences against herbivory contribute to the co-existence of genetically distinct populations of an herbivore, we investigate the interaction between honeysuckle (Lonicera periclymenum) and sympatric specialist and generalist populations of the spider mite Tetranychus urticae. We found that mite folivory induces the production of sticky droplets on honeysuckle, which have a defensive role: they increase mite mortality directly, and potentially indirectly by increasing the arrestment of a predator. We show that droplet induction and the preference to feed on honeysuckle depend on mite genotype, where the generalist avoids this host and the specialist suppresses droplet production. These traits are heritable and dominant in F1 hybrids between generalists and specialists. Selection pressure from honeysuckle and differences in host preference likely reduce the opportunity of mating encounters on this host. We propose that the interplay between selection from host plant defences and ecological barriers to hybridisation contribute to the persistence of genetically distinct populations of a single species in sympatry.

Improving Natural Enemy Selection in Biological Control through Greater Attention to Chemical Ecology and Host-Associated Differentiation of Target Arthropod Pests

Host-associated differentiation (HAD) refers to cases in which genetically distinct populations of a species (e.g., herbivores or natural enemies) preferentially reproduce or feed on different host species. In agroecosystems, HAD often results in unique strains or biotypes of pest species, each attacking different species of crops. However, HAD is not restricted to pest populations, and may cascade to the third trophic level, affecting host selection by natural enemies, and ultimately leading to HAD within natural enemy species. Natural enemy HAD may affect the outcomes of biological control efforts, whether classical, conservation, or augmentative. Here, we explore the potential effects of pest and natural enemy HAD on biological control in agroecosystems, with emphases on current knowledge gaps and implications of HAD for selection of biological control agents. Additionally, given the importance of semiochemicals in mediating interactions between trophic levels, we emphasize the role of chemical ecology in interactions between pests and natural enemies, and suggest areas of consideration for biological control. Overall, we aim to jump-start a conversation concerning the relevance of HAD in biological control by reviewing currently available information on natural enemy HAD, identifying challenges to incorporating HAD considerations into biological control efforts, and proposing future research directions on natural enemy selection and HAD.

The Build-Up of Population Genetic Divergence along the Speciation Continuum during a Recent Adaptive Radiation of Rhagoletis Flies

New species form through the evolution of genetic barriers to gene flow between previously interbreeding populations. The understanding of how speciation proceeds is hampered by our inability to follow cases of incipient speciation through time. Comparative approaches examining different diverging taxa may offer limited inferences, unless they fulfill criteria that make the comparisons relevant. Here, we test for those criteria in a recent adaptive radiation of the Rhagoletis pomonella species group (RPSG) hypothesized to have diverged in sympatry via adaptation to different host fruits. We use a large-scale population genetic survey of 1568 flies across 33 populations to: (1) detect on-going hybridization, (2) determine whether the RPSG is derived from the same proximate ancestor, and (3) examine patterns of clustering and differentiation among sympatric populations. We find that divergence of each in-group RPSG taxon is occurring under current gene flow, that the derived members are nested within the large pool of genetic variation present in hawthorn-infesting populations of R. pomonella, and that sympatric population pairs differ markedly in their degree of genotypic clustering and differentiation across loci. We conclude that the RPSG provides a particularly robust opportunity to make direct comparisons to test hypotheses about how ecological speciation proceeds despite on-going gene flow.

Population genetic structure of the carrot weevil ( Listronotus oregonensis ) in North America

Population genetic studies of insect pests enhance our ability to anticipate problems in agroecosystems, such as pest outbreaks, insecticide resistance, or expansions of the host range. This study focuses on geographic distance and host plant selection as potential determinants of genetic differentiation of the carrot weevil Listronotus oregonensis, a major pest of several apiaceous crops in North America. To undertake genetic studies on this species, we assembled the first complete genome sequence for L. oregonensis. Then, we used both haplotype discrimination with mitochondrial DNA (mtDNA) and a genotyping‐by‐sequencing (GBS) approach to characterize the genetic population structure. A total of 220 individuals were sampled from 17 localities in the provinces of Québec, Ontario, Nova Scotia (Canada), and the state of Ohio (USA). Our results showed significant genetic differences between distant populations across North America, indicating that geographic distance represents an important factor of differentiation for the carrot weevil. Furthermore, the GBS analysis revealed more different clusters than COI analysis between Québec and Nova Scotia populations, suggesting a recent differentiation in the latter province. In contrast, we found no clear evidence of population structure associated with the four cultivated apiaceous plants tested (carrot, parsley, celery, and celeriac) using populations from Québec. This first characterization of the genetic structure of the carrot weevil contributes to a better understanding of the gene flow of the species and helps to adapt local pest management measures to better control this agricultural pest.

Plant pathogen-mediated rapid acclimation of a host-specialized aphid to a non-host plant

Polyphagous aphids often consist of host-specialized lineages, which have greater fitness on their native hosts than on others. The underlying causes are important for understanding of the evolution of diet breadth and host shift of aphids. The cotton-melon aphid Aphis gossypii Glover is extremely polyphagous with many strict host-specialized lineages. Whether and how the lineage specialized on the primary host hibiscus shifts to the secondary host cucumber remains elusive. We found that the hibiscus-specialized lineage suffered high mortality and gave birth to very few nymphs developing into yellow dwarfs on fresh cucumber leaves, and did not inflict any damage symptoms on cucumber plants. The poor performance did not improve with prolonged exposure to cucumber; however, it did significantly improve when the cucumber leaves were pre-infected with a biotrophic phytopathogen Pseudoperonospora cubensis. More importantly, the hibiscus-specialized lineage with two-generation feeding experience on pre-infected cucumber leaves performed as well as the cucumber-specialized lineage did on fresh cucumber leaves, and inflicted typical damage symptoms on intact cucumber plants. Electrical penetration graph (EPG) indicated that the hibiscus-specialized lineage did not ingest phloem sap from fresh cucumber leaves but succeeded in ingesting phloem sap from pre-infected cucumber leaves, which explained the performance improvement of the hibiscus-specialized lineage on pre-infected cucumber leaves. This study revealed a new pathway for the hibiscus-specialized lineage to quickly acclimate to cucumber under the assistance of the phytopathogen. We considered that the short feeding experience on pre-infected cucumber may activate expression of effector genes that are related to specific host utilization. We suggest to identify host-specific effectors by comparing proteomes or/and transcriptomes of the hibiscus-specialized lineage before and after acclimating to cucumber.

Light environment influences mating behaviours during the early stages of divergence in tropical butterflies

Speciation is facilitated when traits under divergent selection also act as mating cues. Fluctuations in sensory conditions can alter signal perception independently of adaptation to the broader sensory environment, but how this fine-scale variation may constrain or promote behavioural isolation has received little attention. The warning patterns of Heliconius butterflies are under selection for aposematism and act as mating cues. Using computer vision, we extracted behavioural data from 1481 h of video footage, for 387 individuals. We show that the putative hybrid species H. heurippa and its close relative H. timareta linaresi differ in their response to divergent warning patterns, but that these differences are strengthened with increased local illuminance. Trials with live individuals reveal low-level assortative mating that is sufficiently explained by differences in visual attraction. Finally, results from hybrid butterflies are consistent with linkage between a major warning pattern gene and the corresponding behaviour, though the differences in behaviour we observe are unlikely to cause rapid reproductive isolation as predicted under a model of hybrid trait speciation. Overall, our results reveal that the contribution of ecological mating cues to reproductive isolation may depend on the immediate sensory conditions during which they are displayed to conspecifics.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

What We Don't Know About Diet-Breadth Evolution in Herbivorous Insects

Half a million species of herbivorous insects have been described. Most of them are diet specialists, using only a few plant species as hosts. Biologists suspect that their specificity is key to their diversity. But why do herbivorous insects tend to be diet specialists? In this review, we catalog a broad range of explanations. We review the evidence for each and suggest lines of research to obtain the evidence we lack. We then draw attention to a second major question, namely how changes in diet breadth affect the rest of a species’ biology. In particular, we know little about how changes in diet breadth feed back on genetic architecture, the population genetic environment, and other aspects of a species’ ecology. Knowing more about how generalists and specialists differ should go a long way toward sorting out potential explanations of specificity, and yield a deeper understanding of herbivorous insect diversity.

The effectiveness of pseudomagic traits in promoting divergence and enhancing local adaptation

"Magic traits," in which the same trait is both under divergent ecological selection and forms the basis of assortative mating, have been sought after due to their supposed unique ability to promote divergence with gene flow. Here, we ask how unique magic traits are, by exploring whether a tightly linked complex of a locus under divergent selection and a locus that acts as a mating cue can mimic a magic trait in its divergence. We find that these "pseudomagic traits" can be very effective in promoting divergence; with tight linkage they are essentially as effective as a magic trait and with loose linkage, and even no linkage, divergence can still be enhanced. Distinguishing between magic and pseudomagic traits in empirical cases may thus not be important when inferring their role in divergence. The ability of divergence in the mating trait to drive divergence in the ecological trait by lowering the effective migration rate, which occurs somewhat even without linkage, is particularly striking; magic traits are typically considered to have the other direction of causality. Our results thus suggest that divergence in a mating trait can at least modestly increase local adaption by allowing more ecological divergence, particularly with tighter linkage.

QTL mapping and GWAS for field kernel water content and kernel dehydration rate before physiological maturity in maize

Kernel water content (KWC) and kernel dehydration rate (KDR) are two main factors affecting maize seed quality and have a decisive influence on the mechanical harvest. It is of great importance to map and mine candidate genes related to KWCs and KDRs before physiological maturity in maize. 120 double-haploid (DH) lines constructed from Si287 with low KWC and JiA512 with high KWC were used as the mapping population. KWCs were measured every 5 days from 10 to 40 days after pollination, and KDRs were calculated. A total of 1702 SNP markers were used to construct a linkage map, with a total length of 1,309.02 cM and an average map distance of 0.77 cM. 10 quantitative trait loci (QTLs) and 27 quantitative trait nucleotides (QTNs) were detected by genome-wide composite interval mapping (GCIM) and multi-locus random-SNP-effect mixed linear model (mrMLM), respectively. One and two QTL hotspot regions were found on Chromosome 3 and 7, respectively. Analysis of the Gene Ontology showed that 2 GO terms of biological processes (BP) were significantly enriched (P ≤ 0.05) and 6 candidate genes were obtained. This study provides theoretical support for marker-assisted breeding of mechanical harvest variety in maize.

Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus

Genomic mapping of the loci associated with phenotypic evolution has revealed genomic "hotspots," or regions of the genome that control multiple phenotypic traits. This clustering of loci has important implications for the speed and maintenance of adaptation and could be due to pleiotropic effects of a single mutation or tight genetic linkage of multiple causative mutations affecting different traits. The threespine stickleback (Gasterosteus aculeatus) is a powerful model for the study of adaptive evolution because the marine ecotype has repeatedly adapted to freshwater environments across the northern hemisphere in the last 12,000 years. Freshwater ecotypes have repeatedly fixed a 16 kilobase haplotype on chromosome IV that contains Ectodysplasin (Eda), a gene known to affect multiple traits, including defensive armor plates, lateral line sensory hair cells, and schooling behavior. Many additional traits have previously been mapped to a larger region of chromosome IV that encompasses the Eda freshwater haplotype. To identify which of these traits specifically map to this adaptive haplotype, we made crosses of rare marine fish heterozygous for the freshwater haplotype in an otherwise marine genetic background. Further, we performed fine-scale association mapping in a fully interbreeding, polymorphic population of freshwater stickleback to disentangle the effects of pleiotropy and linkage on the phenotypes affected by this haplotype. Although we find evidence that linked mutations have small effects on a few phenotypes, a small 1.4-kb region within the first intron of Eda has large effects on three phenotypic traits: lateral plate count, and both the number and patterning of the posterior lateral line neuromasts. Thus, the Eda haplotype is a hotspot of adaptation in stickleback due to both a small, pleiotropic region affecting multiple traits as well as multiple linked mutations affecting additional traits.

Timing alters how a heat shock affects a host-parasitoid interaction

Abiotic factors' effects on species are now well-studied, yet they are still often difficult to predict, especially for strongly interacting species. If these altered abiotic factors and species interactions occur as discrete events in time, such complications may occur because of the events' relative timing. One such discrete abiotic factor is the short-duration, large magnitude increase in temperature called a heat shock. This study investigates how the timing of heat shocks affects the successful attack and reproduction of a parasitoid wasp (Aphidius ervi) attacking its host, the pea aphid (Acyrthosiphon pisum). We tested three relative timings: 1) heat shock before the wasp attacks hosts, 2) heat shock while the wasp is foraging, and 3) heat shock after the wasp has attacked hosts. In each scenario we compared wasp mummy production (pupal stage) with and without a heat shock. Our results showed that a heat shock had the largest effect when it occurred while wasps actively foraged, with fewer mummies produced when exposed to a heat shock compared to the no heat shock control. Follow-up behavioral tests suggest this was caused by wasps becoming inactive during heat shocks. In contrast, when heat shocks were applied three days before or after foraging, we found no difference in mummy production between the heat shock treatment and no heat shock control. These results show the potential importance of timing when considering the ramifications of an altered abiotic factor, especially with relatively discrete abiotic events and interactions.

Intraspecific difference among herbivore lineages and their host-plant specialization drive the strength of trophic cascades

Trophic cascades - the indirect effect of predators on non-adjacent lower trophic levels - are important drivers of the structure and dynamics of ecological communities. However, the influence of intraspecific trait variation on the strength of trophic cascade remains largely unexplored, which limits our understanding of the mechanisms underlying ecological networks. Here we experimentally investigated how intraspecific difference among herbivore lineages specialized on different host plants influences trophic cascade strength in a terrestrial tri-trophic system. We found that the occurrence and strength of the trophic cascade are strongly influenced by herbivores' lineage and host-plant specialization but are not associated with density-dependent effects mediated by the growth rate of herbivore populations. Our findings stress the importance of intraspecific heterogeneities and evolutionary specialization as drivers of trophic cascade strength and underline that intraspecific variation should not be overlooked to decipher the joint influence of evolutionary and ecological factors on the functioning of multi-trophic interactions.

Improved Genome Assembly and Annotation of the Soybean Aphid (Aphis glycines Matsumura)

Aphids are an economically important insect group due to their role as plant disease vectors. Despite this economic impact, genomic resources have only been generated for a small number of aphid species. The soybean aphid (Aphis glycines Matsumura) was the third aphid species to have its genome sequenced and the first to use long-read sequence data. However, version 1 of the soybean aphid genome assembly has low contiguity (contig N50 = 57 Kb, scaffold N50 = 174 Kb), poor representation of conserved genes and the presence of genomic scaffolds likely derived from parasitoid wasp contamination. Here, I use recently developed methods to reassemble the soybean aphid genome. The version 2 genome assembly is highly contiguous, containing half of the genome in only 40 scaffolds (contig N50 = 2.00 Mb, scaffold N50 = 2.51 Mb) and contains 11% more conserved single-copy arthropod genes than version 1. To demonstrate the utility of this improved assembly, I identify a region of conserved synteny between aphids and Drosophila containing members of the Osiris gene family that was split over multiple scaffolds in the original assembly. The improved genome assembly and annotation of A. glycines demonstrates the benefit of applying new methods to old data sets and will provide a useful resource for future comparative genome analysis of aphids.

Genomic signatures of host-associated divergence and adaptation in a coral-eating snail, Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo-Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails' adaptation to their coral hosts. Genome-wide single nucleotide polymorphism (SNP) data from type II restriction site-associated DNA (2b-RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host-associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.

Gene Family Evolution in the Pea Aphid Based on Chromosome-Level Genome Assembly

Genome structural variations, including duplications, deletions, insertions, and inversions, are central in the evolution of eukaryotic genomes. However, structural variations present challenges for high-quality genome assembly, hampering efforts to understand the evolution of gene families and genome architecture. An example is the genome of the pea aphid (Acyrthosiphon pisum) for which the current assembly is composed of thousands of short scaffolds, many of which are known to be misassembled. Here, we present an improved version of the A. pisum genome based on the use of two long-range proximity ligation methods. The new assembly contains four long scaffolds (40-170 Mb), corresponding to the three autosomes and the X chromosome of A. pisum, and encompassing 86% of the new assembly. Assembly accuracy is supported by several quality assessments. Using this assembly, we identify the chromosomal locations and relative ages of duplication events, and the locations of horizontally acquired genes. The improved assembly illuminates the mode of gene family evolution by providing proximity information between paralogs. By estimating nucleotide polymorphism and coverage depth from resequencing data, we determined that many short scaffolds not assembling to chromosomes represent hemizygous regions, which are especially frequent on the highly repetitive X chromosome. Aligning the X-linked aphicarus region, responsible for male wing dimorphism, to the new assembly revealed a 50-kb deletion that cosegregates with the winged male phenotype in some clones. These results show that long-range scaffolding methods can substantially improve assemblies of repetitive genomes and facilitate study of gene family evolution and structural variation.

The Evolution of Immigration Strategies Facilitates Niche Expansion by Divergent Adaptation in a Structured Metapopulation Model

Local adaptation and habitat choice are two key factors that control the distribution and diversification of species. Here we model habitat choice mechanistically as the outcome of dispersal with nonrandom immigration. We consider a structured metapopulation with a continuous distribution of patch types and determine the evolutionarily stable immigration strategy as the function linking patch type to the probability of settling in the patch on encounter. We uncover a novel mechanism whereby coexisting strains that only slightly differ in their local adaptation trait can evolve substantially different immigration strategies. In turn, different habitat use selects for divergent adaptations in the two strains. We propose that the joint evolution of immigration and local adaptation can facilitate diversification and discuss our results in the light of niche conservatism versus niche expansion.

Evolutionary genomics of host plant adaptation: insights from Drosophila

Variation in gene expression in response to the use of alternate host plants can reveal genetic and physiological mechanisms explaining why insect-host relationships vary from host specialism to generalism. Interpreting transcriptome variation relies on well-annotated genomes, making drosophilids valuable model systems, particularly those species with tractable ecological associations. Patterns of whole genome expression and alternate gene splicing in response to growth on different hosts have revealed expression of gene networks of known detoxification genes as well as novel functionally enriched genes of diverse metabolic and structural functions. Integrating trancriptomic responses with fitness differences and levels of phenotypic plasticity in response to alternate hosts will help to reveal the general nature of genotype-phenotype relationships.

Analysis of the extent of synteny and conservation in the gene order in aphids: A first glimpse from the Aphis glycines genome

In the last decade several insect genomes have been sequenced, but for most the chromosomal mapping of the identified scaffolds/annotated genes is not available. The lack of this information makes it difficult to analyse various genetic aspects, including the presence of genome rearrangements and the extent of synteny within and across species. We mapped five multigenic DNA families (major and minor rDNAs, histone gene cluster, esterases and carotenoid desaturases) and seven scaffolds corresponding to 9 Mb of the soybean aphid, Aphis glycines, genome and identified loci spanning the four soybean aphid chromosomes. A comparative analysis of the localization of the annotated A. glycines genes with respect to the peach potato aphid, Myzus persicae, and the fly, Drosophila melanogaster, evidenced a lower degree of synteny between the two aphid species than in the aphid-fly comparison. Only 1.4 genes per syntenic block were observed in aphids in contrast to 2.3 genes per block in flies. This higher chromosomal rearrangement rate in aphids could be explained considering that they possess holocentric chromosomes that can favour the stabilization and inheritance of chromosomal rearrangements. Lastly, our experiments did not detect the presence of chimeric assemblies in the newly available A. glycines biotype 1 genome, differently from what reported in assembled genome of other aphid species, suggesting that chromosomal mapping can be used to ascertain the quality of assembled genomes.

The host range of Aphis gossypii is dependent on aphid genetic background and feeding experience

Background

A polyphagous insect herbivore has a wide range of host plants. However, it has been found that many polyphagous herbivores commonly exhibit a strong preference for a subset of species in their broad host range, and various host biotypes exist in herbivore populations. Nutrition and secondary metabolites in plants affect herbivore preference and performance, but it is still not clear which factors determine the host range and host preference of polyphagous herbivores.

Method

Cotton-melon aphids, Aphis gossypii Glover, collected from cotton and cucumber crops, were used in this study. The genetic backgrounds of these aphids were detected using microsatellite PCR and six genotypes were evaluated. Performance of these six aphid genotypes on excised leaves and plants of cotton and cucumber seedlings were examined through a reciprocal transplant experiment. In order to detect whether the feeding experience on artificial diet would alter aphid host range, the six genotypes of aphids fed on artificial diet for seven days were transferred onto cotton and cucumber leaves, and then their population growth on these two host plants was surveyed.

Results

Aphids from cotton and cucumber plants could not colonize the excised leaves and intact plants of cucumber and cotton seedlings, respectively. All six genotypes of aphids collected from cotton and cucumber plants could survive and produce offspring on artificial diet, which lacked plant secondary metabolites. The feeding experience on the artificial diet did not alter the ability of all six genotypes to use their native host plants. However, after feeding on this artificial diet for seven days, two aphid genotypes from cotton and one from cucumber acquired the ability to use both of the excised leaves from cucumber and cotton plants. The two aphid genotypes from cotton conditioned by the feeding experience on artificial diet and then reared on excised cucumber leaves for >12 generations still maintained the ability to use intact cotton plants but did not establish a population on cucumber plants. However, one cucumber genotype conditioned by artificial diet and then reared on excised cotton leaves could use both the intact cotton and cucumber plants, showing that the expansion of host range was mediated by feeding experience.

Conclusion

Feeding experience on artificial diet induced the expansion of host range of the cucurbit-specialized A. gossypii, and this expansion was genotype-specific. We speculated that feeding on a constant set of host plants in the life cycle of aphids may contribute to the formation of host specialization.

Different phenotypic plastic responses to predators observed among aphid lineages specialized on different host plants

The role of intraspecific variation in the magnitude and direction of plastic responses in ecology and evolution is increasingly recognized. However, the factors underlying intraspecific variation in plastic responses remain largely unexplored, particularly for the hypothesis that the herbivores' phenotypic response to predators might vary amongst lineages associated with different host plants. Here, we tested whether plant-specialized lineages of the pea aphid, Acyrthosiphon pisum, differed in their transgenerational phenotypic response to ladybird predators (i.e., the asexual production of winged offspring by wingless mothers). In a full factorial laboratory experiment, we found that six aphid clonal lineages each specialized either on alfalfa or clover significantly differed in their transgenerational phenotypic response to predators. Some lineages produced an increased number of winged aphids in predator presence while others did not respond. Aphid lineages specialized on alfalfa had stronger phenotypic responses to predators than those specialized on clover. Although we tested only six aphid lineages from two biotypes, our results imply that intraspecific variation in prey phenotypic response of herbivores to predators differs amongst lineages specialized on different host plants. Our findings therefore raise the question of the influence of plant specialization in shaping herbivore phenotypic responses, and highlight the need to consider multi-trophic interactions to understand the causes and consequences of intraspecific variation in complex phenotypic traits.

Inter-chromosomal coupling between vision and pigmentation genes during genomic divergence

Recombination between loci underlying mate choice and ecological traits is a major evolutionary force acting against speciation with gene flow. The evolution of linkage disequilibrium between such loci is therefore a fundamental step in the origin of species. Here, we show that this process can take place in the absence of physical linkage in hamlets-a group of closely related reef fishes from the wider Caribbean that differ essentially in colour pattern and are reproductively isolated through strong visually-based assortative mating. Using full-genome analysis, we identify four narrow genomic intervals that are consistently differentiated among sympatric species in a backdrop of extremely low genomic divergence. These four intervals include genes involved in pigmentation (sox10), axial patterning (hoxc13a), photoreceptor development (casz1) and visual sensitivity (SWS and LWS opsins) that develop islands of long-distance and inter-chromosomal linkage disequilibrium as species diverge. The relatively simple genomic architecture of species differences facilitates the evolution of linkage disequilibrium in the presence of gene flow.

Genetic dissection of assortative mating behavior

The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.

Intrinsic pre-zygotic reproductive isolation of distantly related pea aphid host races

Human activities may weaken or destroy reproductive isolation between young taxa, leading to their fusion with consequences for population and community ecology. Pea aphid host races are adapted to different legume taxa, providing a degree of pre-mating isolation mediated by habitat choice. Yet, all races can feed and reproduce on the broad bean (Vicia faba), a major crop which represents a 'universal host plant', which can promote hybridization between races. Here, we ask if pea aphid host races have reproductive barriers which prevent or reduce gene flow when they co-occur on the universal host plant. We observed mating behaviour, female survival, number of eggs and egg fertilization rates for three types of crosses: among individuals of the same host race, between closely related host races and between distantly related host races. We did not find significant differences in mating behaviour and female survival among the three types of crosses. However, we observed a drastic reduction in the number of eggs laid, and in the number of fertilized eggs, in distant crosses. We conclude that widespread broad bean cultivation in agriculture may predispose closely related-but not distantly related-host races to hybridize, disrupting reproductive isolation between incipient species.

Diet, Gut Microbes and Host Mate Choice: Understanding the significance of microbiome effects on host mate choice requires a case by case evaluation

All organisms live in close association with microbes. However, not all such associations are meaningful in an evolutionary context. Current debate concerns whether hosts and microbes are best described as communities of individuals or as holobionts (selective units of hosts plus their microbes). Recent reports that assortative mating of hosts by diet can be mediated by commensal gut microbes have attracted interest as a potential route to host reproductive isolation (RI). Here, the authors discuss logical problems with this line of argument. The authors briefly review how microbes can affect host mating preferences and evaluate recent findings from fruitflies. Endosymbionts can potentially influence host RI given stable and recurrent co-association of hosts and microbes over evolutionary time. However, observations of co-occurrence of microbes and hosts are ripe for misinterpretation and such associations will rarely represent a meaningful holobiont. A framework in which hosts and their microbes are independent evolutionary units provides the only satisfactory explanation for the observed range of effects and associations.

A Maladaptive Combination of Traits Contributes to the Maintenance of a Drosophila Hybrid Zone

Drosophila teissieri and D. yakuba diverged approximately 3 mya and are thought to share a large, ancestral, African range [1-3]. These species now co-occur in parts of continental Africa and in west Africa on the island of Bioko [1, 4]. While D. yakuba is a human commensal, D. teissieri seems to be associated with Parinari fruits, restricting its range to forests [4-6]. Genome data indicate introgression, despite no evidence of contemporary hybridization. Here we report the discovery of D. yakuba-D. teissieri hybrids at the interface of secondary forests and disturbed, open habitats on Bioko. We demonstrate that hybrids are the F₁ progeny of D. yakuba females and D. teissieri males. At high temperatures like those found on Bioko, D. teissieri females are generally less receptive to mating, and in combination with temperature effects on egg lay and egg-to-adult viability, this decreases the potential for gene flow between female D. teissieri and male D. yakuba relative to the reciprocal cross. Field and laboratory experiments demonstrate that F₁ hybrids have a maladaptive combination of D. yakuba behavior and D. teissieri physiology, generating additional barriers to gene flow. Nevertheless, analysis of introgressed and non-introgressed regions of the genome indicate that, while rare, gene flow is relatively recent. Our observations identify precise intrinsic and extrinsic factors that, along with hybrid male sterility, limit gene flow and maintain these species. These data contribute to a growing body of literature that suggests the Gulf of Guinea may be a hotspot for hybridization.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.

Exposure to predators does not lead to the evolution of larger brains in experimental populations of threespine stickleback

Natural selection is often invoked to explain differences in brain size among vertebrates. However, the particular agents of selection that shape brain size variation remain obscure. Recent studies suggest that predators may select for larger brains because increased cognitive and sensory abilities allow prey to better elude predators. Yet, there is little direct evidence that exposure to predators causes the evolution of larger brains in prey species. We experimentally tested this prediction by exposing families of 1000-2000 F2 hybrid benthic-limnetic threespine stickleback to predators under naturalistic conditions, along with matched controls. After two generations of selection, we found that fish from the predator addition treatment had significantly smaller brains (specifically smaller telencephalons and optic lobes) than fish from the control treatment. After an additional generation of selection, we reared experimental fish in a common environment and found that this difference in brain size was maintained in the offspring of fish from the predator addition treatment. Our results provide direct experimental evidence that (a) predators can indeed drive the evolution of brain size--but not in the fashion commonly expected and (b) that the tools of experimental evolution can be used to the study the evolution of the vertebrate brain.

The Limits to Parapatric Speciation II: Strengthening a Preexisting Genetic Barrier to Gene Flow in Parapatry

By encompassing the whole continuum between allopatric and sympatric scenarios, parapatric speciation includes many potential scenarios for the evolution of new species. Here, we investigate how a genetic barrier to gene flow, that relies on a single postzygotic genetic incompatibility, may further evolve under ongoing migration. We consider a continent island model with three loci involved in pairwise Dobzhansky–Muller incompatibilities (DMIs). Using an analytic approach, we derive the conditions for invasion of a new mutation and its consequences for the strength and stability of the initial genetic barrier. Our results show that the accumulation of genetic incompatibilities in the presence of gene flow is under strong selective constraints. In particular, preexisting incompatibilities do not always facilitate the invasion of further barrier genes. If new mutations do invade, they will often weaken or destroy the barrier rather than strengthening it. We conclude that migration is highly effective at disrupting the so-called “snowball effect”, the accelerated accumulation of DMIs that has been described for allopatric populations en route to reproductive isolation.

Evolution and Genetics of Precocious Burrowing Behavior in Peromyscus Mice

A central challenge in biology is to understand how innate behaviors evolve between closely related species. One way to elucidate how differences arise is to compare the development of behavior in species with distinct adult traits [1]. Here, we report that Peromyscus polionotus is strikingly precocious with regard to burrowing behavior, but not other behaviors, compared to its sister species P. maniculatus. In P. polionotus, burrows were excavated as early as 17 days of age, whereas P. maniculatus did not build burrows until 10 days later. Moreover, the well-known differences in burrow architecture between adults of these species-P. polionotus adults excavate long burrows with an escape tunnel, whereas P. maniculatus dig short, single-tunnel burrows [2-4]-were intact in juvenile burrowers. To test whether this juvenile behavior is influenced by early-life environment, we reciprocally cross-fostered pups of both species. Fostering did not alter the characteristic burrowing behavior of either species, suggesting that these differences are genetic. In backcross hybrids, we show that precocious burrowing and adult tunnel length are genetically correlated and that a P. polionotus allele linked to tunnel length variation in adults is also associated with precocious onset of burrowing in juveniles, suggesting that the same genetic region-either a single gene with pleiotropic effects or linked genes-influences distinct aspects of the same behavior at these two life stages. These results raise the possibility that genetic variants affect behavioral drive (i.e., motivation) to burrow and thereby affect both the developmental timing and adult expression of burrowing behavior.

Genomics of extreme ecological specialists: multiple convergent evolution but no genetic divergence between ecotypes of Maculinea alcon butterflies

Biotic interactions are often acknowledged as catalysers of genetic divergence and eventual explanation of processes driving species richness. We address the question, whether extreme ecological specialization is always associated with lineage sorting, by analysing polymorphisms in morphologically similar ecotypes of the myrmecophilous butterfly Maculinea alcon. The ecotypes occur in either hygric or xeric habitats, use different larval host plants and ant species, but no significant distinctive molecular traits have been revealed so far. We apply genome-wide RAD-sequencing to specimens originating from both habitats across Europe in order to get a view of the potential evolutionary processes at work. Our results confirm that genetic variation is mainly structured geographically but not ecologically - specimens from close localities are more related to each other than populations of each ecotype from distant localities. However, we found two loci for which the association with xeric versus hygric habitats is supported by segregating alleles, suggesting convergent evolution of habitat preference. Thus, ecological divergence between the forms probably does not represent an early stage of speciation, but may result from independent recurring adaptations involving few genes. We discuss the implications of these results for conservation and suggest preserving biotic interactions and main genetic clusters.

Genetic Coupling of Female Mate Choice with Polygenic Ecological Divergence Facilitates Stickleback Speciation

Ecological speciation with gene flow is widespread in nature [1], but it presents a conundrum: how are associations between traits under divergent natural selection and traits that contribute to assortative mating maintained? Theoretical models suggest that genetic mechanisms inhibiting free recombination between loci underlying these two types of traits (hereafter, "genetic coupling") can facilitate speciation [2-4]. Here, we perform a direct test for genetic coupling by mapping both divergent traits and female mate choice in a classic model of ecological speciation: sympatric benthic and limnetic threespine stickleback (Gasterosteus aculeatus). By measuring mate choice in F2 hybrid females, we allowed for recombination between loci underlying assortative mating and those under divergent ecological selection. In semi-natural mating arenas in which females had access to both benthic and limnetic males, we found that F2 females mated with males similar to themselves in body size and shape. In addition, we found two quantitative trait loci (QTLs) associated with female mate choice that also predicted female morphology along the benthic-limnetic trait axis. Furthermore, a polygenic genetic model that explains adaptation to contrasting benthic and limnetic feeding niches [5] also predicted F2 female mate choice. Together, these results provide empirical evidence that genetic coupling of assortative mating with traits under divergent ecological selection helps maintain species in the face of gene flow, despite a polygenic basis for adaptation to divergent environments.

Significant Synteny and Colocalization of Ecologically Relevant Quantitative Trait Loci Within and Across Species of Salmonid Fishes

The organization of functional regions within genomes has important implications for evolutionary potential. Considerable research effort has gone toward identifying the genomic basis of phenotypic traits of interest through quantitative trait loci (QTL) analyses. Less research has assessed the arrangement of QTL in the genome within and across species. To investigate the distribution, extent of colocalization, and the synteny of QTL for ecologically relevant traits, we used a comparative genomic mapping approach within and across a range of salmonid species. We compiled 943 QTL from all available species [lake whitefish (Coregonus clupeaformis), coho salmon (Oncorhynchus kisutch), rainbow trout (O. mykiss), Chinook salmon (O. tshawytscha), Atlantic salmon (Salmo salar), and Arctic charr (Salvelinus alpinus)]. We developed a novel analytical framework for mapping and testing the distribution of these QTL. We found no correlation between QTL density and gene density at the chromosome level but did at the fine-scale. Two chromosomes were significantly enriched for QTL. We found multiple synteny blocks for morphological, life history, and physiological traits across species, but only morphology and physiology had significantly more than expected. Two or three pairs of traits were significantly colocalized in three species (lake whitefish, coho salmon, and rainbow trout). Colocalization and fine-scale synteny suggest genetic linkage between traits within species and a conserved genetic basis across species. However, this pattern was weak overall, with colocalization and synteny being relatively rare. These findings advance our understanding of the role of genomic organization in the renowned ecological and phenotypic variability of salmonid fishes.