Find research datasets worth reusing

Data from: Population structure, genetic connectivity, and adaptation in the Olympia oyster (Ostrea lurida) along the west coast of North America

Effective management of threatened and exploited species requires an understanding of both the genetic connectivity among populations and local adaptation. The Olympia oyster (Ostrea lurida), patchily distributed from Baja California to the central coast of Canada, has a long history of population declines due to anthropogenic stressors. For such coastal marine species, population structure could follow a continuous isolation-by-distance model, contain regional blocks of genetic similarity separated by barriers to gene flow, or be consistent with a null model of no population structure. To distinguish between these hypotheses in O. lurida, 13,424 single-nucleotide polymorphisms (SNPs) were used to characterize rangewide population structure, genetic connectivity, and adaptive divergence. Samples were collected across the species range on the west coast of North America, from southern California to Vancouver Island. A conservative approach for detecting putative loci under selection identified 235 SNPs across 129 GBS loci, which were functionally annotated and analyzed separately from the remaining neutral loci. While strong population structure was observed on a regional scale in both neutral and outlier markers, neutral markers had greater power to detect fine-scale structure. Geographic regions of reduced gene flow aligned with known marine biogeographic barriers, such as Cape Mendocino, Monterey Bay, and the currents around Cape Flattery. The outlier loci identified as under putative selection included genes involved in developmental regulation, sensory information processing, energy metabolism, immune response, and muscle contraction. These loci are excellent candidates for future research and may provide targets for genetic monitoring programs. Beyond specific applications for restoration and management of the Olympia oyster, this study lends to the growing body of evidence for both population structure and adaptive differentiation across a range of marine species exhibiting the potential for panmixia. Computational notebooks are available to facilitate reproducibility and future open-sourced research on the population structure of Ostrea lurida.

opencc-zeroDec 2019View details →

Population structure of five native sheep breeds of Sweden estimated with high density SNP genotypes

Background Native Swedish sheep breeds are part of the North European short-tailed sheep group; characterized in part by their genetic uniqueness. Our objective was to study the population structure of native Swedish sheep. Five breeds were genotyped using the 600 K SNP array. Dalapäls and Klövsjö sheep are from the middle of Sweden; Gotland and Gute sheep from Gotland, an island in the Baltic Sea; and Fjällnäs sheep from northern Sweden. We studied population structure by: principal component analysis (PCA), cluster-based analysis of admixture, and an estimated population tree. Results The analyses of the five Swedish breeds revealed that these breeds are five distinct breeds, while Gute and Gotland are more closely related to each other as seen in all analyses. All breeds had long branch lengths in the population tree indicating they've been subjected to drift. We repeated our analyses using 39 K SNP and including 50 K SNP genotypes from other European and southwestern Asian breeds from the Sheep HapMap project and 600 K SNP genotypes from a dataset of French sheep. Results arranged breeds into five groups: south-west Asia, south-west Europe, central Europe, north Europe and north European short-tailed sheep. Within this last group, Norwegian and Icelandic breeds, Finn and Romanov sheep, Scottish breeds, and Gute and Gotland sheep were more closely related while the remaining Swedish breeds and Ouessant sheep were distinct from all breeds and had longer branches in the population tree. Conclusions We showed population structure of five Swedish breeds and their structure within European and southwestern Asian breeds. Swedish breeds are unique, distinct breeds that have been subjected to drift but group with other north European short-tailed sheep.

Data from: Genetic diagnosis of a rare myrmecochorous species, Plagiorhegma dubium (Berberidaceae): historical genetic bottlenecks and strong spatial structures among populations

Distribution of genetic variation over time and space is relevant to demographic histories, and tightly linked to ecological disturbances as well as evolutionary potential of an organism. Therefore, understanding the pattern of genetic diversity is a primary step in conservation and management projects for rare and threatened plant species. We used 8 microsatellite markers to examine the level of genetic diversity, spatial structure and demographic history of Plagiorhegma dubium, a rare myrmecochorous herb, populations sampled across northeast Asia and Siberia. We found low within-population genetic variation associated with historical bottlenecks. Although pairwise FST values were not much higher than the ones found in similar life form species, STRUCTURE and PCoA revealed a clear broad-scale spatial pattern of genetic structure. Bayesian clustering (best K=6) and PCoA identified three populations that are distinctive from neighboring populations in the Korean peninsula, which suggests potential units for conservation and management plans in Korea. MIGRATE-N and BAYESASS showed that both contemporary (0.003-0.045) and historical migration rates (2e-5-4.6e-4) were low. Our findings provide a good example, where genetic considerations should be integrated for conservation and management plans of rare and threatened species.

Data from: Phoretic dispersal influences parasite population genetic structure

Dispersal is a fundamental component of the life history of most species. Dispersal influences fitness, population dynamics, gene flow, genetic drift, and population genetic structure. Even small differences in dispersal can alter ecological interactions and trigger an evolutionary cascade. Linking such ecological processes with evolutionary patterns is difficult, but can be done in the proper comparative context. Here we investigate how differences in phoretic dispersal influence the population genetic structure of two different parasites of the same host species. We focus on two species of host-specific feather lice (Phthiraptera: Ischnocera) that co-occur on feral rock pigeons (Columba livia). Although these lice are ecologically very similar, "wing lice" (Columbicola columbae) disperse phoretically by "hitchhiking" on pigeon flies (Diptera: Hippoboscidae), while "body lice" (Campanulotes compar) do not. Differences in the phoretic dispersal of these species are thought to underlie observed differences in host specificity, as well as the degree of host-parasite cospeciation. These ecological and macroevolutionary patterns suggest that body lice should exhibit more genetic differentiation than wing lice. We tested this prediction among lice on individual birds and among lice on birds from three pigeon flocks. We found higher levels of genetic differentiation in body lice compared to wing lice at two spatial scales. Our results indicate that differences in phoretic dispersal can explain microevolutionary differences in population genetic structure and are consistent with macroevolutionary differences in the degree of host-parasite cospeciation.

opencc-zeroJan 2020View details →

Genetic analyses reveal population structure and recent decline in leopards (Panthera pardus fusca) across Indian subcontinent

Background Large carnivores maintain the stability and functioning of ecosystems. Currently, many carnivore species face declining population sizes due to natural and anthropogenic pressures. The leopard, Panthera pardus, is probably the most widely distributed and highly adaptable large felid globally, still persisting in most of its historic range. However, we lack subspecies-level data on country or regional scale on population trends, as ecological monitoring approaches are difficult to apply on such wide-ranging species. We used genetic data from leopards sampled across the Indian subcontinent to investigate population structure and patterns of demographic decline. Methods We collected faecal samples from the Terai-Arc landscape of north India and identified 56 unique individuals using a panel of 13 microsatellite markers. We merged this data with already available 143 leopard individuals and assessed genetic structure at country scale. Subsequently, we investigated the demographic history of each identified subpopulations and compared genetic decline analyses with countrywide local extinction probabilities. Results Our genetic analyses revealed four distinct subpopulations corresponding to Western Ghats, Deccan Plateau-Semi Arid, Shivalik and Terai region of the north Indian landscape, each with high genetic variation. Coalescent simulations with microsatellite loci revealed a possibly human-induced 75-90% population decline between ∼120-200 years ago across India. Population-specific estimates of genetic decline are in concordance with ecological estimates of local extinction probabilities in these subpopulations obtained from occupancy modeling of the historic and current distribution of leopards in India. Conclusions Our results confirm the population decline of a widely distributed, adaptable large carnivore. We re-iterate the relevance of indirect genetic methods for such species in conjunction with occupancy assessment and recommend that detailed, landscape-level ecological studies on leopard populations are critical to future conservation efforts. Our approaches and inference are relevant to other widely distributed, seemingly unaffected carnivores such as the leopard.

Data from: Intraspecific DNA contamination distorts subtle population structure in a marine fish: decontamination of herring samples before restriction-site associated (RAD) sequencing and its effects on population genetic statistics

Wild specimens are often collected in challenging field conditions, where samples may be contaminated with the DNA of conspecific individuals. This contamination can result in false genotype calls, which are difficult to detect, but may also cause inaccurate estimates of heterozygosity, allele frequencies, and genetic differentiation. Marine broadcast spawners are especially problematic, because population genetic differentiation is low and samples are often collected in bulk and sometimes from active spawning aggregations. Here, we used contaminated and clean Pacific herring (Clupea pallasi) samples to test (i) the efficacy of bleach decontamination, (ii) the effect of decontamination on RAD genotypes, and (iii) the consequences of contaminated samples on population genetic analyses. We collected fin tissue samples from actively spawning (and thus contaminated) wild herring and non-spawning (uncontaminated) herring. Samples were soaked for 10 minutes in bleach or left untreated, and extracted DNA was used to prepare DNA libraries using a restriction-site associated DNA (RAD) approach. Our results demonstrate that intraspecific DNA contamination affects patterns of individual and population variability, causes an excess of heterozygotes, and biases estimates of population structure. Bleach decontamination was effective at removing intraspecific DNA contamination and compatible with RAD sequencing, producing high-quality sequences, reproducible genotypes, and low levels of missing data. Although sperm contamination may be specific to broadcast spawners, intraspecific contamination of samples may be common and difficult to detect from high-throughput sequencing data, and can impact downstream analyses.

opencc-zeroDec 2015View details →

Data from: Population genetic structure between Yap and Palau for the coral Acropora hyacinthus

Information on connectivity is becoming increasingly in demand as marine protected areas are being designed as an integral part of a network to protect marine resources at the ecosystem level. Larval dispersal and population structure, however, remain very difficult to assess. Here, we tested the predictions of a detailed oceanographic connectivity model of larval dispersal and coral recruitment within Palau and between Palau and Yap, which was developed to support the review of the existing network of marine protected areas in Palau. We used high throughput microsatellite genotyping of the coral Acropora hyacinthus to characterize population genetic structure. Pairwise F′ST values between Palau and Yap (0.10), Palau and Ngulu (0.09) and Yap and Ngulu (0.09) were all significant and similar to pairwise F′ST values of sites within Palau (0.02–0.12) and within Yap (0.02–0.09) highlighting structure at island scale and indicating that recruitment may be even more localized than previously anticipated. A bottleneck test did not reveal any signs of a founder effect between Yap and Palau. Overall, the data supports the idea that recovery of A. hyacinthus in Palau did not come exclusively from a single source but most likely came from a combination of areas, including sites within Palau. In light of these results there seems to be very little connectivity around the barrier reef and management recommendation would be to increase the number or the size of MPAs within Palau.

opencc-zeroOct 2019View details →

Population structure in landrace barley (Hordeum vulgare L.) during the late 19th century crop failures in Fennoscandia

Agricultural disasters and the subsequent need for supply of relief seed can be expected to influence the genetic composition of crop plant populations. The consequences of disasters and seed relief have, however, rarely been studied since specimens sampled before the events are seldomly available. A series of crop failures struck northern Fennoscandia (Norway, Sweden and Finland) during the second half of the 19th century. In order to assess population genetic dynamics of landrace barley (Hordeum vulgare), and consequences of crop failure and possible seed relief during this time period, we genotyped seeds from 16 historical accessions originating from two time periods spanning the period of repeated crop failure. Reliable identification of genetic structuring is highly dependent on sampling regimes and detecting fine-scale geographic or temporal differentiation requires large sample sizes. The robustness of the results under different sampling regimes was evaluated by analyzing subsets of the data and an artificially pooled dataset. The results led to the conclusion that six individuals per accession were insufficient for reliable detection of the observed genetic structure. We found that population structure among the data was best explained by collection year of accessions, rather than geographic origin. The correlation with collection year indicated a change in genetic composition of landrace barley in the area after repeated crop failures, likely a consequence of introgression of relief seed in local populations. Identical genotypes were found to be shared among some accessions, suggesting founder effects and local seed exchange along known routes for trade and cultural exchange.

opencc-zeroDec 2012View details →

Data from: Variation in population structure and dynamics of montane forest tree species in Ethiopia guide priorities for conservation and research

The greatest extent of Afromontane environments in the world is found in Ethiopia. These areas support exceptional biodiversity, but forest cover and ecological integrity have declined sharply in recent decades. Conservation and management efforts are hampered in part by an inadequate understanding of the basic ecology of major tree species. We investigated population structure and inferred population dynamics from size frequency distributions of 22 forest tree species encountered in montane forests of Ethiopia. We collected new empirical data from four sites in the Bale Mountains, where some of the country's most extensive and least disturbed forests remain, and conducted a systematic review and analysis of all such studies that reported population structure for one or more of these species in Ethiopia. Thirteen widespread montane tree species showed a reverse-J size distribution, indicating a relatively stable population structure. Six other species had size-frequency distributions that indicate episodic recruitment and/or removal of certain size classes. Specific causes of these patterns are uncertain: they may involve timber harvesting, herbivory, fire, or natural disturbances, but patterns were inconsistent and locality dependent. For three other tree species, existing data are inadequate for any interpretation of population structure and dynamics. A species of particular conservation concern that emerged from this analysis was Hagenia abyssinica, which was found in all areas to consist only of larger individuals with no recent recruitment. For management and conservation purposes, the species in most urgent need of new research are those with inadequate or inconsistent data, and H. abyssinica..

Data from: Inferring continuous and discrete population genetic structure across space

A classic problem in population genetics is the characterization of discrete population structure in the presence of continuous patterns of genetic differentiation. Especially when sampling is discontinuous, the use of clustering or assignment methods may incorrectly ascribe differentiation due to continuous processes (e.g., geographic isolation by distance) to discrete processes, such as geographic, ecological, or reproductive barriers between populations. This reflects a shortcoming of current methods for inferring and visualizing population structure when applied to genetic data deriving from geographically distributed populations. Here, we present a statistical framework for the simultaneous inference of continuous and discrete patterns of population structure. The method estimates ancestry proportions for each sample from a set of two-dimensional population layers, and, within each layer, estimates a rate at which relatedness decays with distance. This thereby explicitly addresses the "clines versus clusters" problem in modeling population genetic variation, and remedies some of the overfitting to which nonspatial models are prone. The method produces useful descriptions of structure in genetic relatedness in situations where separated, geographically distributed populations interact, as after a range expansion or secondary contact. We demonstrate the utility of this approach using simulations and by applying it to empirical datasets of poplars and black bears in North America.

Data from: Genetic structure of the Painted Bunting and its implications for conservation of migratory populations

The Painted Bunting Passerina ciris is a Neotropical songbird which breeds primarily in the United States during the summer and migrates to Mexico, Central America, southern Florida, and the Caribbean over the winter. Male Painted Buntings are brightly coloured, which makes them highly sought after as pets, particularly in Mexico, Central America and Europe. We used short sequence repeats (microsatellite DNA) to investigate the population genetic structure of the Painted Bunting and its implications in conservation management of migratory populations. We found a detectable level of population differentiation as revealed by pairwise FST and RST comparisons and Bayesian clustering analyses, with strong support for differentiation between eastern and western Painted Buntings (e.g. Oklahoma and Georgia FST = 0.1; P = 0.005; RST = 0.18; P = 0.04) in accordance with previous mitochondrial DNA analysis. We recovered additional support for two sub‐groups within the western clade. While linking migrant songbirds captured outside of the United States to their breeding populations remains a challenge, we show that natural levels of population genetic differentiation can be detected via microsatellite DNA markers and exploited in migratory connectivity studies. We also demonstrate the potential utility of our low‐cost markers for population identification of birds recovered from the pet trade by screening a small subset of samples (n = 5) collected as part of wildlife tracking. We discuss the implications of our results for future efforts to understand patterns of population decline in Painted Buntings more generally, as well as how we might expand this methodology to combat illegal pet‐trade activity in this and other songbird species.

opencc-zeroOct 2023View details →

Of Mojave milkweed and mirrors: The population genomic structure of a species impacted by solar energy development

A rapid renewable energy transition has facilitated the development of large, ground‐mounted solar energy facilities worldwide. Deserts, and other sensitive aridland ecosystems, are the second most common land‐cover type for solar energy development globally. Thus, it is necessary to understand existing diversity within environmentally sensitive desert plant populations to understand spatiotemporal effects of solar energy siting and design. Overall, few population genomic studies of desert plants exist, and much of their biology is unknown. To help fill this knowledge gap, we sampled Mojave milkweed (Asclepias nyctaginifolia) in and around the Ivanpah Solar Electric Generating Station (ISEGS) in the Mojave Desert of California to understand the species' population structure, standing genetic variation, and how that intersects with solar development. We performed Restriction‐site Associated Sequencing (RADseq) and discovered 9942 single nucleotide polymorphisms (SNPs). Using these data, we found clear population structure over small spatial scales, suggesting each site sampled comprised a genetically distinct population of Mojave milkweed. While mowing, in lieu of blading, the vegetation across the solar energy facility's footprint prevented the immediate loss of the ISEGS Mojave milkweed population, we show that the effects of land‐cover change, especially those impacting desert washes, may impact long‐term genetic diversity and persistence. Potential implications of this include a risk of overall loss of genetic diversity, or even hastened extirpation. These findings highlight the need to consider the genetic diversity of impacted species when predicting the impact and necessary conservation measures of large‐scale land‐cover changes on species with small population sizes.

opencc-zeroNov 2023View details →

Data from: Relict stands of Central European oaks: unravelling autochthony and genetic structure based on a multi-population study

Central European white oaks expanded rapidly after the last glacial period and reached their current distribution range during the early Holocene. They have been an important resource of timber, fuelwood and animal feed for humans, who actively promoted their presence in forests and other landscape types at least since the early historical times. Besides stands with intensive management, putatively relict populations of three native oak species can be found on unproductive sites with restricted accessibility. Here, we apply chloroplast and nuclear microsatellite markers in order to address the autochthony of relict and managed stands and compare the spatial distribution of genetic variation between them. Based on data from more than 150 populations, we demonstrate that oak autochthony was preserved throughout historical times which is likely the result of traditional silvicultural treatment. This is supported by the fact that the spatial pattern of chloroplast haplotype distribution still reflects the post-glacial recolonization in both relict and old managed stands. We observed significant admixture of haplotypes only in stands established after the Second World War, which is attributable to the transfer of reproductive material used for afforestation. In terms of nuclear genetic variation, we observed marked differences among species. Quercus pubescens exhibited a pronounced genetic structure. Genetic drift and limited gene flow among its small and isolated populations in our study area might have contributed to this pattern. Varying extent of genetic introgression with other sympatric oak species could offer an additional explanation. On the contrary, the gene pools of Q. petraea and Q. robur are highly homogenous, displaying only weak isolation-by-distance. We found no significant differences of genetic diversity and differentiation between relict and managed stands. This suggests that seed transfer mostly occurred within our study area, even in those stands established in post-war times, verifying previous findings that point out limited human interference. We recommend consideration of population genetic structure for gene conservation, with a finer resolution of gene conservation units needed for Q. pubescens due to its spatial genetic structure. Both relict and old managed stands, species-pure or mixed, are suitable for conservation, as they host autochthonous gene pools. Coppice-with-standard management could contribute to preservation of autochthony. In the face of climate change, it is also important to maintain the evolutionary potential of the stands, by facilitating generative reproduction and allowing for hybridization in mixed stands.

opencc-zeroJan 2024View details →

Microgeographic population structuring in a genus of California trapdoor spiders and discovery of an enigmatic new species (Euctenizidae: Promyrmekiaphila korematsui sp. nov.)

The recognition and delineation of cryptic species remains a perplexing problem in systematics, evolution, and species delimitation. Once recognized as such, cryptic species complexes provide fertile ground for studying genetic divergence within the context of phenotypic and ecological divergence (or lack thereof). Herein we document the discovery of a new cryptic species of trapdoor spider, Promyrmekiaphila korematsui sp. nov. Using subgenomic data obtained via target enrichment, we document the phylogeography of the California endemic genus Promyrmekiaphila and its constituent species, which also includes P. clathrata and P. winnemem. Based on these data we show a pattern of strong geographic structuring among populations but cannot entirely discount recent gene flow among populations that are parapatric, particularly for deeply diverged lineages within P. clathrata. The genetic data, in addition to revealing a new undescribed species, also allude to a pattern of potential phenotypic differentiation where species likely come into contact. Alternatively, phenotypic cohesion among genetically divergent P. clathrata lineages suggests that some level of gene flow is ongoing or occurred in the recent past. Despite considerable field collection efforts over many years, additional sampling in potential zones of contact for both species and lineages is needed to completely resolve the dynamics of divergence in Promyrmekiaphila at the population-species interface.

Calling structural variants with confidence from short-read data in wild bird populations

Comprehensive characterisation of structural variation in natural populations has only become feasible in the last decade. To investigate the population genomic nature of structural variation (SV), reproducible and high-confidence SV callsets are first required. We created a population-scale reference of the genome-wide landscape of structural variation across 33 Nordic house sparrows (Passer domesticus) individuals. To produce a consensus callset across all samples using short-read data, we compare heuristic-based quality filtering and visual curation (Samplot/PlotCritic and Samplot-ML) approaches. We demonstrate that curation of SVs is important for reducing putative false positives and that the time invested in this step outweighs the potential costs of analysing short-read discovered SV datasets that include many potential false positives. We find that even a lenient manual curation strategy (e.g. applied by a single curator) can reduce the proportion of putative false positives by up to 80%, thus enriching the proportion of high-confidence variants. Crucially, in applying a lenient manual curation strategy with a single curator, nearly all (>99%) variants rejected as putative false positives were also classified as such by a more stringent curation strategy using three additional curators. Furthermore, variants rejected by manual curation failed to reflect the expected population structure from SNPs, whereas variants passing curation did. Combining heuristic-based quality-filtering with rapid manual curation of structural variants in short-read data can therefore become a time- and cost-effective first step for functional and population genomic studies requiring high-confidence SV callsets.

Population structure and connectivity among coastal and freshwater Kelp Gull (Larus dominicanus) populations from Patagonia

The genetic identification of significant evolutionary units and information on their connectivity can be used to design effective management and conservation plans. Despite having high dispersal capacity, several seabird species show population structure due to both abiotic and biotic barriers to gene flow. The Kelp Gull is the most abundant species of gull in the southern hemisphere. In Argentina it reproduces in both marine and freshwater environments, with more than 100,000 pairs following a metapopulation dynamic across 140 colonies in the Atlantic coast of Patagonia. However, little is known about the demography and connectivity of inland populations. We aim to provide information on the connectivity of the largest freshwater colonies (those from Nahuel Huapi Lake) with the closest Pacific and Atlantic populations to evaluate if these freshwater colonies are being subsidized by the larger coastal populations. We sampled three geographic regions (Nahuel Huapi Lake and the Atlantic and Pacific coasts) and employed a reduced-representation genomic approach to genotype individuals for single-nucleotide polymorphisms (SNPs). We found, using clustering and phylogenetic analyses, that there are three genetic groups, each corresponding to one of our sampled regions. Individuals from marine environments are more closely related to each other than to those from Nahuel Huapi Lake, indicating that the latter population constitutes the first freshwater Kelp Gull colony to be identified as a significant evolutionary unit in Patagonia.

Data from: Population structure and species delimitation in the Wehrle's salamander complex

Species are the fundamental unit of biodiversity studies. However, many species complexes are difficult to delimit, especially those characterized by complicated patterns of population structure. Salamanders in the family Plethodontidae often form species by slowly fragmenting across a landscape over space and time. They thus provide many examples of species complexes in which gradual Darwinian evolution has resulted in multiple units of varying degrees of differentiation, including incompletely separated lineages. Here we report on a molecular systematic investigation of woodland salamanders in the Plethodon wehrlei group, which has recently been split from two species into five. To quantify patterns of genetic variation, we collected genetic samples from 24 individuals from 20 populations, including all species and representing a carefully selected subset of previous work. From these samples, we obtained genomic data using anchored hybrid enrichment, which resulted in 319 loci averaging 1300 base pairs in length. Biallelic single nucleotide polymorphisms (SNPs) were randomly selected from 316 of these loci for some analyses. We examined patterns of genetic structure using PCA, DAPC, FEEMS, and STRUCTURE, and found that all of the recognized species formed genetic clusters; however, P. wehrlei and P. punctatus were relatively weakly differentiated, and STRUCTURE identified three separate clusters within P. jacksoni. Species trees inferred using wASTRAL, BPP, and TreeMix all recovered the same topology, with P. dixi sister to the other taxa, which included a northern clade (P. wehrlei, P. punctatus, P. pauleyi) and a southern clade (P. jacksoni, with three separate groups). TreeMix only inferred one gene flow event. We evaluated the candidate species using BPP and the genealogical divergence index (gdi). While BPP delimited all candidate species with strong support (all posterior probabilities = 1.0), the gdi only strongly supported P. dixi and P. pauleyi, both of which have only been recently described. We discuss the difficult problem of species delimitation in groups that form species via range fragmentation. We also provide a vision for future research with the aim of better testing and diagnosing the species diversity within the P. wehrlei group.

Data from: Range-wide genetic analysis of an endangered bumble bee (Bombus affinis) reveals population structure, isolation by distance, and low colony abundance

Declines in bumblebee species ranges and abundances are documented across multiple continents and have prompted the need for research to aid species recovery and conservation. The rusty patched bumblebee (Bombus affinis) is the first federally-listed bumblebee species in North America. We conducted a range-wide population genetics study of B. affinis from across all extant conservation units to inform conservation efforts. To understand the species' vulnerability and help establish recovery targets, we examined population structure, patterns of genetic diversity, and population differentiation. Additionally, we conducted site-level analysis of colony abundance to inform prioritizing areas for conservation, translocation, and other recovery actions. We find substantial evidence of population structuring along an east-to-west gradient. Putative populations show evidence of isolation by distance, high inbreeding coefficients, and a range wide male diploidy rate of ~15%. Our results suggest the Appalachians represents a genetically distinct cluster with high levels of private alleles and substantial differentiation from the rest of the extant range. Site-level analyses suggest low colony abundance estimates for B. affinis compared to similar datasets of stable, co-occurring species. These results lend genetic support to trends from observational studies suggesting B. affinis has undergone a recent decline and exhibits substantial spatial structure. The low colony abundances observed here suggest caution in overinterpreting the stability of populations even where B. affinis is reliably detected interannually. These results help delineate informed management units, provide context for the potential risks of translocation programs, and can help set clear recovery targets for this and other threatened bumblebee species.

Data from: Landscape heterogeneity drives population structure in four western bumble bee species

Bumble bees are critical pollinators in wild, agricultural, and urban ecosystems—providing the necessary ecological services for food and crop production. In western North America, mountain ranges have high bumble bee species richness. However, as climate change increases temperatures and restricts montane populations to higher elevational spaces, their ability to disperse and maintain genetic diversity decreases. This genetic isolation could lead to the extirpation of local pollinator communities and an overall loss of pollinators. We analyzed the genetic diversity of four broadly sympatric species of bumble bees across the Rocky and Cascade Mountains of western North America to assess habitat isolation's impact on population genetic structure. We expected species restricted to higher elevation habitats to display higher population structure and less genetic diversity across the landscape. We sampled approximately 150 bees per species from seven to eight sites across each species' range. We genotyped bees with an average of 10 loci and used FST and Bayesian Structure analysis to quantify population differentiation. Using isolation by distance and isolation by resistance analyses, species with both narrow and broad habitat suitability requirements showed evidence of habitat suitability restricting gene flow. Although each species showed varying degrees of genetic structure and gene flow, knowing how habitat heterogeneity drives genetic structure and isolation can help guide conservation efforts and determine focal regions for bumble bee conservation in the face of climate change.