Urban land use limits regional bumble bee gene flow

Potential declines in native pollinator communities and increased reliance on pollinator‐dependent crops have raised concerns about native pollinator conservation and dispersal across human‐altered landscapes. Bumble bees are one of the most effective native pollinators and are often the first to be extirpated in human‐altered habitats, yet little is known about how bumble bees move across fine spatial scales and what landscapes promote or limit their gene flow. In this study, we examine regional genetic differentiation and fine‐scale relatedness patterns of the yellow‐faced bumble bee, Bombus vosnesenskii, to investigate how current and historic habitat composition impact gene flow. We conducted our study across a landscape mosaic of natural, agricultural and urban/suburban habitats, and we show that B. vosnesenskii exhibits low but significant levels of differentiation across the study system (F_ST = 0.019, D_est = 0.049). Most importantly, we reveal significant relationships between pairwise F_ST and resistance models created from contemporary land use maps. Specifically, B. vosnesenskii gene flow is most limited by commercial, industrial and transportation‐related impervious cover. Finally, our fine‐scale analysis reveals significant but declining relatedness between individuals at the 1–9 km spatial scale, most likely due to local queen dispersal. Overall, our results indicate that B. vosnesenskii exhibits considerable local dispersal and that regional gene flow is significantly limited by impervious cover associated with urbanization.     

Contemporary human‐altered landscapes and oceanic barriers reduce bumble bee gene flow

Much of the world's terrestrial landscapes are being altered by humans in the form of agriculture, urbanization and pastoral systems, with major implications for biodiversity. Bumble bees are one of the most effective pollinators in both natural and cultivated landscapes, but are often the first to be extirpated in human‐altered habitats. Yet, little is known about the role of natural and human‐altered habitats in promoting or limiting bumble bee gene flow. In this study, I closely examine the genetic structure of the yellow‐faced bumble bee, Bombus vosnesenskii, across the southwestern US coast and find strong evidence that natural oceanic barriers, as well as contemporary human‐altered habitats, limit bee gene flow. Heterozygosity and allelic richness were lower in island populations, while private allelic richness was higher in island populations compared to mainland populations. Genetic differentiation, measured for three indices across the 1000 km study region, was significantly greater than the null expectation (F_ST = 0.041, F’_ST = 0.044 and D_est = 0.155) and correlated with geographic distance. Furthermore, genetic differentiation patterns were most strongly correlated with contemporary (2011) not past (2006, 2001) resistance maps calibrated for high dispersal limitation over oceans, impervious habitat and croplands. Despite the incorporation of dramatic elevation gradients, the analyses reveal that oceans and contemporary human land use, not mountains, are the primary dispersal barriers for B. vosnesenskii gene flow. These findings reinforce the importance of maintaining corridors of suitable habitat across the distribution range of native pollinators to promote their persistence and safeguard their ability to provide essential pollination services.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

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.     

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.     

Interannual bumble bee abundance is driven by indirect climate effects on floral resource phenology

Climate change can influence consumer populations both directly, by affecting survival and reproduction, and indirectly, by altering resources. However, little is known about the relative importance of direct and indirect effects, particularly for species important to ecosystem functioning, like pollinators. We used structural equation modelling to test the importance of direct and indirect (via floral resources) climate effects on the interannual abundance of three subalpine bumble bee species. In addition, we used long‐term data to examine how climate and floral resources have changed over time. Over 8 years, bee abundances were driven primarily by the indirect effects of climate on the temporal distribution of floral resources. Over 43 years, aspects of floral phenology changed in ways that indicate species‐specific effects on bees. Our study suggests that climate‐driven alterations in floral resource phenology can play a critical role in governing bee population responses to global change.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Reduced abundance and earlier collection of bumble bee workers under intensive cultivation of a mass‐flowering prairie crop

One of the most commonly seeded crops in Canada is canola, a cultivar of oilseed rape (Brassica napus). As a mass‐flowering crop grown intensively throughout the Canadian Prairies, canola has the potential to influence pollinator success across tens of thousands of square kilometers of cropland. Bumble bees (Bombus sp.) are efficient pollinators of many types of native and crop plants. We measured the influence of this mass‐flowering crop on the abundance and phenology of bumble bees, and on another species of social bee (a sweat bee; Halictus rubicundus), by continuously deploying traps at different levels of canola cultivation intensity, spanning the start and end of canola bloom. Queen bumble bees were more abundant in areas with more canola cover, indicating that this crop is attractive to queens. However, bumble bee workers were significantly fewer in these locations later in the season, suggesting reduced colony success. The median collection dates of workers of three bumble bee species were earlier near canola fields, suggesting a dynamic response of colonies to the increased floral resources. Different species experienced this shift to different extents. The sweat bee was not affected by canola cultivation intensity. Our findings suggest that mass‐flowering crops such as canola are attractive to bumble bee queens and therefore may lead to higher rates of colony establishment, but also that colonies established near this crop may be less successful. We propose that the effect on bumble bees can be mitigated by spacing the crop more evenly with respect to alternate floral resources.     

Population genomics fits the bill: genetics of adaptive beak variation in Darwin's finches

Darwin's finches are an iconic case of adaptive radiation. The size and shape of their beaks are key adaptive traits related to trophic niche that vary among species and evolve rapidly when the food supply changes. Building on recent studies, a paper in this issue of Molecular Ecology (Chaves et al. 2016 ) investigates the genomic basis of beak size variation in sympatric populations of three species of ground finch (Geospiza) by performing a Genome‐wide association study using RAD‐seq data. The authors find that variation in a small number of markers can explain a substantial proportion of variation in beak size. Some of these markers are in genomic regions that have previously been implicated in beak size variation in Darwin's finches, whereas other markers have not, suggesting both conservation and divergence in the genetic basis of morphological evolution. Overall, the study confirms that loci of large effect are involved in beak size variation, which helps to explain the high heritability and rapid response to selection of this trait. The independent identification of regions containing HMGA2 and DLK1 loci in a GWAS makes them prime targets for functional studies. The study also shows that under the right conditions, RAD‐seq can be a viable alternative to genome sequencing for GWAS in wild vertebrate populations.     

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.     

Mito‐nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern “Rocky Mountain''’ color form with ferruginous mid‐abdominal segments (B. m. melanopygus) and a southern “Pacific''’ form with black mid‐abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color‐controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito‐nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.     

Physiological and regulatory underpinnings of geographic variation in reptilian cold tolerance across a latitudinal cline

Understanding the mechanisms that produce variation in thermal performance is a key component to investigating climatic effects on evolution and adaptation. However, disentangling the effects of local adaptation and phenotypic plasticity in shaping patterns of geographic variation in natural populations can prove challenging. Additionally, the physiological mechanisms that cause organismal dysfunction at extreme temperatures are still largely under debate. Using the green anole, Anolis carolinensis, we integrate measures of cold tolerance (CT_min), standard metabolic rate, heart size, blood lactate concentration and RNAseq data from liver tissue to investigate geographic variation in cold tolerance and its underlying mechanisms along a latitudinal cline. We found significant effects of thermal acclimation and latitude of origin on variation in cold tolerance. Increased cold tolerance correlates with decreased rates of oxygen consumption and blood lactate concentration (a proxy for oxygen limitation), suggesting elevated performance is associated with improved oxygen economy during cold exposure. Consistent with these results, co‐expression modules associated with blood lactate concentration are enriched for functions associated with blood circulation, coagulation and clotting. Expression of these modules correlates with thermal acclimation and latitude of origin. Our findings support the oxygen and capacity‐limited thermal tolerance hypothesis as a potential contributor to variation in reptilian cold tolerance. Moreover, differences in gene expression suggest regulation of the blood coagulation cascade may play an important role in reptilian cold tolerance and may be the target of natural selection in populations inhabiting colder environments.     

Comparative genomics analysis in Prunoideae to identify biologically relevant polymorphisms

Prunus is an economically important genus with a wide range of physiological and biological variability. Using the peach genome as a reference, sequencing reads from four almond accessions and one sweet cherry cultivar were used for comparative analysis of these three Prunus species. Reference mapping enabled the identification of many biological relevant polymorphisms within the individuals. Examining the depth of the polymorphisms and the overall scaffold coverage, we identified many potentially interesting regions including hundreds of small scaffolds with no coverage from any individual. Non‐sense mutations account for about 70 000 of the 13 million identified single nucleotide polymorphisms (SNPs). Blast2GO analyses on these non‐sense SNPs revealed several interesting results. First, non‐sense SNPs were not evenly distributed across all gene ontology terms. Specifically, in comparison with peach, sweet cherry is found to have non‐sense SNPs in two 1‐aminocyclopropane‐1‐carboxylate synthase (ACS) genes and two 1‐aminocyclopropane‐1‐carboxylate oxidase (ACO) genes. These polymorphisms may be at the root of the nonclimacteric ripening of sweet cherry. A set of candidate genes associated with bitterness in almond were identified by comparing sweet and bitter almond sequences. To the best of our knowledge, this is the first report in plants of non‐sense SNP abundance in a genus being linked to specific GO terms.     

Molecular diversity and landscape genomics of the crop wild relative Triticum urartu across the Fertile Crescent

Modern plant breeding can benefit from the allelic variation that exists in natural populations of crop wild relatives that evolved under natural selection in varying pedoclimatic conditions. In this study, next‐generation sequencing was used to generate 1.3 million genome‐wide single nucleotide polymorphisms (SNPs) on ex situ collections of Triticum urartu L., the wild donor of the A^u subgenome of modern wheat. A set of 75 511 high‐quality SNPs were retained to describe 298 T. urartu accessions collected throughout the Fertile Crescent. Triticum urartu showed a complex pattern of genetic diversity, with two main genetic groups distributed sequentially from west to east. The incorporation of geographical information on sampling points showed that genetic diversity was correlated to the geographical distance (R² = 0.19) separating samples from Jordan and Lebanon, from Syria and southern Turkey, and from eastern Turkey, Iran and Iraq. The wild emmer genome was used to derive the physical positions of SNPs on the seven chromosomes of the A^u subgenome, allowing us to describe a relatively slow decay of linkage disequilibrium in the collection. Outlier loci were described on the basis of the geographic distribution of the T. urartu accessions, identifying a hotspot of directional selection on chromosome 4A. Bioclimatic variation was derived from grid data and related to allelic variation using a genome‐wide association approach, identifying several marker–environment associations (MEAs). Fifty‐seven MEAs were associated with altitude and temperature measures while 358 were associated with rainfall measures. The most significant MEAs and outlier loci were used to identify genomic loci with adaptive potential (some already reported in wheat), including dormancy and frost resistance loci. We advocate the application of genomics and landscape genomics on ex situ collections of crop wild relatives to efficiently identify promising alleles and genetic materials for incorporation into modern crop breeding.     

Functional genomic analysis of corals from natural CO2‐seeps reveals core molecular responses involved in acclimatization to ocean acidification

Little is known about the potential for acclimatization or adaptation of corals to ocean acidification and even less about the molecular mechanisms underpinning these processes. Here, we examine global gene expression patterns in corals and their intracellular algal symbionts from two replicate population pairs in Papua New Guinea that have undergone long‐term acclimatization to natural variation in pCO₂. In the coral host, only 61 genes were differentially expressed in response to pCO₂ environment, but the pattern of change was highly consistent between replicate populations, likely reflecting the core expression homeostasis response to ocean acidification. Functional annotations highlight lipid metabolism and a change in the stress response capacity of corals as key parts of this process. Specifically, constitutive downregulation of molecular chaperones was observed, which may impact response to combined climate change‐related stressors. Elevated CO₂ has been hypothesized to benefit photosynthetic organisms but expression changes of in hospite Symbiodinium in response to acidification were greater and less consistent among reef populations. This population‐specific response suggests hosts may need to adapt not only to an acidified environment, but also to changes in their Symbiodinium populations that may not be consistent among environments, adding another challenging dimension to the physiological process of coping with climate change.     

Genomic analysis of a cardinalfish with larval homing potential reveals genetic admixture in the Okinawa Islands

Discrepancies between potential and observed dispersal distances of reef fish indicate the need for a better understanding of the influence of larval behaviour on recruitment and dispersal. Population genetic studies can provide insight on the degree to which populations are connected, and the development of restriction site‐associated sequencing (RAD‐Seq) methods has made such studies of nonmodel organisms more accessible. We applied double‐digest RAD‐Seq methods to test for population differentiation in the coral reef‐dwelling cardinalfish, Siphamia tubifer, which based on behavioural studies, have the potential to use navigational cues to return to natal reefs. Analysis of 11,836 SNPs from fish collected at coral reefs in Okinawa, Japan, from eleven locations over 3 years reveals little genetic differentiation between groups of S. tubifer at spatial scales from 2 to 140 km and between years at one location: pairwise F_ST values were between 0.0116 and 0.0214. These results suggest that the Kuroshio Current largely influences larval dispersal in the region, and in contrast to expectations based on studies of other cardinalfishes, there is no evidence of population structure for S. tubifer at the spatial scales examined. However, analyses of outlier loci putatively under selection reveal patterns of temporal differentiation that indicate high population turnover and variable larval supply from divergent source populations between years. These findings highlight the need for more studies of fishes across various geographic regions that also examine temporal patterns of genetic differentiation to better understand the potential connections between early life‐history traits and connectivity of reef fish populations.