SELECTION AND GENOMIC DIFFERENTIATION DURING ECOLOGICAL SPECIATION: ISOLATING THE CONTRIBUTIONS OF HOST ASSOCIATION VIA A COMPARATIVE GENOME SCAN OF NEOCHLAMISUS BEBBIANAE LEAF BEETLES

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.

Why are there not more herbivorous insect species?

Insect species richness is estimated to exceed three million species, of which roughly half is herbivorous. Despite the vast number of species and varied life histories, the proportion of herbivorous species among plant-consuming organisms is lower than it could be due to constraints that impose limits to their diversification. These include ecological factors, such as vague interspecific competition; anatomical and physiological limits, such as neural limits and inability of handling a wide range of plant allelochemicals; phylogenetic constraints, like niche conservatism; and most importantly, a low level of concerted genetic variation necessary to a phyletic conversion. It is suggested that diversification ultimately depends on what we call the intrinsic trend of diversification of the insect genome. In support of the above, we survey the major types of host-specificity, the mechanisms and constraints of host specialization, possible pathways of speciation, and hypotheses concerning insect diversification.

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.

Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control

Cryptic taxa have often been observed in the form of host‐associated species that diverged as the result of adaptation to alternate host plants. Untangling cryptic diversity in species complexes that encompass invasive species is a mandatory task for pest management. Moreover, investigating the evolutionary history of a species complex may help to understand the drivers of their diversification. The mealybug Hypogeococcus pungens was believed to be a polyphagous species from South America and has been reported as a pest devastating native cacti in Puerto Rico, also threatening cactus diversity in the Caribbean and North America. There is neither certainty about the identity of the pest nor the source population from South America. Recent studies pointed to substantial genetic differentiation among local populations, suggesting that H. pungens is a species complex. In this study, we used a combination of genome‐wide SNPs and mtDNA variation to investigate species diversity within H. pungens sensu lato to establish host plant ranges of each one of the putative members of the complex, to evaluate whether the pattern of host plant association drove diversification in the species complex, and to determine the source population of the Puerto Rican cactus pest. Our results suggested that H. pungens comprises at least five different species, each one strongly associated with specific host plants. We also established that the Puerto Rican cactus pest derives from southeastern Brazilian mealybugs. This is an important achievement because it will help to design reliable strategies for biological control using natural enemies of the pest from its native range.

Machine-learning-based detection of adaptive divergence of the stream mayfly Ephemera strigata populations

Adaptive divergence is a key mechanism shaping the genetic variation of natural populations. A central question linking ecology with evolutionary biology is how spatial environmental heterogeneity can lead to adaptive divergence among local populations within a species. In this study, using a genome scan approach to detect candidate loci under selection, we examined adaptive divergence of the stream mayfly Ephemera strigata in the Natori River Basin in northeastern Japan. We applied a new machine-learning method (i.e., random forest) besides traditional distance-based redundancy analysis (dbRDA) to examine relationships between environmental factors and adaptive divergence at non-neutral loci. Spatial autocorrelation analysis based on neutral loci was employed to examine the dispersal ability of this species. We conclude the following: (a) E. strigata show altitudinal adaptive divergence among the populations in the Natori River Basin; (b) random forest showed higher resolution for detecting adaptive divergence than traditional statistical analysis; and (c) separating all markers into neutral and non-neutral loci could provide full insight into parameters such as genetic diversity, local adaptation, and dispersal ability.

Divergence with gene flow is driven by local adaptation to temperature and soil phosphorus concentration in teosinte subspecies (Zea mays parviglumis and Zea mays mexicana)

Patterns of genomic divergence between hybridizing taxa can be heterogeneous along the genome. Both differential introgression and local adaptation may contribute to this pattern. Here, we analysed two teosinte subspecies, Zea mays ssp. parviglumis and ssp. mexicana, to test whether their divergence has occurred in the face of gene flow and to infer which environmental variables have been important drivers of their ecological differentiation. We generated 9,780 DArTseqTM SNPs for 47 populations, and used an additional data set containing 33,454 MaizeSNP50 SNPs for 49 populations. With these data, we inferred features of demographic history and performed genome wide scans to determine the number of outlier SNPs associated with climate and soil variables. The two data sets indicate that divergence has occurred or been maintained despite continuous gene flow and/or secondary contact. Most of the significant SNP associations were to temperature and to phosphorus concentration in the soil. A large proportion of these candidate SNPs were located in regions of high differentiation that had been identified previously as putative inversions. We therefore propose that genomic differentiation in teosintes has occurred by a process of adaptive divergence, with putative inversions contributing to reduced gene flow between locally adapted populations.

First molecular phylogeny of Agrilus (Coleoptera: Buprestidae), the largest genus on Earth, with DNA barcode database for forestry pest diagnostics

All more than 3000 species of Agrilus beetles are phytophagous and some cause economically significant damage to trees and shrubs. Facilitated by international trade, Agrilus species regularly invade new countries and continents. This necessitates a rapid identification of Agrilus species, as the first step for subsequent protective measures. This study provides the first DNA reference library for ~100 Agrilus species from the Northern Hemisphere based on three mitochondrial markers: cox1-5' (DNA barcode fragment), cox1-3', and rrnL. All 329 Agrilus records available in the Barcode of Life Database format, including specimen images and geo data, are released through a public dataset 'Agrilus1 329' available at: dx.doi.org/10.5883/DS-AGRILUS1. All Agrilus species were identified using adult morphology and by using molecular phylogenetic trees, as well as distance- and tree-based algorithms. Most DNA-based species limits agree well with the morphology-based identification. Our results include cases of high intraspecific variability and multiple species para- and polyphyly. DNA barcoding is a powerful species identification tool in Agrilus, although it frequently fails to recover morphologically-delimited Agrilus species-group. Even though the current three-gene database covers only ~3% of the known Agrilus diversity, it contains representatives of all principal lineages from the Northern Hemisphere and represents the most extensive dataset built for DNA-delimited species identification within this genus so far. Molecular data analyses can rapidly and cost-effectively identify an unknown sample, including immature stages and/or non-native taxa, or species not yet formally named.

Small and genetically highly structured populations in a long-legged bee, Rediviva longimanus, as inferred by pooled RAD-seq

Adaptation to local host plants may impact a pollinator's population genetic structure by reducing gene flow and driving population genetic differentiation, representing an early stage of ecological speciation. South African Rediviva longimanus bees exhibit elongated forelegs, a bizarre adaptation for collecting oil from floral spurs of their Diascia hosts. Furthermore, R. longimanus foreleg length (FLL) differs significantly among populations, which has been hypothesised to result from selection imposed by inter-population variation in Diascia floral spur length. Here, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to investigate the population genetic structure of R. longimanus and to test if phenotypic differences in FLL translate into increased genetic differentiation (i) between R. longimanus populations and (ii) between phenotypes across populations. We also inferred the effects of demographic processes on population genetic structure and tested for genetic markers underpinning local adaptation. RESULTS: Populations showed marked genetic differentiation (average FST = 0.165), though differentiation was not statistically associated with differences between populations in FLL. All populations exhibited very low genetic diversity and were inferred to have gone through recent bottleneck events, suggesting extremely low effective population sizes. Genetic differentiation between samples pooled by leg length (short versus long) rather than by population of origin was even higher (FST = 0.260) than between populations, suggesting reduced interbreeding between long and short-legged individuals. Signatures of selection were detected in 1119 (3.8%) of a total of 29,721 SNP markers, CONCLUSIONS: Populations of R. longimanus appear to be small, bottlenecked and isolated. Though we could not detect the effect of local adaptation (FLL in response to floral spurs of host plants) on population genetic differentiation, short and long legged bees appeared to be partially differentiated, suggesting incipient ecological speciation. To test this hypothesis, greater resolution through the use of individual-based whole-genome analyses is now needed to quantify the degree of reproductive isolation between long and short legged bees between and even within populations.

Genomic islands of divergence in the Yellow Tang and the Brushtail Tang Surgeonfishes

The current ease of obtaining thousands of molecular markers challenges the notion that full phylogenetic concordance, as proposed by phylogenetic species concepts, is a requirement for defining species delimitations. Indeed, the presence of genomic islands of divergence, which may be the cause, or in some cases the consequence, of speciation, precludes concordance. Here, we explore this issue using thousands of RAD markers on two sister species of surgeonfishes (Teleostei: Acanthuridae), Zebrasoma flavescens and Z. scopas, and several populations within each species. Species are readily distinguished based on their colors (solid yellow and solid brown, respectively), yet populations and species are neither distinguishable using mitochondrial markers (cytochrome c oxidase 1), nor using 5193 SNPs (pairwise Φst = 0.034). In contrast, when using outlier loci, some of them presumably under selection, species delimitations, and strong population structure follow recognized taxonomic positions (pairwise Φst = 0.326). Species and population delimitation differences based on neutral and selected markers are likely due to local adaptation, thus being consistent with the idea that these genomic islands of divergence arose as a consequence of isolation. These findings, which are not unique, raise the question of a potentially important pathway of divergence based on local adaptation that is only evident when looking at thousands of loci.

Heritable gene expression differences between lake and stream stickleback include both parallel and antiparallel components

The repeated phenotypic patterns that characterize populations undergoing parallel evolution provide support for a deterministic role of adaptation by natural selection. Determining the level of parallelism also at the genetic level is thus central to our understanding of how natural selection works. Many studies have looked for repeated genomic patterns in natural populations, but work on gene expression is less common. The studies that have examined gene expression have found some support for parallelism, but those studies almost always used samples collected from the wild that potentially confounds the effects of plasticity with heritable differences. Here we use two independent pairs of lake and stream threespine stickleback (Gasterosteus aculeatus) raised in common garden conditions to assess both parallel and antiparallel (that is, similar versus different directions of lake–stream expression divergence in the two watersheds) heritable gene expression differences as measured by total RNA sequencing. We find that more genes than expected by chance show either parallel (22 genes, 0.18% of expressed genes) or antiparallel (24 genes, 0.20% of expressed genes) lake–stream expression differences. These results correspond well with previous genomic studies in stickleback ecotype pairs that found similar levels of parallelism. We suggest that parallelism might be similarly constrained at the genomic and transcriptomic levels.

Parallel Patterns of Host-Specific Morphology and Genetic Admixture in Sister Lineages of a Commensal Barnacle

Symbiotic relationships are often species specific, allowing symbionts to adapt to their host environments. Host generalists, on the other hand, have to cope with diverse environments. One coping strategy is phenotypic plasticity, defined by the presence of host-specific phenotypes in the absence of genetic differentiation. Recent work indicates that such host-specific phenotypic plasticity is present in the West Pacific lineage of the commensal barnacle Chelonibia testudinaria (Linnaeus, 1758). We investigated genetic and morphological host-specific structure in the genetically distinct Atlantic sister lineage of C. testudinaria. We collected adult C. testudinaria from loggerhead sea turtles, horseshoe crabs, and blue crabs along the eastern U.S. coast between Delaware and Florida and in the Gulf of Mexico off Mississippi. We find that shell morphology, especially shell thickness, is host specific and comparable in similar host species between the Atlantic and West Pacific lineages. We did not detect significant genetic differentiation related to host species when analyzing data from 11 nuclear microsatellite loci and mitochondrial sequence data, which is comparable to findings for the Pacific lineage. The most parsimonious explanation for these parallel patterns between distinct lineages of C. testudinaria is that C. testudinaria maintained phenotypic plasticity since the lineages diverged 4-5 mya.

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow - a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co-occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment-length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host-associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.

Evidence for ecological speciation via a host shift in the holly leaf miner, Phytomyza glabricola (Diptera: Agromyzidae)

Evolutionary radiations have been well documented in plants and insects, and natural selection may often underly these radiations. If radiations are adaptive, the diversity of species could be due to ecological speciation in these lineages. Agromyzid flies exhibit patterns of repeated host-associated radiations. We investigated whether host-associated population divergence and evidence of divergent selection exist in the leaf miner Phytomyza glabricola on its sympatric host plants, the holly species, Ilex coriacea and I. glabra. Using AFLPs and nuclear sequence data, we found substantial genetic divergence between host-associated populations of these flies throughout their geographic range. Genome scans using the AFLP data identified 13 loci under divergent selection, consistent with processes of ecological speciation. EF-1α data suggest that I. glabra is the original host of P. glabricola and that I. coriacea is the novel host, but the AFLP data are ambiguous with regard to directionality of the host shift.

Patterns of Genome-Wide Variation in Glossina fuscipes fuscipes Tsetse Flies from Uganda

The tsetse fly Glossina fuscipes fuscipes (Gff) is the insect vector of the two forms of Human African Trypanosomiasis (HAT) that exist in Uganda. Understanding Gff population dynamics, and the underlying genetics of epidemiologically relevant phenotypes is key to reducing disease transmission. Using ddRAD sequence technology, complemented with whole-genome sequencing, we developed a panel of ∼73,000 single-nucleotide polymorphisms (SNPs) distributed across the Gff genome that can be used for population genomics and to perform genome-wide-association studies. We used these markers to estimate genomic patterns of linkage disequilibrium (LD) in Gff, and used the information, in combination with outlier-locus detection tests, to identify candidate regions of the genome under selection. LD in individual populations decays to half of its maximum value (r(2) max/2) between 1359 and 2429 bp. The overall LD estimated for the species reaches r(2) max/2 at 708 bp, an order of magnitude slower than in Drosophila Using 53 infected (Trypanosoma spp.) and uninfected flies from four genetically distinct Ugandan populations adapted to different environmental conditions, we were able to identify SNPs associated with the infection status of the fly and local environmental adaptation. The extent of LD in Gff likely facilitated the detection of loci under selection, despite the small sample size. Furthermore, it is probable that LD in the regions identified is much higher than the average genomic LD due to strong selection. Our results show that even modest sample sizes can reveal significant genetic associations in this species, which has implications for future studies given the difficulties of collecting field specimens with contrasting phenotypes for association analysis.

Recurrent specialization on a toxic fruit in an island Drosophila population

Recurrent specialization on similar host plants offers a unique opportunity to unravel the evolutionary and genetic mechanisms underlying dietary shifts. Recent studies have focused on ecological races belonging to the same species, but it is hard in many cases to untangle the role of adaptive introgression versus distinct mutations in facilitating recurrent evolution. We discovered on the island of Mayotte a population of the generalist fly Drosophila yakuba that is strictly associated with noni (Morinda citrifolia). This case strongly resembles Drosophila sechellia, a genetically isolated insular relative of D. yakuba whose intensely studied specialization on toxic noni fruits has always been considered a unique event in insect evolution. Experiments revealed that unlike mainland D. yakuba strains, Mayotte flies showed strong olfactory attraction and significant toxin tolerance to noni. Island females strongly discriminated against mainland males, suggesting that dietary adaptation has been accompanied by partial reproductive isolation. Population genomic analysis indicated a recent colonization (∼29 kya), at a time when year-round noni fruits may have presented a predictable resource on the small island, with ongoing migration after colonization. This relatively recent time scale allowed us to search for putatively adaptive loci based on genetic variation. Strong signals of genetic differentiation were found for several detoxification genes, including a major toxin tolerance locus in D. sechellia Our results suggest that recurrent evolution on a toxic resource can involve similar historical events and common genetic bases, and they establish an important genetic system for the study of early stages of ecological specialization and speciation.

Genomic scanning using AFLP to detect loci under selection in the moss Funaria hygrometrica along a climate gradient in the Sierra Nevada Mountains, Spain

The common cord moss Funaria hygrometrica has a worldwide distribution and thrives in a wide variety of environments. Here, we studied the genetic diversity in F. hygrometrica along an abiotic gradient in the Mediterranean high mountain of Sierra Nevada (Spain) using a genome scan method. Eighty-four samples from 17 locations from 24 to 2700 m were fingerprinted based on their amplified fragment length polymorphism (AFLP) banding pattern. Using PCA and Bayesian inference we found that the genetic diversity was structured in three or four clusters, respectively. Using a genome scan method we identified 13 outlier loci, which showed a signature of positive selection. Partial Mantel tests were performed between the Euclidean distance matrices of geographic and climatic variables, versus the pair-wise genetic distance of the AFLP dataset and AFLP-positive outliers dataset. AFLP-positive outlier data were significantly correlated with the gradient of the climatic variables, suggesting adaptive variation among populations of F. hygrometrica along the Sierra Nevada Mountains. We highlight the additional analyses necessary to identify the nature of these loci, and their biological role in the adaptation process.

Strong premating reproductive isolation drives incipient speciation in Mimulus aurantiacus

Determining which forms of reproductive isolation have the biggest impact on the process of divergence is a major goal of speciation research. These barriers are often divided into those that affect the potential for hybridization (premating isolation), and those that occur after mating (postmating isolation), and much debate has surrounded the relative importance of these categories. Within the species Mimulus aurantiacus, red- and yellow-flowered ecotypes occur in the southwest corner of California, and a hybrid zone occurs where their ranges overlap. We show that premating barriers are exclusively responsible for isolation in this system, with both ecogeographic and pollinator isolation contributing significantly to total isolation. Postmating forms of reproductive isolation have little or no impact on gene flow, indicating that hybrids likely contribute to introgression at neutral loci. Analysis of molecular variation across thousands of restriction-site associated DNA sequencing (RAD-seq) markers reveals that the genomes of these taxa are largely undifferentiated. However, structure analysis shows that these taxa are distinguishable genetically, likely due to the impact of loci underlying differentiated adaptive phenotypes. These data exhibit the power of divergent natural selection to maintain highly differentiated phenotypes in the face of gene flow during the early stages of speciation.

The roles of geography and founder effects in promoting host-associated differentiation in the generalist bogus yucca moth Prodoxus decipiens

There is ample evidence that host shifts in plant-feeding insects have been instrumental in generating the enormous diversity of insects. Changes in host use can cause host-associated differentiation (HAD) among populations that may lead to reproductive isolation and eventual speciation. The importance of geography in facilitating this process remains controversial. We examined the geographic context of HAD in the wide-ranging generalist yucca moth Prodoxus decipiens. Previous work demonstrated HAD among sympatric moth populations feeding on two different Yucca species occurring on the barrier islands of North Carolina, USA. We assessed the genetic structure of P. decipiens across its entire geographic and host range to determine whether HAD is widespread in this generalist herbivore. Population genetic analyses of microsatellite and mtDNA sequence data across the entire range showed genetic structuring with respect to host use and geography. In particular, genetic differentiation was relatively strong between mainland populations and those on the barrier islands of North Carolina. Finer scale analyses, however, among sympatric populations using different host plant species only showed significant clustering based on host use for populations on the barrier islands. Mainland populations did not form population clusters based on host plant use. Reduced genetic diversity in the barrier island populations, especially on the derived host, suggests that founder effects may have been instrumental in facilitating HAD. In general, results suggest that the interplay of local adaptation, geography and demography can determine the tempo of HAD. We argue that future studies should include comprehensive surveys across a wide range of environmental and geographic conditions to elucidate the contribution of various processes to HAD.

Widespread signals of convergent adaptation to high altitude in Asia and america

Living at high altitude is one of the most difficult challenges that humans had to cope with during their evolution. Whereas several genomic studies have revealed some of the genetic bases of adaptations in Tibetan, Andean, and Ethiopian populations, relatively little evidence of convergent evolution to altitude in different continents has accumulated. This lack of evidence can be due to truly different evolutionary responses, but it can also be due to the low power of former studies that have mainly focused on populations from a single geographical region or performed separate analyses on multiple pairs of populations to avoid problems linked to shared histories between some populations. We introduce here a hierarchical Bayesian method to detect local adaptation that can deal with complex demographic histories. Our method can identify selection occurring at different scales, as well as convergent adaptation in different regions. We apply our approach to the analysis of a large SNP data set from low- and high-altitude human populations from America and Asia. The simultaneous analysis of these two geographic areas allows us to identify several candidate genome regions for altitudinal selection, and we show that convergent evolution among continents has been quite common. In addition to identifying several genes and biological processes involved in high-altitude adaptation, we identify two specific biological pathways that could have evolved in both continents to counter toxic effects induced by hypoxia.

Genetic analysis of population differentiation and adaptation in Leuciscus waleckii

Demographic events and natural selection both influence animal phenotypic and genetic variation; exploring the effects of demography and selection on population divergence is of great significance in evolutionary biology. To uncover the causes behind the patterns of genetic differentiation and adaptation among six populations of Leuciscus waleckii from Dali Basin (two populations, alkaline vs. freshwater) and Amur Basin (four populations, freshwater rivers vs. alkaline lake), a set of 21 unlinked polymorphic microsatellite markers and two mitochondrial DNA sequences (Cytb and D-loop) were applied to examine whether populations from different environments or habitats have distinct genetic differentiation and whether alkalinity is the major factor that caused population divergence. Bayesian analysis and principal component analysis as well as haplotype network analysis showed that these populations are primarily divided into two groups, which are congruent with geographic separation but not inconsistent with the habitat environment (alkalinity). Using three different approaches, outlier detection indicated that one locus, HLJYL017, may be under directional selection and involved in local adaptation processes. Overall, this study suggested that demographic events and selection of local environmental conditions including of alkalinity are jointly responsible for population divergence. These findings constitute an important step towards the understanding of the genetic basis of differentiation and adaptation, as well as towards the conservation of L. waleckii.

Reproductive isolation between host races of Phytomyza glabricola on Ilex coriacea and I. glabra

Recently diverged taxa often show discordance in genetic divergence among genomic loci, where some loci show strong divergence and others show none at all. Genetic studies alone cannot distinguish among the possible mechanisms but experimental studies on other aspects of divergence may provide guidance in the inference of causes of observed discordances. In this study, we used no-choice mating trials to test for the presence of reproductive isolation between host races of the leaf-mining fly, Phytomyza glabricola on its two holly host species, Ilex coriacea and I. glabra. These trials inform our effort to determine the cause of significant differences in the degree of divergence of nuclear and mitochondrial loci of flies collected from the two host plants. We present evidence of reproductive isolation between host races in a controlled greenhouse setting: significantly more mate pairs consisting of flies from the same host plant species produced offspring than inter-host mate pairs, which produced no offspring. We also tested whether the presence of the natal or non-natal host plant affects reproductive success. Flies collected from I. coriacea were more likely to produce offspring when in the presence of the natal host, whereas the presence or absence of either the natal or non-natal host had no effect on flies collected from I. glabra. The results indicate discordant patterns of nuclear and mitochondrial divergence among host races of P. glabricola are likely due to incomplete lineage sorting, and the host races may be well on their way to becoming biological species.

Natural selection and neutral evolution jointly drive population divergence between alpine and lowland ecotypes of the allopolyploid plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074-0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (FST  = 0.041-0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.

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.

Niche differentiation and colonization of a novel environment by an asexual parasitic wasp

How do asexual taxa become adapted to a diversity of environments, and how do they persist despite changing environmental conditions? These questions are linked by their mutual focus on the relationship between genetic variation, which is often limited in asexuals, and the ability to respond to environmental variation. Asexual taxa originating from a single ancestor present a unique opportunity to assess rates of phenotypic and genetic change when access to new genetic variation is limited to mutation. Diachasma muliebre is an asexual Hymenopteran wasp that is geographically and genetically isolated from all sexual relatives. D. muliebre attack larvae of the western cherry fruit fly (Rhagoletis indifferens), which in turn feed inside bitter cherry fruit (Prunus emarginata) in August and September. R. indifferens has recently colonized a new host plant with an earlier fruiting phenology (June/July), domesticated sweet cherries (P. avium), and D. muliebre has followed its host into this temporally earlier niche. We tested three hypotheses: 1) that all D. muliebre lineages originate from a single asexual ancestor; 2) that different D. muliebre lineages (as defined by unique mtDNA haplotypes) have differentiated on their ancestral host in an important life-history trait, eclosion timing; and 3) that early-eclosing lineages have preferentially colonized the new sweet cherry niche. We find that mitochondrial COI and microsatellite data provide strong support for a single ancestral origin for all lineages. Furthermore, COI sequencing revealed five mitochondrial haplotypes among D. muliebre, and individual wasps possessing one distinctive mitochondrial haplotype (haplotype II) eclosed as reproductive adults significantly earlier than wasps with all other haplotypes. In addition, this early-eclosing lineage of D. muliebre is one of two lineages that have colonized the P. avium habitat, consistent with the preferential colonization hypothesis. These data suggest that D. muliebre has evolved adaptive phenotypic variation despite limited genetic variation, and that this variation has subsequently allowed an expansion of some wasps into a novel habitat. The D. muliebre system may allow for in-depth study of adaptation and long-term persistence of asexual taxa.

Detecting the footprints of divergent selection in oaks with linked markers

Genome scans are increasingly used to study ecological speciation, providing a useful genome-wide perspective on divergent selection in the presence of gene flow. Here, we compare current approaches to detect footprints of divergent selection in closely related species. We analyzed 192 individuals from two interfertile European temperate oak species using 30 nuclear microsatellites from eight linkage groups. These markers present little intraspecific differentiation and can be used in combination to assign individual genotypes to species. We first show that different outlier detection tests give somewhat different results, possibly due to model constraints. Second, using linkage information for these markers, we further characterize the signature of divergent selection in the presence of gene flow. In particular, we show that recombination estimates for regions with outlier markers are lower than those for a control region, in line with a prediction from ecological speciation theory. Most importantly, we show that analyses at the haplotype level can distinguish between truly divergent (bi-directional) selection and positive selection in one of the two species, offering a new and improved method for characterizing the speciation process.

Automated masking of AFLP markers improves reliability of phylogenetic analyses

The amplified fragment length polymorphisms (AFLP) method has become an attractive tool in phylogenetics due to the ease with which large numbers of characters can be generated. In contrast to sequence-based phylogenetic approaches, AFLP data consist of anonymous multilocus markers. However, potential artificial amplifications or amplification failures of fragments contained in the AFLP data set will reduce AFLP reliability especially in phylogenetic inferences. In the present study, we introduce a new automated scoring approach, called “AMARE” (AFLP MAtrix REduction). The approach is based on replicates and makes marker selection dependent on marker reproducibility to control for scoring errors. To demonstrate the effectiveness of our approach we record error rate estimations, resolution scores, PCoA and stemminess calculations. As in general the true tree (i.e. the species phylogeny) is not known, we tested AMARE with empirical, already published AFLP data sets, and compared tree topologies of different AMARE generated character matrices to existing phylogenetic trees and/or other independent sources such as morphological and geographical data. It turns out that the selection of masked character matrices with highest resolution scores gave similar or even better phylogenetic results than the original AFLP data sets.

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.

Clonality despite sex: the evolution of host-associated sexual neighborhoods in the pathogenic fungus Penicillium marneffei

Molecular genetic approaches typically detect recombination in microbes regardless of assumed asexuality. However, genetic data have shown the AIDS-associated pathogen Penicillium marneffei to have extensive spatial genetic structure at local and regional scales, and although there has been some genetic evidence that a sexual cycle is possible, this haploid fungus is thought to be genetically, as well as morphologically, asexual in nature because of its highly clonal population structure. Here we use comparative genomics, experimental mixed-genotype infections, and population genetic data to elucidate the role of recombination in natural populations of P. marneffei. Genome wide comparisons reveal that all the genes required for meiosis are present in P. marneffei, mating type genes are arranged in a similar manner to that found in other heterothallic fungi, and there is evidence of a putatively meiosis-specific mutational process. Experiments suggest that recombination between isolates of compatible mating types may occur during mammal infection. Population genetic data from 34 isolates from bamboo rats in India, Thailand and Vietnam, and 273 isolates from humans in China, India, Thailand, and Vietnam show that recombination is most likely to occur across spatially and genetically limited distances in natural populations resulting in highly clonal population structure yet sexually reproducing populations. Predicted distributions of three different spatial genetic clusters within P. marneffei overlap with three different bamboo rat host distributions suggesting that recombination within hosts may act to maintain population barriers within P. marneffei.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.

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.

Two methods to easily obtain nucleotide sequences from AFLP loci of interest

Genome scans based on anonymous Amplified Fragment Length Polymorphism (AFLP) markers scattered throughout the genome are becoming an increasingly popular approach to study the genetic basis of adaptation and speciation in natural populations. A shortcoming of this approach is that despite its efficiency to detect signatures of selection, it can hardly help pinpoint the specific genomic region(s), gene(s), or mutation(s) targeted by selection. Here, we present two methods to be undertaken after performing an AFLP-based genome scan to easily obtain the sequences of AFLP loci detected as outliers by population genomics approaches. The first one is based on the gel excision of the target AFLP fragment, after simplification of the AFLP fingerprint and separation of the fragments by migration. The second one is a combination of classical AFLP protocol and 454 pyrosequencing.

Does a shift in host plants trigger speciation in the Alpine leaf beetle Oreina speciosissima (Coleoptera, Chrysomelidae)?

BACKGROUND

Within the Coleoptera, the largest order in the animal kingdom, the exclusively herbivorous Chrysomelidae are recognized as one of the most species rich beetle families. The evolutionary processes that have fueled radiation into the more than thirty-five thousand currently recognized leaf beetle species remain partly unresolved. The prominent role of leaf beetles in the insect world, their omnipresence across all terrestrial biomes and their economic importance as common agricultural pest organisms make this family particularly interesting for studying the mechanisms that drive diversification. Here we specifically focus on two ecotypes of the alpine leaf beetle Oreina speciosissima (Scop.), which have been shown to exhibit morphological differences in male genitalia roughly corresponding to the subspecies Oreina speciosissima sensu stricto and Oreina speciosissima troglodytes. In general the two ecotypes segregate along an elevation gradient and by host plants: Oreina speciosissima sensu stricto colonizes high forb vegetation at low altitude and Oreina speciosissima troglodytes is found in stone run vegetation at higher elevations. Both host plants and leaf beetles have a patchy geographical distribution. Through use of gene sequencing and genome fingerprinting (AFLP) we analyzed the genetic structure and habitat use of Oreina speciosissima populations from the Swiss Alps to examine whether the two ecotypes have a genetic basis. By investigating a wide range of altitudes and focusing on the structuring effect of habitat types, we aim to provide answers regarding the factors that drive adaptive radiation in this phytophagous leaf beetle.

RESULTS

While little phylogenetic resolution was observed based on the sequencing of four DNA regions, the topology and clustering resulting from AFLP genotyping grouped specimens according to their habitat, mostly defined by plant associations. A few specimens with intermediate morphologies clustered with one of the two ecotypes or formed separate clusters consistent with habitat differences. These results were discussed in an ecological speciation framework.

CONCLUSIONS

The question of whether this case of ecological differentiation occurred in sympatry or allopatry remains open. Still, the observed pattern points towards ongoing divergence between the two ecotypes which is likely driven by a recent shift in host plant use.

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae)

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.