Subdivision and haplotype structure in natural populations of Arabidopsis lyrata

High levels of inbreeding are expected to cause a strong reduction in levels of genetic variability, effective recombination rates and in adaptation compared with related outcrossing populations. We examined patterns of DNA polymorphism at five nuclear loci and one chloroplast locus within and between four populations of the outcrossing plant Arabidopsis lyrata, a close relative of the highly self-fertilizing model species A. thaliana. The observed patterns are compared with species-wide polymorphism at orthologous loci, as well as within- and between-population patterns at other studied loci in A. thaliana. In addition to evidence for much higher average within-population diversity, species-wide levels of silent polymorphism are generally higher in A. lyrata than in A. thaliana, unlike the results from a previous study of the ADH locus. However, polymorphism is also low in the North American A. lyrata subspecies lyrata compared with the European subspecies petraea, suggesting either a population bottleneck in North American populations or recent admixture involving diverged European populations. Differentiation between the two subspecies is strong, although there are few fixed differences, suggesting that their isolation is recent. Estimates of intralocus recombination rates and analysis of haplotype structure in European A. lyrata populations indicate lower recombination than predicted based on the variability together with physical recombination rates estimated from A. thaliana. This may be due to strong population subdivision, or to recent departures from demographic equilibrium such as a bottleneck or population admixture. Alternatively, there may be consistently lower recombination rates in the outcrossing species. In contrast, estimates of recombination rates from species-wide samples of A. thaliana are close to the values expected assuming a high rate of self-fertilization. Complex population histories in both A. thaliana and A. lyrata complicate theoretical predictions and empirical tests of the effects of inbreeding on polymorphism and molecular evolution.     

Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.     

GENETIC LEAKAGE AFTER ADAPTIVE AND NONADAPTIVE DIVERGENCE IN THE ENSATINA ESCHSCHOLTZII RING SPECIES

The salamander Ensatina eschscholtzii is an example of a ring species in which extant intermediate stages of terminal forms have a nearly continuous range, offering replicated interactions at several stages of divergence. We employ a greatly expanded allozyme database and individual-based analyses to separate the effects of divergence time and gene flow to evaluate how gradual divergence of populations around the ring contributes to the development of reproductive isolation. Despite the high degree of genetic ( D ≤ 0.39) and ecomorphological divergence observed in secondary contacts around the ring, reproductive isolation or rare hybridization is observed only at the terminus of the ring. Instead, in the secondary contacts sampled around the ring, hybrids are common and reproductively successful, enabling genetic leakage between parental genomes and the potential for genetic merger. Nevertheless, genetic admixture is geographically broad (<100 km) only in contacts between ecomorphologically similar populations (within subspecies). When divergence is accompanied by alternative patterns of adaptive divergence (between subspecies), zones of intergradation are narrower and affect populations only locally (>8 km). Diversification and consequent genetic interactions in Ensatina reveal a continuum between populations, ecological races, and species, where polytypic traits and high genetic differentiation are maintained without reproductive isolation.     

DNA barcoding will often fail to discover new animal species over broad parameter space

With increasing force, genetic divergence of mitochondrial DNA (mtDNA) is being argued as the primary tool for discovery of animal species. Two thresholds of single-gene divergence have been proposed: reciprocal monophyly, and 10 times greater genetic divergence between than within species (the "10x rule"). To explore quantitatively the utility of each approach, we couple neutral coalescent theory and the classical Bateson-Dobzhansky-Muller (BDM) model of speciation. The joint stochastic dynamics of these two processes demonstrate that both thresholds fail to "discover" many reproductively isolated lineages under a single incompatibility BDM model, especially when BDM loci have been subject to divergent selection. Only when populations have been isolated for > 4 million generations did these thresholds achieve error rates of < 10% under our model that incorporates variable population sizes. The high error rate evident in simulations is corroborated with six empirical data sets. These properties suggest that single-gene, high-throughput approaches to discovering new animal species will bias large-scale biodiversity surveys, particularly toward missing reproductively isolated lineages that have emerged by divergent selection or other mechanisms that accelerate reproductive isolation. Because single-gene thresholds for species discovery can result in substantial error at recent divergence times, they will misrepresent the correspondence between recently isolated populations and reproductively isolated lineages (= species).     

A Widespread Chromosomal Inversion Polymorphism Contributes to a Major Life-History Transition,Local Adaptation,and Reproductive Isolation

The role of chromosomal inversions in adaptation and speciation is controversial. Historically, inversions were thought to contribute to these processes either by directly causing hybrid sterility or by facilitating the maintenance of co-adapted gene complexes. Because inversions suppress recombination when heterozygous, a recently proposed local adaptation mechanism predicts that they will spread if they capture alleles at multiple loci involved in divergent adaptation to contrasting environments. Many empirical studies have found inversion polymorphisms linked to putatively adaptive phenotypes or distributed along environmental clines. However, direct involvement of an inversion in local adaptation and consequent ecological reproductive isolation has not to our knowledge been demonstrated in nature. In this study, we discovered that a chromosomal inversion polymorphism is geographically widespread, and we test the extent to which it contributes to adaptation and reproductive isolation under natural field conditions. Replicated crosses between the prezygotically reproductively isolated annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, revealed that alternative chromosomal inversion arrangements are associated with life-history divergence over thousands of kilometers across North America. The inversion polymorphism affected adaptive flowering time divergence and other morphological traits in all replicated crosses between four pairs of annual and perennial populations. To determine if the inversion contributes to adaptation and reproductive isolation in natural populations, we conducted a novel reciprocal transplant experiment involving outbred lines, where alternative arrangements of the inversion were reciprocally introgressed into the genetic backgrounds of each ecotype. Our results demonstrate for the first time in nature the contribution of an inversion to adaptation, an annual/perennial life-history shift, and multiple reproductive isolating barriers. These results are consistent with the local adaptation mechanism being responsible for the distribution of the two inversion arrangements across the geographic range of M. guttatus and that locally adaptive inversion effects contribute directly to reproductive isolation. Such a mechanism may be partially responsible for the observation that closely related species often differ by multiple chromosomal rearrangements.     

Costs and benefits of admixture between foreign genotypes and local populations in the field

Admixture is the hybridization between populations within one species. It can increase plant fitness and population viability by alleviating inbreeding depression and increasing genetic diversity. However, populations are often adapted to their local environments and admixture with distant populations could break down local adaptation by diluting the locally adapted genomes. Thus, admixed genotypes might be selected against and be outcompeted by locally adapted genotypes in the local environments. To investigate the costs and benefits of admixture, we compared the performance of admixed and within‐population F1 and F2 generations of the European plant Lythrum salicaria in a reciprocal transplant experiment at three European field sites over a 2‐year period. Despite strong differences between site and plant populations for most of the measured traits, including herbivory, we found limited evidence for local adaptation. The effects of admixture depended on experimental site and plant population, and were positive for some traits. Plant growth and fruit production of some populations increased in admixed offspring and this was strongest with larger parental distances. These effects were only detected in two of our three sites. Our results show that, in the absence of local adaptation, admixture may boost plant performance, and that this is particularly apparent in stressful environments. We suggest that admixture between foreign and local genotypes can potentially be considered in nature conservation to restore populations and/or increase population viability, especially in small inbred or maladapted populations.     

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

Nuclear-cytoplasmic interactions reduce male fertility in hybrids of Arabidopsis lyrata subspecies

We examined the level of postzygotic reproductive isolation in F(1) and F(2) hybrids of reciprocal crosses between the Arabidopsis lyrata subspecies lyrata (North American) and petraea (European). Our main results are: first, the percentage of fertile pollen was significantly reduced in the F(1) and F(2) compared to the parental populations. Second, mean pollen fertility differed markedly between reciprocal crosses: 84% in the F(2) with ssp. lyrata cytoplasm and 61% in the F(2) with ssp. petraea cytoplasm. Third, 17% of the F(2) with ssp. petraea cytoplasm showed male sterility (produced less than 30 pollen grains in our subsample). The hybrids were female fertile. We used QTL mapping to find the genomic regions that determine pollen fertility and that restore cytoplasmic male sterility (CMS). In the F(2) with ssp. lyrata cytoplasm, an epistatic pair of QTLs was detected. In the reciprocal F(2) progeny, four QTLs demonstrated within-population polymorphism for hybrid male sterility. In addition, in the F(2) with ssp. petraea cytoplasm, there was a strong male fertility restorer locus on chromosome 2 where a cluster of CMS restorer gene-related PPR genes have been found in A. lyrata. Our results underline the importance of cytonuclear interactions in understanding genetics of the early stages of speciation.     

Local adaptation across a climatic gradient despite small effective population size in the rare sapphire rockcress

When assigning conservation priorities in endangered species, two common management strategies seek to protect remnant populations that (i) are the most genetically divergent or (ii) possess the highest diversity at neutral genetic markers. These two approaches assume that variation in molecular markers reflects variation in ecologically important traits and ignore the possibility of local adaptation among populations that show little divergence or variation at marker loci. Using common garden experiments, we demonstrate that populations of the rare endemic plant Arabis fecunda are physiologically adapted to the local microclimate. Local adaptation occurs despite (i) the absence of divergence at almost all marker loci and (ii) very small effective population sizes, as evidenced by extremely low levels of allozyme and DNA sequence polymorphism. Our results provide empirical evidence that setting conservation priorities based exclusively on molecular marker diversity may lead to the loss of locally adapted populations.     

Contrasting patterns of nucleotide polymorphism at the alcohol dehydrogenase locus in the outcrossing Arabidopsis lyrata and the selfing Arabidopsis thaliana

Nucleotide variation at the alcohol dehydrogenase locus (Adh) was studied in the outcrossing Arabidopsis lyrata, a close relative of the selfing Arabidopsis thaliana. Overall, estimated nucleotide diversity in the North American ssp. lyrata and two European ssp. petraea populations was 0.0038, lower than the corresponding specieswide estimate for A. thaliana at the same set of nucleotide sites. The distribution of segregating sites across the gene differed between the two species. Estimated sequence diversity within an A. lyrata population with a large sample size (0.0023) was much higher than has previously been observed for A. thaliana. This North American population has an excess of sites at intermediate frequencies compared with neutral expectation (Tajima's D = 2.3, P < 0.005), suggestive of linked balancing selection or a recent population bottleneck. In contrast, an excess of rare polymorphisms has been found in A. thaliana. Polymorphism within A. lyrata and divergence from A. thaliana appear to be correlated across the Adh gene sequence. The geographic distribution of polymorphism was quite different from that of A. thaliana, for which earlier studies of several genes found low within-population nucleotide site polymorphism and no overall continental differentiation of variation despite large differences in site frequencies between local populations. Differences between the outcrossing A. lyrata and the selfing A. thaliana reflect the impact of differences in mating system and the influence of bottlenecks in A. thaliana during rapid colonization on DNA sequence polymorphism. The influence of additional variability-reducing mechanisms, such as background selection or hitchhiking, may not be discernible.     

MULTIPLE MODES OF SPECIATION INVOLVED IN THE PARALLEL EVOLUTION OF SYMPATRIC MORPHOTYPES OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS,SALMONIDAE)

We performed a phylogenetic analysis of mtDNA variation among seven sympatric pairs of dwarf and normal morphotypes of whitefish from northern Québec and the St. John River drainage to address three questions relevant to understanding their radiation. Are all sympatric pairs reproductively isolated? Do phylogenetic analyses confirm that sympatric whitefish morphotypes found in eastern North America represent the outcome of polyphyletic evolutionary events? If so, did all sympatric pairs from the St. John River drainage originate from the same scenario of allopatric divergence and secondary contact? The hypothesis of genetic differentiation was supported for all sympatric pairs from the St. John River drainage, whereas lack of mtDNA diversity precluded any test of reproductive isolation for northern Québec populations. Patterns of mtDNA variation confirmed that dwarf and normal morphotypes evolved in parallel among independent, yet closely related, lineages, thus providing indirect evidence for the role of natural selection in promoting phenotypic radiation in whitefish. Patterns of mtDNA diversity among sympatric pairs of the St. John River indicated a complex picture of whitefish evolution that implied sympatric divergence and multiple allopatric divergence/secondary contact events on a small geographic scale. These results suggests that ecological opportunities, namely trophic niche availability, may promote population divergence in whitefish.     

Coalescent patterns for chromosomal inversions in divergent populations

Chromosomal inversions allow genetic divergence of locally adapted populations by reducing recombination between chromosomes with different arrangements. Divergence between populations (or hybridization between species) is expected to leave signatures in the neutral genetic diversity of the inverted region. Quantitative expectations for these patterns, however, have not been obtained. Here, we develop coalescent models of neutral sites linked to an inversion polymorphism in two locally adapted populations. We consider two scenarios of local adaptation: selection on the inversion breakpoints and selection on alleles inside the inversion. We find that ancient inversion polymorphisms cause genetic diversity to depart dramatically from neutral expectations. Other situations, however, lead to patterns that may be difficult to detect; important determinants are the age of the inversion and the rate of gene flux between arrangements. We also study inversions under genetic drift, finding that they produce patterns similar to locally adapted inversions of intermediate age. Our results are consistent with empirical observations, and provide the foundation for quantitative analyses of the roles that inversions have played in speciation.     

Reproductive capacity of varroa destructor in four different honey bee subspecies

Varroa tolerance as a consequence of host immunity may contribute substantially to reduce worldwide colony declines. Therefore, special breeding programs were established and varroa surviving populations investigated to understand mechanisms behind this adaptation. The aim of this study was to investigate the reproductive capacity in the three most common subspecies of the European honey bee (Carnica, Mellifera, Ligustica) and the F2 generation of a varroa surviving population, to identify if managed host populations possibly have adapted over time already. Both, singly infested drone and worker brood were assessed to determine fertility and fecundity of varroa foundresses in their respective group. We found neither parameter to be significantly different within the four subspecies, demonstrating that no adaptations have occurred in terms of the reproductive success of Varroa destructor. In all groups mother mites reproduce equally successful and are potentially able to cause detrimental damage to their host when not being treated sufficiently. The data further suggests that a population once varroa tolerant does not necessarily inherit this trait to following generations after the F1, which could be of particular interest when selecting populations for resistance breeding. Reasons and consequences are discussed.     

Restricted gene flow and local adaptation highlight the vulnerability of high‐latitude reefs to rapid environmental change

Global climate change poses a serious threat to the future health of coral reef ecosystems. This calls for management strategies that are focused on maximizing the evolutionary potential of coral reefs. Fundamental to this is an accurate understanding of the spatial genetic structure in dominant reef‐building coral species. In this study, we apply a genotyping‐by‐sequencing approach to investigate genome‐wide patterns of genetic diversity, gene flow, and local adaptation in a reef‐building coral, Pocillopora damicornis, across 10 degrees of latitude and a transition from temperate to tropical waters. We identified strong patterns of differentiation and reduced genetic diversity in high‐latitude populations. In addition, genome‐wide scans for selection identified a number of outlier loci putatively under directional selection with homology to proteins previously known to be involved in heat tolerance in corals and associated with processes such as photoprotection, protein degradation, and immunity. This study provides genomic evidence for both restricted gene flow and local adaptation in a widely distributed coral species, and highlights the potential vulnerability of leading‐edge populations to rapid environmental change as they are locally adapted, reproductively isolated, and have reduced levels of genetic diversity.     

Detecting adaptive trait loci in nonmodel systems: divergence or admixture mapping?

Mapping adaptive trait loci (ATL) underlying ecological divergence is an essential step towards understanding the processes that generate phenotypic diversity. Technological advances have made it possible to sequence exomes in nonmodel systems, providing an efficient means of analysing functional genetic variants. Divergence scans of genetic markers for outlier loci, or ‘divergence mapping’, have been used to map locally adapted genes, but this approach is likely to be underpowered when background divergence is elevated. Genotype–phenotype association tests in admixed populations, or ‘admixture mapping’, may provide a useful approach for mapping locally adapted loci when neutral divergence is high. To determine the power and limits of divergence mapping, we simulated exomes containing a single ATL across two parental populations of varying neutral divergence, estimated divergence and quantified the power to identify the ATL. We found that divergence mapping had very high power when background F_ST is <0.2, but decreased dramatically above this level. To evaluate the utility of admixture mapping, we simulated exomes from admixed populations, then simulated phenotypes, conducted genotype–phenotype association tests and found that even two generations of random mating after admixture could provide high mapping power in scenarios where pure divergence mapping was ineffective (F_ST = 0.35). Moreover, admixture mapping had high power across all levels of divergence after 20 generations since admixture. Together with high‐throughput exome sequencing, admixture mapping could be used to map ATL in systems such as Heliconius butterflies or Gryllus crickets when experimental design and analytical approach are chosen accordingly.     

MOLECULAR AND REPRODUCTIVE CHARACTERIZATION OF SIBLING SPECIES IN THE EUROPEAN EARWIG (FORFICULA AURICULARIA)

The European earwig (Forficula auricularia) was formerly thought to present a mosaic of populations differing in their reproductive biology. We show that it is comprised of two as yet unrecognized sibling species. The molecular divergence between the two species, for a 627-bp amplified fragment overlapping the COI and COII mitochondrial loci, is six times larger than intraspecific variation. A species with two clutches a year lives predominantly in lowland and oceanic European habitats. A species with one clutch a year—except in the Mediterranean area where it has two clutches—lives predominantly in highland and continental European habitats. They both invaded North America during the 20th century, respectively, from the west and the east coasts, with no apparent mixing of their populations. The two species can occur in sympatry in Europe and are reproductively isolated by nearly complete failure to produce F₁ hybrids.     

Thermal adaptation and ecological speciation

Ecological speciation is defined as the emergence of reproductive isolation as a direct or indirect consequence of divergent ecological adaptation. Several empirical examples of ecological speciation have been reported in the literature which very often involve adaptation to biotic resources. In this review, we investigate whether adaptation to different thermal habitats could also promote speciation and try to assess the importance of such processes in nature. Our survey of the literature identified 16 animal and plant systems where divergent thermal adaptation may underlie (partial) reproductive isolation between populations or may allow the stable coexistence of sibling taxa. In many of the systems, the differentially adapted populations have a parapatric distribution along an environmental gradient. Isolation often involves extrinsic selection against locally maladapted parental or hybrid genotypes, and additional pre- or postzygotic barriers may be important. Together, the identified examples strongly suggest that divergent selection between thermal environments is often strong enough to maintain a bimodal genotype distribution upon secondary contact. What is less clear from the available data is whether it can also be strong enough to allow ecological speciation in the face of gene flow through reinforcement-like processes. It is possible that intrinsic features of thermal gradients or the genetic basis of thermal adaptation make such reinforcement-like processes unlikely but it is equally possible that pertinent systems are understudied. Overall, our literature survey highlights (once again) the dearth of studies that investigate similar incipient species along the continuum from initial divergence to full reproductive isolation and studies that investigate all possible reproductive barriers in a given system.     

Genetic hitchhiking associated with life history divergence and colonization of North America in the European corn borer moth

A primary goal for evolutionary biology is to reveal the genetic basis for adaptive evolution and reproductive isolation. Using Z and E pheromone strains the European corn borer (ECB) moth, I address this problem through multilocus analyses of DNA polymorphism. I find that the locus Triose phosphate isomerase (Tpi) is a statistically significant outlier in coalescent simulations of demographic histories of population divergence, including strict allopatric isolation, restricted migration, secondary contact, and population growth or decline. This result corroborates a previous QTL study that identified the Tpi chromosomal region as a repository for gene(s) contributing to divergence in life history. Patterns of nucleotide polymorphism at Tpi suggest a recent selective sweep and genetic hitchhiking associated with colonization of North America from Europe ~200 generations ago. These results indicate that gene genealogies initially diverge during speciation because of selective sweeps, but differential introgression may play a role in the maintenance of differentiation for sympatric populations.     

Estimating the impact of divergent mating phenology between residents and migrants on the potential for gene flow

Gene flow between populations can allow the spread of beneficial alleles and genetic diversity between populations, with importance to conservation, invasion biology, and agriculture. Levels of gene flow between populations vary not only with distance, but also with divergence in reproductive phenology. Since phenology is often locally adapted, arriving migrants may be reproductively out of synch with residents, which can depress realized gene flow. In flowering plants, the potential impact of phenological divergence on hybridization between populations can be predicted from overlap in flowering schedules—the daily count of flowers capable of pollen exchange—between a resident and migrant population. The accuracy of this prospective hybridization estimate, based on parental phenotypes, rests upon the assumptions of unbiased pollen transfer between resident and migrant active flowers. We tested the impact of phenological divergence on resident–migrant mating frequencies in experiments that mimicked a single large gene flow event. We first prospectively estimated mating frequencies two lines of Brassica rapaselected or early and late flowering. We then estimated realized mating frequencies retrospectively through progeny testing. The two estimates strongly agreed in a greenhouse experiment, where procedures ensured saturating, unbiased pollination. Under natural pollination in the field, the rate of resident–migrant mating, was lower than estimated by phenological divergence alone, although prospective and retrospective estimates were correlated. In both experiments, differences between residents and migrants in flowering schedule shape led to asymmetric hybridization. Results suggest that a prospective estimate of hybridization based on mating schedules can be a useful, although imperfect, tool for evaluating potential gene flow. They also illustrate the impact of mating phenology on the magnitude and symmetry of reproductive isolation.     

Genomic analysis of differentiation between soil types reveals candidate genes for local adaptation in Arabidopsis lyrata

Serpentine soil, which is naturally high in heavy metal content and has low calcium to magnesium ratios, comprises a difficult environment for most plants. An impressive number of species are endemic to serpentine, and a wide range of non-endemic plant taxa have been shown to be locally adapted to these soils. Locating genomic polymorphisms which are differentiated between serpentine and non-serpentine populations would provide candidate loci for serpentine adaptation. We have used the Arabidopsis thaliana tiling array, which has 2.85 million probes throughout the genome, to measure genetic differentiation between populations of Arabidopsis lyrata growing on granitic soils and those growing on serpentinic soils. The significant overrepresentation of genes involved in ion transport and other functions provides a starting point for investigating the molecular basis of adaptation to soil ion content, water retention, and other ecologically and economically important variables. One gene in particular, calcium-exchanger 7, appears to be an excellent candidate gene for adaptation to low CaratioMg ratio in A. lyrata.