Beekeeping practices and geographic distance, not land use, drive gene flow across tropical bees

Effects of landscape, resource use, and body size on genetic structure in bee populations

Quantifying genetic structure and levels of genetic variation are fundamentally important to predicting the ability of populations to persist in human-altered landscapes and adapt to future environmental changes. Genetic structure reflects the dispersal of individuals over generations, which can be mediated by species-level traits or environmental factors. Dispersal distances are commonly positively associated with body size and negatively associated with the amount of degraded habitat between sites, motivating the investigation of these potential drivers of dispersal concomitantly. We quantified genetic structure and genetic variability within populations of seven bee species from the genus Euglossa across fragmented landscapes. We genotyped bees at SNP loci and tested the following predictions: (1) deforested areas restrict gene flow; (2) larger species have lower genetic structure; (3) species with greater resource specialization have higher genetic structure; and (4) sites surrounded by more intact habitat have higher genetic diversity. Contrasting with previous work on bees, we found no associations between body size and genetic structure. Genetic structure was higher for species with greater resource specialization, and the amount of intact habitat between or surrounding sites was positively associated with parameters reflecting gene flow and genetic diversity. These results challenge the dominant paradigm that individuals of larger species disperse farther, and they suggest that landscape and resource requirements are important factors mediating dispersal.

Insights from Melipona bicolor hybrid genome assembly: a stingless bee genome with chromosome-level scaffold

BACKGROUND

The highly eusocial stingless bees are crucial pollinators of native and agricultural ecosystems. Nevertheless, genomic studies within this bee tribe remain scarce. We present the genome assembly of the stingless bee Melipona bicolor. This bee is a remarkable exception to the typical single-queen colony structure, since in this species, multiple queens may coexist and share reproductive duties, resulting in genetically diverse colonies with weak kinship connections. As the only known genuinely polygynous bee, M. bicolor's genome provides a valuable resource for investigating sociality beyond kin selection.

RESULTS

The genome was assembled employing a hybrid approach combining short and long reads, resulting in 241 contigs spanning 259 Mb (N50 of 6.2 Mb and 97.5% complete BUSCOs). Comparative analyses shed light on some evolutionary aspects of stingless bee genomics, including multiple chromosomal rearrangements in Melipona. Additionally, we explored the evolution of venom genes in M. bicolor and other stingless bees, revealing that, apart from two genes, the conserved repertoire of venom components remains under purifying selection in this clade.

CONCLUSION

This study advances our understanding of stingless bee genomics, contributing to the conservation efforts of these vital pollinators and offering insights into the evolutionary mechanisms driving their unique adaptations.

Limited introgression from non-native commercial strains and signatures of adaptation in the key pollinator Bombus terrestris

Insect pollination is fundamental for natural ecosystems and agricultural crops. The bumblebee species Bombus terrestris has become a popular choice for commercial crop pollination worldwide due to its effectiveness and ease of mass rearing. Bumblebee colonies are mass produced for the pollination of more than 20 crops and imported into over 50 countries including countries outside their native ranges, and the risk of invasion by commercial non-native bumblebees is considered an emerging issue for global conservation and biological diversity. Here, we use genome-wide data from seven wild populations close to and far from farms using commercial colonies, as well as commercial populations, to investigate the implications of utilizing commercial bumblebee subspecies in the UK. We find evidence for generally low levels of introgression between commercial and wild bees, with higher admixture proportions in the bees occurring close to farms. We identify genomic regions putatively involved in local and global adaptation, and genes in locally adaptive regions were found to be enriched for functions related to taste receptor activity, oxidoreductase activity, fatty acid and lipid biosynthetic processes. Despite more than 30 years of bumblebee colony importation into the UK, we observe low impact on the genetic integrity of local B. terrestris populations, but we highlight that even limited introgression might negatively affect locally adapted populations.

Shifting range in a stingless bee leads to pre-mating reproductive interference between species

Human-induced shifts in species’ ranges can increase contact between closely related species and lead to reproductive interference. In Australia, climate change and trade in stingless bee colonies is increasing the range overlap of two cryptic species: Tetragonula carbonaria and T. hockingsi. To investigate reproductive interactions between these species, we validated a diagnostic-PCR test based on the mitochondrial gene COI to ID field specimens to species. We then assessed the likelihood of reproductive interference in four ways. First, we imaged the male genitalia of each species and found no evidence of reproductive character displacement. Second, we assessed species composition of mating aggregations in an area of sympatry (Southeast Queensland) and confirmed that some males join the mating aggregations of interspecific colonies. Third, we translocated T. hockingsi colonies into the southern range of T. carbonaria (Sydney) and tracked their ability to requeen. These translocated colonies attracted mating aggregations comprised almost entirely of interspecific males, but never formed hybrid colonies; instead, queens either mated with their brothers, or the colony failed to requeen at all. Finally, we presented T. carbonaria males with either conspecific or interspecific virgin queens and found that males attempted to mate only with their own species’ queens. In all, we conclude that reproductive barriers between these species are complete with respect to “short-range” mating cues, but not for “long-range” mate attraction cues. Our study highlights that hive movements can increase some forms of pre-mating reproductive interference between managed bee species, even where the species do not actually mate or hybridize.

Evidence for morphological and genetic structuring of Plebeia flavocincta (Apidae: Meliponini) populations in Northeast Brazil

Tropical dry forests are increasingly threatened by human activities. In Northeast Brazil, the Caatinga, an area of tropical dry forests surrounded by wetlands, is considered globally unique among these habitats. In this region, the stingless bee Plebeia flavocincta is found in a variety of environmental, ecological, and demographic conditions. We aimed to characterize P. flavocincta populations within its natural range through wing geometric morphometrics and mitochondrial DNA (mtDNA) analyses. The characterization of population variability can clarify whether the species is morphologically and genetically diverse and whether populations are morphologically and genetically structured. We analyzed 673 samples by wing morphometry and 75 by cytochrome-b assays. Our results revealed P. flavocincta is genetically and morphologically diverse and populations are morphologically and genetically structured. Despite the differentiation between the two most morphologically distant populations, we verified a large overlap of morphological variation between all populations. The genetic analysis showed that the haplotypes were geographically structured into six clusters, four of which were located in coastal areas, and the remaining two in the drier inland region. The characterization of P. flavocincta populations is an important step toward decision-making in programs for the protection, management, and sustainable use of this stingless bee in local breeding efforts.

Males Are Capable of Long-Distance Dispersal in a Social Bee

Pollinator conservation is aided by knowledge of dispersal behavior, which shapes gene flow and population structure. In many bees, dispersal is thought to be male-biased, and males’ movements may be critical to maintaining gene flow in disturbed and fragmented habitats. Yet male bee movements are challenging to track directly and male dispersal ability remains poorly understood in most species. Here, we combine field manipulations and models to assess male dispersal ability in a stingless bee (Tetragonula carbonaria). We placed colonies with virgin queens at varying distances apart (1–48 km), genotyped the males that gathered at mating aggregations outside each colony, and used pairwise sibship assignment to determine the distribution of likely brothers across aggregations. We then compared simulations of male dispersal to our observed distributions and found best-fit models when males dispersed an average of 2–3 km (>2-fold female flight ranges), and maximum of 20 km (30-fold female flight ranges). Our data supports the view that male bee dispersal can facilitate gene flow over long-distances, and thus play a key role in bee populations’ resilience to habitat loss and fragmentation. In addition, we show that the number of families contributing to male aggregations can be used to estimate local stingless bee colony densities, allowing population monitoring of these important tropical pollinators.

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

Genetic structure across urban and agricultural landscapes reveals evidence of resource specialization and philopatry in the Eastern carpenter bee, Xylocopa virginica L

Human activity continues to impact global ecosystems, often by altering the habitat suitability, persistence, and movement of native species. It is thus critical to examine the population genetic structure of key ecosystemservice providers across human-altered landscapes to provide insight into the forces that limit wildlife persistence and movement across multiple spatial scales. While some studies have documented declines of bee pollinators as a result of human-mediated habitat alteration, others suggest that some bee species may benefit from altered land use due to increased food or nesting resource availability; however, detailed population and dispersal studies have been lacking. We investigated the population genetic structure of the Eastern carpenter bee, Xylocopa virginica, across 14 sites spanning more than 450 km, including dense urban areas and intensive agricultural habitat. X. virginica is a large bee which constructs nests in natural and human-associated wooden substrates, and is hypothesized to disperse broadly across urbanizing areas. Using 10 microsatellite loci, we detected significant genetic isolation by geographic distance and significant isolation by land use, where urban and cultivated landscapes were most conducive to gene flow. This is one of the first population genetic analyses to provide evidence of enhanced insect dispersal in human-altered areas as compared to semi-natural landscapes. We found moderate levels of regional-scale population structure across the study system (G'ST = 0.146) and substantial co-ancestry between the sampling regions, where co-ancestry patterns align with major human transportation corridors, suggesting that human-mediated movement may be influencing regional dispersal processes. Additionally, we found a signature of strong site-level philopatry where our analyses revealed significant levels of high genetic relatedness at very fine scales (<1 km), surprising given X. virginica's large body size. These results provide unique genetic evidence that insects can simultaneously exhibit substantial regional dispersal as well as high local nesting fidelity in landscapes dominated by human activity.

Minimal genetic differentiation of the malaria vector Nyssorhynchus darlingi associated with forest cover level in Amazonian Brazil

The relationship between deforestation and malaria in Amazonian Brazil is complex, and a deeper understanding of this relationship is required to inform effective control measures in this region. Here, we are particularly interested in characterizing the impact of land use and land cover change on the genetics of the major regional vector of malaria, Nyssorhynchus darlingi (Root). We used nextera-tagmented, Reductively Amplified DNA (nextRAD) genotyping-by-sequencing to genotype 164 Ny. darlingi collected from 16 collection sites with divergent forest cover levels in seven municipalities in four municipality groups that span the state of Amazonas in northwestern Amazonian Brazil: São Gabriel da Cachoeira, Presidente Figueiredo, four municipalities in the area around Cruzeiro do Sul, and Lábrea. Using a dataset of 5,561 Single Nucleotide Polymorphisms (SNPs), we investigated the genetic structure of these Ny. darlingi populations with a combination of model- and non-model-based analyses. We identified weak to moderate genetic differentiation among the four municipality groups. There was no evidence for microgeographic genetic structure of Ny. darlingi among forest cover levels within the municipality groups, indicating that there may be gene flow across areas of these municipalities with different degrees of deforestation. Additionally, we conducted an environmental association analysis using two outlier detection methods to determine whether individual SNPs were associated with forest cover level without affecting overall population genetic structure. We identified 14 outlier SNPs, and investigated functions associated with their proximal genes, which could be further characterized in future studies.

Population genomics of Bombus terrestris reveals high but unstructured genetic diversity in a potential glacial refugium

Ongoing climate change is expected to cause an increase in temperature and a reduction of precipitation levels in the Mediterranean region, which might cause changes in many species distributions. These effects negatively influence species gene pools, decreasing genetic variability and adaptive potential. Here, we use mitochondrial DNA and RADseq to analyse population genetic structure and genetic diversity of the bumblebee species Bombus terrestris (subspecies Bombus terrestris lusitanicus), in the Iberian Peninsula. Although this subspecies shows a panmictic pattern of population structure across Iberia and beyond, we found differentiation between subspecies B. t. lusitanicus and B. t. africanus, probably caused by the existence of barriers to gene flow between Iberia and North Africa. Furthermore, the results revealed that the Iberian Peninsula harbours a large fraction of B. terrestris intraspecific genetic variation, with the highest number of mitochondrial haplotypes found when compared with any other region in Europe studied so far, suggesting a potential role for the Iberian Peninsula as a glacial refugium. Our findings strengthen the idea that Iberia is a very important source of diversity for the global genetic pool of this species, because rare alleles might play a role in population resilience against human- or climate-mediated changes.

Habitat Loss Does Not Always Entail Negative Genetic Consequences

Although habitat loss has large, consistently negative effects on biodiversity, its genetic consequences are not yet fully understood. This is because measuring the genetic consequences of habitat loss requires accounting for major methodological limitations like the confounding effect of habitat fragmentation, historical processes underpinning genetic differentiation, time-lags between the onset of disturbances and genetic outcomes, and the need for large numbers of samples, genetic markers, and replicated landscapes to ensure sufficient statistical power. In this paper we overcame all these challenges to assess the genetic consequences of extreme habitat loss driven by mining in two herbs endemic to Amazonian savannas. Relying on genotyping-by-sequencing of hundreds of individuals collected across two mining landscapes, we identified thousands of neutral and independent single-nucleotide polymorphisms (SNPs) in each species and used these to evaluate population structure, genetic diversity, and gene flow. Since open-pit mining in our study region rarely involves habitat fragmentation, we were able to assess the independent effect of habitat loss. We also accounted for the underlying population structure when assessing landscape effects on genetic diversity and gene flow, examined the sensitivity of our analyses to the resolution of spatial data, and used annual species and cross-year analyses to minimize and quantify possible time-lag effects. We found that both species are remarkably resilient, as genetic diversity and gene flow patterns were unaffected by habitat loss. Whereas historical habitat amount was found to influence inbreeding; heterozygosity and inbreeding were not affected by habitat loss in either species, and gene flow was mainly influenced by geographic distance, pre-mining land cover, and local climate. Our study demonstrates that it is not possible to generalize about the genetic consequences of habitat loss, and implies that future conservation efforts need to consider species-specific genetic information.

Defining Pollinator Health: A Holistic Approach Based on Ecological, Genetic, and Physiological Factors

Evidence for global bee population declines has catalyzed a rapidly evolving area of research that aims to identify the causal factors and to effectively assess the status of pollinator populations. The term pollinator health emerged through efforts to understand causes of bee decline and colony losses, but it lacks a formal definition. In this review, we propose a definition for pollinator health and synthesize the available literature on the application of standardized biomarkers to assess health at the individual, colony, and population levels. We focus on biomarkers in honey bees, a model species, but extrapolate the potential application of these approaches to monitor the health status of wild bee populations. Biomarker-guided health measures can inform beekeeper management decisions, wild bee conservation efforts, and environmental policies. We conclude by addressing challenges to pollinator health from a One Health perspective that emphasizes the interplay between environmental quality and human, animal, and bee health.

Landscape genomics to the rescue of a tropical bee threatened by habitat loss and climate change

Habitat degradation and climate change are currently threatening wild pollinators, compromising their ability to provide pollination services to wild and cultivated plants. Landscape genomics offers powerful tools to assess the influence of landscape modifications on genetic diversity and functional connectivity, and to identify adaptations to local environmental conditions that could facilitate future bee survival. Here, we assessed range-wide patterns of genetic structure, genetic diversity, gene flow, and local adaptation in the stingless bee Melipona subnitida, a tropical pollinator of key biological and economic importance inhabiting one of the driest and hottest regions of South America. Our results reveal four genetic clusters across the species' full distribution range. All populations were found to be under a mutation-drift equilibrium, and genetic diversity was not influenced by the amount of reminiscent natural habitats. However, genetic relatedness was spatially autocorrelated and isolation by landscape resistance explained range-wide relatedness patterns better than isolation by geographic distance, contradicting earlier findings for stingless bees. Specifically, gene flow was enhanced by increased thermal stability, higher forest cover, lower elevations, and less corrugated terrains. Finally, we detected genomic signatures of adaptation to temperature, precipitation, and forest cover, spatially distributed in latitudinal and altitudinal patterns. Taken together, our findings shed important light on the life history of M. subnitida and highlight the role of regions with large thermal fluctuations, deforested areas, and mountain ranges as dispersal barriers. Conservation actions such as restricting long-distance colony transportation, preserving local adaptations, and improving the connectivity between highlands and lowlands are likely to assure future pollination services.

A trait-based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life-history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species-level life-history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G'_ST , and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G'_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.

Mitochondrial genome characterization of Melipona bicolor: Insights from the control region and gene expression data

The stingless bee Melipona bicolor is the only bee in which true polygyny occurs. Its mitochondrial genome was first sequenced in 2008, but it was incomplete and no information about its transcription was known. We combined short and long reads of M. bicolor DNA with RNASeq data to obtain insights about mitochondrial evolution and gene expression in bees. The complete genome has 15,001 bp, including a control region of 255 bp that contains all conserved structures described in honeybees with the highest AT content reported so far for bees (98.1%), displaying a compact but functional region. Gene expression control is similar to other insects however unusual patterns of expression may suggest the existence of different isoforms for the mitochondrially encoded 12S rRNA. Results reveal unique and shared features of the mitochondrial genome in terms of sequence evolution and gene expression making M. bicolor an interesting model to study mitochondrial genomic evolution.

Circuit-theory applications to connectivity science and conservation

Conservation practitioners have long recognized ecological connectivity as a global priority for preserving biodiversity and ecosystem function. In the early years of conservation science, ecologists extended principles of island biogeography to assess connectivity based on source patch proximity and other metrics derived from binary maps of habitat. From 2006 to 2008, the late Brad McRae introduced circuit theory as an alternative approach to model gene flow and the dispersal or movement routes of organisms. He posited concepts and metrics from electrical circuit theory as a robust way to quantify movement across multiple possible paths in a landscape, not just a single least-cost path or corridor. Circuit theory offers many theoretical, conceptual, and practical linkages to conservation science. We reviewed 459 recent studies citing circuit theory or the open-source software Circuitscape. We focused on applications of circuit theory to the science and practice of connectivity conservation, including topics in landscape and population genetics, movement and dispersal paths of organisms, anthropogenic barriers to connectivity, fire behavior, water flow, and ecosystem services. Circuit theory is likely to have an effect on conservation science and practitioners through improved insights into landscape dynamics, animal movement, and habitat-use studies and through the development of new software tools for data analysis and visualization. The influence of circuit theory on conservation comes from the theoretical basis and elegance of the approach and the powerful collaborations and active user community that have emerged. Circuit theory provides a springboard for ecological understanding and will remain an important conservation tool for researchers and practitioners around the globe.

Everything you always wanted to know about gene flow in tropical landscapes (but were afraid to ask)

The bulk of the world’s biodiversity is found in tropical regions, which are increasingly threatened by the human-led degradation of natural habitats. Yet, little is known about tropical biodiversity responses to habitat loss and fragmentation. Here we review all available literature assessing landscape effects on gene flow in tropical species, aiming to help unravel the factors underpinning functional connectivity in the tropics. We map and classify studies by focus species, the molecular markers employed, statistical approaches to assess landscape effects on gene flow, and the evaluated landscape and environmental variables. We then compare qualitatively and quantitatively landscape effects on gene flow across species and units of analysis. We found 69 articles assessing landscape effects on gene flow in tropical organisms, most of which were published in the last five years, were concentrated in the Americas, and focused on amphibians or mammals. Most studies employed population-level approaches, microsatellites were the preferred type of markers, and Mantel and partial Mantel tests the most common statistical approaches used. While elevation, land cover and forest cover were the most common gene flow predictors assessed, habitat suitability was found to be a common predictor of gene flow. A third of all surveyed studies explicitly assessed the effect of habitat degradation, but only 14 of these detected a reduced gene flow with increasing habitat loss. Elevation was responsible for most significant microsatellite-based isolation by resistance effects and a single study reported significant isolation by non-forested areas in an ant. Our study reveals important knowledge gaps on the study of landscape effects on gene flow in tropical organisms, and provides useful guidelines on how to fill them.

Functional connectivity in replicated urban landscapes in the land snail (Cornu aspersum)

Urban areas are highly fragmented and thereby exert strong constraints on individual dispersal. Despite this, some species manage to persist in urban areas, such as the garden snail, Cornu aspersum, which is common in cityscapes despite its low mobility. Using landscape genetic approaches, we combined study area replication and multiscale analysis to determine how landscape composition, configuration and connectivity influence snail dispersal across urban areas. At the overall landscape scale, areas with a high percentage of roads decreased genetic differentiation between populations. At the population scale, genetic differentiation was positively linked with building surface, the proportion of borders where wooded patches and roads appeared side by side and the proportion of borders combining wooded patches and other impervious areas. Analyses based on pairwise genetic distances validated the isolation-by-distance and isolation-by-resistance models for this land snail, with an equal fit to least-cost paths and circuit-theory-based models. Each of the 12 landscapes analysed separately yielded specific relations to environmental features, whereas analyses integrating all replicates highlighted general common effects. Our results suggest that urban transport infrastructures facilitate passive snail dispersal. At a local scale, corresponding to active dispersal, unfavourable habitats (wooded and impervious areas) isolate populations. This work upholds the use of replicated landscapes to increase the generalizability of landscape genetics results and shows how multiscale analyses provide insight into scale-dependent processes.

Landscape Genomic Conservation Assessment of a Narrow-Endemic and a Widespread Morning Glory From Amazonian Savannas

Although genetic diversity ultimately determines the ability of organisms to adapt to environmental changes, conservation assessments like the widely used International Union for Conservation of Nature (IUCN) Red List Criteria do not explicitly consider genetic information. Including a genetic dimension into the IUCN Red List Criteria would greatly enhance conservation efforts, because the demographic parameters traditionally considered are poor predictors of the evolutionary resilience of natural populations to global change. Here we perform the first genomic assessment of genetic diversity, gene flow, and patterns of local adaptation in tropical plant species belonging to different IUCN Red List Categories. Employing RAD-sequencing we identified tens of thousands of single-nucleotide polymorphisms in an endangered narrow-endemic and a least concern widespread morning glory (Convolvulaceae) from Amazonian savannas, a highly threatened and under-protected tropical ecosystem. Our results reveal greater genetic diversity and less spatial genetic structure in the endangered species. Whereas terrain roughness affected gene flow in both species, forested and mining areas were found to hinder gene flow in the endangered plant. Finally we implemented environmental association tests and genome scans for selection, and identified a higher proportion of candidate adaptive loci in the widespread species. These mainly contained genes related to pathogen resistance and physiological adaptations to life in nutrient-limited environments. Our study emphasizes that IUCN Red List Criteria do not always prioritize species with low genetic diversity or whose genetic variation is being affected by habitat loss and fragmentation, and calls for the inclusion of genetic information into conservation assessments. More generally, our study exemplifies how landscape genomic tools can be employed to assess the status, threats and adaptive responses of imperiled biodiversity.