Fine‐scale spatial genetic structure of common and declining bumble bees across an agricultural landscape

Land‐use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land‐use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine‐scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri‐environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (F_IS = 0.01–0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri‐environment scheme conservation measures to facilitate fine‐scale gene flow by creating a more even distribution of suitable habitats across landscapes.     

Population genomics reveals a candidate gene involved in bumble bee pigmentation

Variation in bumble bee color patterns is well‐documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red‐banded to black‐banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome‐wide differentiation between red‐ and black‐banded forms. Here, we instead focus on the closely related black‐banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase‐like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.     

Local bumble bee decline linked to recovery of honey bees,drought effects on floral resources

Time series of abundances are critical for understanding how abiotic factors and species interactions affect population dynamics, but are rarely linked with experiments and also scarce for bee pollinators. This gap is important given concerns about declines in some bee species. I monitored honey bee (Apis mellifera) and bumble bee (Bombus spp.) foragers in coastal California from 1999, when feral A. mellifera populations were low due to Varroa destructor, until 2014. Apis mellifera increased substantially, except between 2006 and 2011, coinciding with declines in managed populations. Increases in A. mellifera strongly correlated with declines in Bombus and reduced diet overlap between them, suggesting resource competition consistent with past experimental results. Lower Bombus numbers also correlated with diminished floral resources. Declines in floral abundances were associated with drought and reduced spring rainfall. These results illustrate how competition with an introduced species may interact with climate to drive local decline of native pollinators.     

Distance models as a tool for modelling detection probability and density of native bumblebees

Effective monitoring of native bee populations requires accurate estimates of population size and relative abundance among habitats. Current bee survey methods, such as netting or pan trapping, may be adequate for a variety of study objectives but are limited by a failure to account for imperfect detection. Biases due to imperfect detection could result in inaccurate abundance estimates or erroneous insights about the response of bees to different environments. To gauge the potential biases of currently employed survey methods, we compared abundance estimates of bumblebees (Bombus spp.) derived from hierarchical distance sampling models (HDS) to bumblebee counts collected from fixed‐area net surveys (“net counts”) and fixed‐width transect counts (“transect counts”) at 47 early‐successional forest patches in Pennsylvania. Our HDS models indicated that detection probabilities of Bombus spp. were imperfect and varied with survey‐ and site‐covariates. Despite being conspicuous, Bombus spp. were not reliably detected beyond 5 m. Habitat associations of Bombus spp. density were similar across methods, but the strength of association with shrub cover differed between HDS and net counts. Additionally, net counts suggested sites with more grass hosted higher Bombus spp. densities whereas HDS suggested that grass cover was associated with higher detection probability but not Bombus spp. density. Density estimates generated from net counts and transect counts were 80%–89% lower than estimates generated from distance sampling. Our findings suggest that distance modelling provides a reliable method to assess Bombus spp. density and habitat associations, while accounting for imperfect detection caused by distance from observer, vegetation structure, and survey covariates. However, detection/non‐detection data collected via point‐counts, line‐transects and distance sampling for Bombus spp. are unlikely to yield species‐specific density estimates unless individuals can be identified by sight, without capture. Our results will be useful for informing the design of monitoring programs for Bombus spp. and other pollinators.     

Landscape genomic signatures indicate reduced gene flow and forest‐associated adaptive divergence in an endangered neotropical turtle

Human‐induced transformations of ecosystems usually result in fragmented populations subject to increased extinction risk. Fragmentation is also often associated with novel environmental heterogeneity, which in combination with restricted gene flow may increase the opportunity for local adaptation. To manage at‐risk populations in these landscapes, it is important to understand how gene flow is changing, and how populations respond to habitat loss. We conducted a landscape genomics analysis using Restriction‐site Associated DNA sequencing to investigate the evolutionary response of the critically endangered Dahl's Toad‐headed turtle (Mesoclemmys dahli) to severe habitat modification. The species has lost almost all of its natural habitat in the southwestern part of its range and about 70% in the northeast. Based on least cost path analysis across different resistance surfaces for 3,211 SNPs, we found that the landscape matrix is restricting gene flow, causing the fragmentation of the species into at least six populations. Genome scans and allele‐environment association analyses indicate that the population fragments in the deforested grasslands of the southwest are adaptively different from those in the more forested northeast. Populations in areas with no forest had low levels of adaptive genetic diversity and the fixation of ancestrally‐polymorphic SNPs, consistent with directional selection in this novel environment. Our results suggest that this forest‐stream specialist is adapting to pond‐grassland conditions, but it is also suffering from negative consequences of habitat loss, including genetic erosion, isolation, small effective population sizes, and inbreeding. We recommend gene flow restoration via genetic rescue to counteract these threats, and provide guidance for this strategy.     

Population structure and gene flow in the global pest,Helicoverpa armigera

Helicoverpa armigera is a major agricultural pest that is distributed across Europe, Asia, Africa and Australasia. This species is hypothesized to have spread to the Americas 1.5 million years ago, founding a population that is at present, a distinct species, Helicoverpa zea. In 2013, H. armigera was confirmed to have re‐entered South America via Brazil and subsequently spread. The source of the recent incursion is unknown and population structure in H. armigera is poorly resolved, but a basic understanding would highlight potential biosecurity failures and determine the recent evolutionary history of region‐specific lineages. Here, we integrate several end points derived from high‐throughput sequencing to assess gene flow in H. armigera and H. zea from populations across six continents. We first assemble mitochondrial genomes to demonstrate the phylogenetic relationship of H. armigera with other Heliothine species and the lack of distinction between populations. We subsequently use de novo genotyping‐by‐sequencing and whole‐genome sequences aligned to bacterial artificial chromosomes, to assess levels of admixture. Primarily, we find that Brazilian H. armigera are derived from diverse source populations, with strong signals of gene flow from European populations, as well as prevalent signals of Asian and African ancestry. We also demonstrate a potential field‐caught hybrid between H. armigera and H. zea, and are able to provide genomic support for the presence of the H. armigera conferta subspecies in Australasia. While structure among the bulk of populations remains unresolved, we present distinctions that are pertinent to future investigations as well as to the biosecurity threat posed by H. armigera.     

Revisiting comparisons of genetic diversity in stable and declining species: assessing genome‐wide polymorphism in North American bumble bees using RAD sequencing

Genetic variation is of key importance for a species’ evolutionary potential, and its estimation is a major component of conservation studies. New DNA sequencing technologies have enabled the analysis of large portions of the genome in nonmodel species, promising highly accurate estimates of such population genetic parameters. Restriction site‐associated DNA sequencing (RADseq) is used to analyse thousands of variants in the bumble bee species Bombus impatiens, which is common, and Bombus pensylvanicus, which is in decline. Previous microsatellite‐based analyses have shown that gene diversity is lower in the declining B. pensylvanicus than in B. impatiens. RADseq nucleotide diversities appear much more similar in the two species. Both species exhibit allele frequencies consistent with historical population expansions. Differences in diversity observed at microsatellites thus do not appear to have arisen from long‐term differences in population size and are either recent in origin or may result from mutational processes. Additional research is needed to explain these discrepancies and to investigate the best ways to integrate next‐generation sequencing data and more traditional molecular markers in studies of genetic diversity.     

The genomics of adaptation,divergence and speciation: a congealing theory

In this issue, Flaxman et al. ( 2014 ) report the results of sophisticated whole‐genome simulations of speciation with gene flow, enhancing our understanding of the process by building on previous single‐locus, multilocus and analytical works. Their findings provide us with new insights about how genomes can diverge and the importance of statistical and chromosomal linkage in facilitating reproductive isolation. The authors characterize the conditions under which, even with high gene flow and weak divergent selection, reproductive isolation between populations can occur due to the emergent stochastic process of genomewide congealing, where numerous statistically or physically linked loci of small effect allow selection to limit effective migration rates. The initial congealing event can occur within a broad range conditions, and once initiated, the self‐reinforcing process leads to rapid divergence and ultimately two reproductively isolated populations. Flaxman et al.'s ( 2014 ) work is a valuable contribution to our understanding of speciation with gene flow and in making a more predictive field of evolutionary genomics and speciation.     

Parasitoids as drivers of symbiont diversity in an insect host

Immune systems have repeatedly diversified in response to parasite diversity. Many animals have outsourced part of their immune defence to defensive symbionts, which should be affected by similar evolutionary pressures as the host’s own immune system. Protective symbionts provide efficient and specific protection and respond to changing selection pressure by parasites. Here we use the aphid Aphis fabae, its protective symbiont Hamiltonella defensa, and its parasitoid Lysiphlebus fabarum to test whether parasite diversity can maintain diversity in protective symbionts. We exposed aphid populations with the same initial symbiont composition to parasitoid populations that differed in their diversity. As expected, single parasitoid genotypes mostly favoured a single symbiont that was most protective against that particular parasitoid, while multiple symbionts persisted in aphids exposed to more diverse parasitoid populations, which in turn affected aphid population density and rates of parasitism. Parasite diversity may be crucial to maintaining symbiont diversity in nature.     

Climate structures genetic variation across a species' elevation range: a test of range limits hypotheses

Gene flow may influence the formation of species range limits, and yet little is known about the patterns of gene flow with respect to environmental gradients or proximity to range limits. With rapid environmental change, it is especially important to understand patterns of gene flow to inform conservation efforts. Here we investigate the species range of the selfing, annual plant, Mimulus laciniatus, in the California Sierra Nevada. We assessed genetic variation, gene flow, and population abundance across the entire elevation‐based climate range. Contrary to expectations, within‐population plant density increased towards both climate limits. Mean genetic diversity of edge populations was equivalent to central populations; however, all edge populations exhibited less genetic diversity than neighbouring interior populations. Genetic differentiation was fairly consistent and moderate among all populations, and no directional signals of contemporary gene flow were detected between central and peripheral elevations. Elevation‐driven gene flow (isolation by environment), but not isolation by distance, was found across the species range. These findings were the same towards high‐ and low‐elevation range limits and were inconsistent with two common centre‐edge hypotheses invoked for the formation of species range limits: (i) decreasing habitat quality and population size; (ii) swamping gene flow from large, central populations. This pattern demonstrates that climate, but not centre‐edge dynamics, is an important range‐wide factor structuring M. laciniatus populations. To our knowledge, this is the first empirical study to relate environmental patterns of gene flow to range limits hypotheses. Similar investigations across a wide variety of taxa and life histories are needed.     

Genetic divergence and isolation by thermal environment in geothermal populations of an aquatic invertebrate

Temperature is one of the most influential forces of natural selection impacting all biological levels. In the face of increasing global temperatures, studies over small geographic scales allowing investigations on the effects of gene flow are of great value for understanding thermal adaptation. Here, we investigated genetic population structure in the freshwater gastropod Radix balthica originating from contrasting thermal habitats in three areas of geothermal activity in Iceland. Snails from 32 sites were genotyped at 208 AFLP loci. Five AFLPs were identified as putatively under divergent selection in Lake Mývatn, a geothermal lake with an almost 20 °C difference in mean temperature across a distance of a few kilometres. In four of these loci, variation across all study populations was correlated with temperature. We found significant population structure in neutral markers both within and between the areas. Cluster analysis using neutral markers classified the sites mainly by geography, whereas analyses using markers under selection differentiated the sites based on temperature. Isolation by distance was stronger in the neutral than in the outlier loci. Pairwise differences based on outlier F_ST were significantly correlated with temperature at different spatial scales, even after correcting for geographic distance or neutral pairwise F_ST differences. In general, genetic variation decreased with increasing environmental temperature, possibly suggesting that natural selection had reduced the genetic diversity in the warm origin sites. Our results emphasize the influence of environmental temperature on the genetic structure of populations and suggest local thermal adaptation in these geothermal habitats.     

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.     

Genomic differentiation and patterns of gene flow between two long‐tailed tit species (Aegithalos)

Patterns of heterogeneous genomic differentiation have been well documented between closely related species, with some highly differentiated genomic regions (“genomic differentiation islands”) spread throughout the genome. Differential levels of gene flow are proposed to account for this pattern, as genomic differentiation islands are suggested to be resistant to gene flow. Recent studies have also suggested that genomic differentiation islands could be explained by linked selection acting on genomic regions with low recombination rates. Here, we investigate genomic differentiation and gene‐flow patterns for autosomes using RAD‐seq data between two closely related species of long‐tailed tits (Aegithalos bonvaloti and A. fuliginosus) in both allopatric and contact zone populations. The results confirm recent or ongoing gene flow between these two species. However, there is little evidence that the genomic regions that were found to be highly differentiated between the contact zone populations are resistant to gene flow, suggesting that differential levels of gene flow is not the cause of the heterogeneous genomic differentiation. Linked selection may be the cause of genomic differentiation islands between the allopatric populations with no or very limited gene flow, but this could not account for the heterogeneous genomic differentiation between the contact zone populations, which show evidence of recent or ongoing gene flow.     

Demographic and environmental drivers of metagenomic viral diversity in vampire bats

Viruses infect all forms of life and play critical roles as agents of disease, drivers of biochemical cycles and sources of genetic diversity for their hosts. Our understanding of viral diversity derives primarily from comparisons among host species, precluding insight into how intraspecific variation in host ecology affects viral communities or how predictable viral communities are across populations. Here we test spatial, demographic and environmental hypotheses explaining viral richness and community composition across populations of common vampire bats, which occur in diverse habitats of North, Central and South America. We demonstrate marked variation in viral communities that was not consistently predicted by a null model of declining community similarity with increasing spatial or genetic distances separating populations. We also find no evidence that larger bat colonies host greater viral diversity. Instead, viral diversity follows an elevational gradient, is enriched by juvenile‐biased age structure, and declines with local anthropogenic food resources as measured by livestock density. Our results establish the value of linking the modern influx of metagenomic sequence data with comparative ecology, reveal that snapshot views of viral diversity are unlikely to be representative at the species level, and affirm existing ecological theories that link host ecology not only to single pathogen dynamics but also to viral communities.     

Evolution of resistance to single and combined floral phytochemicals by a bumble bee parasite

Repeated exposure to inhibitory compounds can drive the evolution of resistance, which weakens chemical defence against antagonists. Floral phytochemicals in nectar and pollen have antimicrobial properties that can ameliorate infection in pollinators, but evolved resistance among parasites could diminish the medicinal efficacy of phytochemicals. However, multicompound blends, which occur in nectar and pollen, present simultaneous chemical challenges that may slow resistance evolution. We assessed evolution of resistance by the common bumble bee gut parasite Crithidia bombi to two floral phytochemicals, singly and combined, over 6 weeks (~100 generations) of chronic exposure. Resistance of C. bombi increased under single and combined phytochemical exposure, without any associated costs of reduced growth under phytochemical‐free conditions. After 6 weeks’ exposure, phytochemical concentrations that initially inhibited growth by > 50%, and exceeded concentrations in floral nectar, had minimal effects on evolved parasite lines. Unexpectedly, the phytochemical combination did not impede resistance evolution compared to single compounds. These results demonstrate that repeated phytochemical exposure, which could occur in homogeneous floral landscapes or with therapeutic phytochemical treatment of managed hives, can cause rapid evolution of resistance in pollinator parasites. We discuss possible explanations for submaximal phytochemical resistance in natural populations. Evolved resistance could diminish the antiparasitic value of phytochemical ingestion, weakening an important natural defence against infection.     

Genome divergence and diversification within a geographic mosaic of coevolution

Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome‐wide divergence was slight between ecotypes (F_ST = 0.011–0.035), we found evidence of relative genetic differentiation (as measured by F_ST) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome‐wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource‐based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000 years.     

Hybridization,polyploidy and clonality influence geographic patterns of diversity and salt tolerance in the model halophyte seashore paspalum (Paspalum vaginatum)

In plant species, variation in levels of clonality, ploidy and interspecific hybridization can interact to influence geographic patterns of genetic diversity. These factors commonly vary in plants that specialize on saline habitats (halophytes) and may play a role in how they adapt to salinity variation across their range. One such halophyte is the turfgrass and emerging genomic model system seashore paspalum (Paspalum vaginatum Swartz). To investigate how clonal propagation, ploidy variation, and interspecific hybridization vary across ecotypes and local salinity levels in wild P. vaginatum, we employed genotyping‐by‐sequencing, cpDNA sequencing and flow cytometry in 218 accessions representing > 170 wild collections from throughout the coastal southern United States plus USDA germplasm. We found that the two morphologically distinct ecotypes of P. vaginatum differ in their adaptive strategies. The fine‐textured ecotype is diploid and appears to reproduce in the wild both sexually and by clonal propagation; in contrast, the coarse‐textured ecotype consists largely of clonally‐propagating triploid and diploid genotypes. The coarse‐textured ecotype appears to be derived from hybridization between fine‐textured P. vaginatum and an unidentified Paspalum species. These clonally propagating hybrid genotypes are more broadly distributed than clonal fine‐textured genotypes and may represent a transition to a more generalist adaptive strategy. Additionally, the triploid genotypes vary in whether they carry one or two copies of the P. vaginatum subgenome, indicating multiple evolutionary origins. This variation in subgenome composition shows associations with local ocean salinity levels across the sampled populations and may play a role in local adaptation.     

The role of immigration and local adaptation on fine‐scale genotypic and phenotypic population divergence in a less mobile passerine

Dispersal and local patterns of adaptation play a major role on the ecological and evolutionary trajectory of natural populations. In this study, we employ a combination of genetic (25 microsatellite markers) and field‐based information (seven study years) to analyse the impact of immigration and local patterns of adaptation in two nearby (< 7 km) blue tit (Cyanistes caeruleus) populations. We used genetic assignment analyses to identify immigrant individuals and found that dispersal rate is female‐biased (72%). Data on lifetime reproductive success indicated that immigrant females produced fewer local recruits than their philopatric counterparts whereas immigrant males recruited more offspring than those that remained in their natal location. In spite of the considerably higher immigration rates of females, our results indicate that, in absolute terms, their demographic and genetic impact in the receiving populations is lower than that in immigrant males. Immigrants often brought novel alleles into the studied populations and a high proportion of them were transmitted to their recruits, indicating that the genetic impact of immigrants is not ephemeral. Although only a few kilometres apart, the two study populations were genetically differentiated and showed strong divergence in different phenotypic and life‐history traits. An almost absent inter‐population dispersal, together with the fact that both populations receive immigrants from different source populations, is probably the main cause of the observed pattern of genetic differentiation. However, phenotypic differentiation (P_ST) for all the studied traits greatly exceeded neutral genetic differentiation (F_ST), indicating that divergent natural selection is the prevailing factor determining the evolutionary trajectory of these populations. Our study highlights the importance of integrating individual‐ and population‐based approaches to obtain a comprehensive view about the role of dispersal and natural selection on structuring the genotypic and phenotypic characteristics of natural populations.     

Climate and geochemistry as drivers of eucalypt diversification in Australia

Eucalypts cover most of Australia. Here, we investigate the relative contribution of climate and geochemistry to the distribution and diversity of eucalypts. Using geostatistics, we estimate major element concentrations, pH, and electrical conductivity at sites where eucalypts have been recorded. We compare the median predicted geochemistry and reported substrate for individual species that appear associated with extreme conditions; this provides a partial evaluation of the predictions. We generate a site‐by‐species matrix by aggregating observations to the centroids of 100‐km‐wide grid cells, calculate diversity indices, and use numerical ecology methods (ordination, variation partitioning) to investigate the ecology of eucalypts and their response to climatic and geochemical gradients. We find that β‐diversity coincides with variations in climatic and geochemical patterns. Climate and geochemistry together account for less than half of the variation in eucalypt species assemblages across Australia but for greater than 80% in areas of high species richness. Climate is more important than geochemistry in explaining eucalypts species distribution and change in assemblages across Australia as a whole but there are correlations between the two sets of environmental variables. Many individual eucalypt species and entire taxonomic sections (Aromatica, Longistylus of subgenus Eucalyptus, Dumaria, and Liberivalvae of subgenus Symphyomyrtus) have distributions affected strongly by geochemistry. We conclude that eucalypt diversity is driven by steep geochemical gradients that have arisen as climate patterns have fluctuated over Australia over the Cenozoic, generally aridifying since the Miocene. The diversification of eucalypts across Australia is thus an excellent example of co‐evolution of landscapes and biota in space and time and challenges accepted notions of macroecology.