Landscape genetics of a tropical rescue pollinator

Stingless bee Trigona spinipes (Hymenoptera: Apidae) behavior on chayote flowers (Sechium edule)

This experiment aimed to evaluate the forage species and their effect on fruit production of the chayote crop. For this, the culture was under observation in the first ten minutes at each time, from 7:00 am to 6:00 pm, with three repetitions, in each year studied (1994, 2001, 2009 and 2016). With the exception of 2009, Trigona spinipes stingless bees were frequent and constant insects in these flowers between 8:00 am and 1:00 pm, visiting mainly male flowers for nectar collection and due to foraging behavior this stingless bee may be considered a pollinator of the chayote.

Biomonitoring via DNA metabarcoding and light microscopy of bee pollen in rainforest transformation landscapes of Sumatra

BACKGROUND

Intense conversion of tropical forests into agricultural systems contributes to habitat loss and the decline of ecosystem functions. Plant-pollinator interactions buffer the process of forest fragmentation, ensuring gene flow across isolated patches of forests by pollen transfer. In this study, we identified the composition of pollen grains stored in pot-pollen of stingless bees, Tetragonula laeviceps, via dual-locus DNA metabarcoding (ITS2 and rbcL) and light microscopy, and compared the taxonomic coverage of pollen sampled in distinct land-use systems categorized in four levels of management intensity (forest, shrub, rubber, and oil palm) for landscape characterization.

RESULTS

Plant composition differed significantly between DNA metabarcoding and light microscopy. The overlap in the plant families identified via light microscopy and DNA metabarcoding techniques was low and ranged from 22.6 to 27.8%. Taxonomic assignments showed a dominance of pollen from bee-pollinated plants, including oil-bearing crops such as the introduced species Elaeis guineensis (Arecaceae) as one of the predominant taxa in the pollen samples across all four land-use types. Native plant families Moraceae, Euphorbiaceae, and Cannabaceae appeared in high proportion in the analyzed pollen material. One-way ANOVA (p > 0.05), PERMANOVA (R² values range from 0.14003 to 0.17684, for all tests p-value > 0.5), and NMDS (stress values ranging from 0.1515 to 0.1859) indicated a lack of differentiation between the species composition and diversity of pollen type in the four distinct land-use types, supporting the influx of pollen from adjacent areas.

CONCLUSIONS

Stingless bees collected pollen from a variety of agricultural crops, weeds, and wild plants. Plant composition detected at the family level from the pollen samples likely reflects the plant composition at the landscape level rather than the plot level. In our study, the plant diversity in pollen from colonies installed in land-use systems with distinct levels of forest transformation was highly homogeneous, reflecting a large influx of pollen transported by stingless bees through distinct land-use types. Dual-locus approach applied in metabarcoding studies and visual pollen identification showed great differences in the detection of the plant community, therefore a combination of both methods is recommended for performing biodiversity assessments via pollen identification.

What Can Restoration Do for Bee Communities? An Example in the Atlantic Rainforest in Paraná State, Southern Brazil

Biodiversity conservation is a challenge for today. Studies regarding different ecosystems have become extremely important for understanding communities and promoting strategies for conservation, especially those involving forest restoration as strategy to reverse biodiversity loss. We compared bee diversity indices between three remnants of Atlantic rainforest in southern Brazil and four adjacent areas that were reforested after shale mining, and which are now under different restoration levels. Seven sites were monitored for over 5 years (2011-2016), by sampling bees directly on flowers or in flight using an entomological net, with 400 collected individuals/site/year. Bee species composition differed between post-mining ages and between sites. In all, we sampled 14.185 specimens and 236 bee species. The introduced Africanized Apis mellifera Linnaeus was the most frequent and abundant species, followed by Trigona spinipes and Psaenythia bergii Holmberg. Among habitats, the reforested area in initial phase showed lower richness and diversity in relation all others sample sites. Conversely, all indices were higher in forest remnants, middle phase II, and advanced phase reforested areas, reinforcing the importance of reforestation for conservation, notably endangered species, such as Oxytrigona sp., Schwarziana quadripunctata (Lepeletier), and the solitary species of the genus Hylaeus, all found in the restored areas. These results represent an important contribution for understanding the recovery of the bee fauna in restored mining habitats. The dataset reveals an interesting response in areas that were mined for shale extraction and are now undergoing different levels of restoration, suggesting that older reforested habitats have a higher probability of having a fully recovered bee community.

Knowledge gaps hamper understanding the relationship between fragmentation and biodiversity loss: the case of Atlantic Forest fruit-feeding butterflies

Background

A key challenge for conservation biology in the Neotropics is to understand how deforestation affects biodiversity at various levels of landscape fragmentation. Addressing this challenge requires expanding the coverage of known biodiversity data, which remain to date restricted to a few well-surveyed regions. Here, we assess the sampling coverage and biases in biodiversity data on fruit-feeding butterflies at the Brazilian Atlantic Forest, discussing their effect on our understanding of the relationship between forest fragmentation and biodiversity at a large-scale. We hypothesize that sampling effort is biased towards large and connected fragments, which occur jointly in space at the Atlantic forest.

Methods

We used a comprehensive dataset of Atlantic Forest fruit-feeding butterfly communities to test for sampling biases towards specific geographical areas, climate conditions and landscape configurations.

Results

We found a pattern of geographical aggregation of sampling sites, independently of scale, and a strong sampling bias towards large and connected forest fragments, located near cities and roads. Sampling gaps are particularly acute in small and disconnected forest fragments and rare climate conditions. In contrast, currently available data can provide a fair picture of fruit-feeding butterfly communities in large and connected Atlantic Forest remnants.

Discussion

Biased data hamper the inference of the functional relationship between deforestation and biodiversity at a large-scale, since they are geographically clustered and have sampling gaps in small and disconnected fragments. These data are useful to inform decision-makers regarding conservation efforts to curb biodiversity loss in the Atlantic Forest. Thus, we suggest to expand sampling effort to small and disconnected forest fragments, which would allow more accurate evaluations of the effects of landscape modification.

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.

Reduction in the potential distribution of bumble bees (Apidae: Bombus) in Mesoamerica under different climate change scenarios: Conservation implications

Bumble bees are an ecologically and economically important group of pollinating insects worldwide. Global climate change is predicted to affect bumble bee ecology including habitat suitability and geographic distribution. Our study aims to estimate the impact of projected climate change on 18 Mesoamerican bumble bee species. We used ecological niche modeling (ENM) using current and future climate emissions scenarios (representative concentration pathway 4.5, 6.0, and 8.5) and models (CCSM4, HadGEM2-AO, and MIROC-ESM-CHEM). Regardless of the scenario and model applied, our results suggest that all bumble bee species are predicted to undergo a reduction in their potential distribution and habitat suitability due to projected climate change. ENMs based on low emission scenarios predict a distribution loss ranging from 7% to 67% depending on the species for the year 2050. Furthermore, we discovered that the reduction of bumble bee geographic range shape will be more evident at the margins of their distribution. The reduction of suitable habitat is predicted to be accompanied by a 100-500 m upslope change in altitude and 1-581 km shift away from the current geographic centroid of a species' distribution. On average, protected natural areas in Mesoamerica cover ~14% of each species' current potential distribution, and this proportion is predicted to increase to ~23% in the high emission climate change scenarios. Our models predict that climate change will reduce Mesoamerican bumble bee habitat suitability, especially for rare species, by reducing their potential distribution ranges and suitability. The small proportion of current and future potential distribution falling in protected natural areas suggests that such areas will likely have marginal contribution to bumble bee habitat conservation. Our results have the capacity to inform stakeholders in designing effective landscape management for bumble bees, which may include developing restoration plans for montane pine oak forests habitats and native flowering plants.

Type and intensity of surrounding human land use, not local environment, shape genetic structure of a native grassland plant

Landscape heterogeneity can shape genetic structure and functional connectivity of populations. When this heterogeneity imposes variable costs of moving across the landscape, populations can be structured according to a pattern of "isolation by resistance" (IBR). At the same time, divergent local environmental filters can limit gene flow, creating an alternative pattern of "isolation by environment" (IBE). Here, we evaluate IBR and IBE in the insect-pollinated, biennial plant Sabatia angularis (L.) Pursh (Gentianaceae) across serpentine grasslands in the fragmented landscape of SE Pennsylvania, USA using ~4500 neutral SNP loci. Specifically, we test the extent to which radical alteration of the landscape matrix by humans has fundamentally altered the cost of movement, imprinting a pattern of IBR dictated by land use type and intensity, and the potential for IBE in relation to a gradient of heavy metal concentrations found in serpentine soil. We reveal a strong signal of IBR and a weak signal of IBE across sites, indicating the greater importance of the landscape matrix in shaping genetic structure of S. angularis populations in the study region. Based on Circuitscape and least cost path approaches, we find that both low- and high-intensity urbanization resist gene flow by orders of magnitude greater than "natural" habitats, although resistance to low-intensity urbanization weakens at larger spatial scales. While cropland presents a substantially lower barrier than urban development, cumulative human land use surrounding populations predicts within-population genetic diversity and inbreeding in S. angularis. Our results emphasize the role of forest buffers and corridors in facilitating gene flow between serpentine grassland patches and averting local extinction of plant populations.

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.

The use and misuse of regression models in landscape genetic analyses

The field of landscape genetics has been rapidly evolving, adopting and adapting analytical frameworks to address research questions. Current studies are increasingly using regression-based frameworks to infer the individual contributions of landscape and habitat variables on genetic differentiation. This paper outlines appropriate and inappropriate uses of multiple regression for these purposes, and demonstrates through simulation the limitations of different analytical frameworks for making correct inference. Of particular concern are recent studies seeking to explain genetic differences by fitting regression models with effective distance variables calculated independently on separate landscape resistance surfaces. When moving across the landscape, organisms cannot respond independently and uniquely to habitat and landscape features. Analyses seeking to understand how landscape features affect gene flow should model a single conductance or resistance surface as a parameterized function of relevant spatial covariates, and estimate the values of these parameters by linking a single set of resistance distances to observed genetic dissimilarity via a loss function. While this loss function may involve a regression-like step, the associated nuisance parameters are not interpretable in terms of organismal movement and should not be conflated with what is actually of interest: the mapping between spatial covariates and conductance/resistance. The growth and evolution of landscape genetics as a field has been rapid and exciting. It is the goal of this paper to highlight past missteps and demonstrate limitations of current approaches to ensure that future use of regression models will appropriately consider the process being modeled, which will provide clarity to model interpretation.

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.

Patterns of genetic differentiation in Colorado potato beetle correlate with contemporary, not historic, potato land cover

Changing landscape heterogeneity can influence connectivity and alter genetic variation in local populations, but there can be a lag between ecological change and evolutionary responses. Temporal lag effects might be acute in agroecosystems, where land cover has changed substantially in the last two centuries. Here, we evaluate how patterns of an insect pest's genetic differentiation are related to past and present agricultural land cover change over a 150-year period. We quantified change in the amount of potato, Solanum tuberosum L., land cover since 1850 using county-level agricultural census reports, obtained allele frequency data from 7,408 single-nucleotide polymorphism loci, and compared effects of historic and contemporary landscape connectivity on genetic differentiation of Colorado potato beetle, Leptinotarsa decemlineata Say, in two agricultural landscapes in the United States. We found that potato land cover peaked in Wisconsin in the early 1900s, followed by rapid decline and spatial concentration, whereas it increased in amount and extent in the Columbia Basin of Oregon and Washington beginning in the 1960s. In both landscapes, we found small effect sizes of landscape resistance on genetic differentiation, but a 20× to 1,000× larger effect of contemporary relative to historic landscape resistances. Demographic analyses suggest population size trajectories were largely consistent among regions and therefore are not likely to have differentially impacted the observed patterns of population structure in each region. Weak landscape genetic associations might instead be related to the coarse resolution of our historical land cover data. Despite rapid changes in agricultural landscapes over the last two centuries, genetic differentiation among L. decemlineata populations appears to reflect ongoing landscape change. The historical landscape genetic framework employed in this study is broadly applicable to other agricultural pests and might reveal general responses of pests to agricultural land-use change.

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.

The effect of habitat fragmentation on the bee visitor assemblages of three Australian tropical rainforest tree species

Tropical forest loss and fragmentation can change bee community dynamics and potentially interrupt plant-pollinator relationships. While bee community responses to forest fragmentation have been investigated in a number of tropical regions, no studies have focused on this topic in Australia. In this study, we examine taxonomic and functional diversity of bees visiting flowers of three tree species across small and large rainforest fragments in Australian tropical landscapes. We found lower taxonomic diversity of bees visiting flowers of trees in small rainforest fragments compared with large forest fragments and show that bee species in small fragments were subsets of species in larger fragments. Bees visiting trees in small fragments also had higher mean body sizes than those in larger fragments, suggesting that small-sized bees may be less likely to persist in small fragments. Lastly, we found reductions in the abundance of eusocial stingless bees visiting flowers in small fragments compared to large fragments. These results suggest that pollinator visits to native trees living in small tropical forest remnants may be reduced, which may in turn impact on a range of processes, potentially including forest regeneration and diversity maintenance in small forest remnants in Australian tropical countryside landscapes.

Conserving relics from ancient underground worlds: assessing the influence of cave and landscape features on obligate iron cave dwellers from the Eastern Amazon

The degradation of subterranean habitats is believed to represent a serious threat for the conservation of obligate subterranean dwellers (troglobites), many of which are short-range endemics. However, while the factors influencing cave biodiversity remain largely unknown, the influence of the surrounding landscape and patterns of subterranean connectivity of terrestrial troglobitic communities have never been systematically assessed. Using spatial statistics to analyze the most comprehensive speleological database yet available for tropical caves, we first assess the influence of iron cave characteristics and the surrounding landscape on troglobitic communities from the Eastern Amazon. We then determine the spatial pattern of troglobitic community composition, species richness, phylogenetic diversity, and the occurrence of frequent troglobitic species, and finally quantify how different landscape features influence the connectivity between caves. Our results reveal the key importance of habitat amount, guano, water, lithology, geomorphology, and elevation in shaping iron cave troglobitic communities. While mining within 250 m from the caves influenced species composition, increasing agricultural land cover within 50 m from the caves reduced species richness and phylogenetic diversity. Troglobitic species composition, species richness, phylogenetic diversity, and the occurrence of frequent troglobites showed spatial autocorrelation for up to 40 km. Finally, our results suggest that the conservation of cave clusters should be prioritized, as geographic distance was the main factor determining connectivity between troglobitic communities. Overall, our work sheds important light onto one of the most overlooked terrestrial ecosystems, and highlights the need to shift conservation efforts from individual caves to subterranean habitats as a whole.

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.

Adding landscape genetics and individual traits to the ecosystem function paradigm reveals the importance of species functional breadth

Animal pollination mediates both reproduction and gene flow for the majority of plant species across the globe. However, past functional studies have focused largely on seed production; although useful, this focus on seed set does not provide information regarding species-specific contributions to pollen-mediated gene flow. Here we quantify pollen dispersal for individual pollinator species across more than 690 ha of tropical forest. Specifically, we examine visitation, seed production, and pollen-dispersal ability for the entire pollinator community of a common tropical tree using a series of individual-based pollinator-exclusion experiments followed by molecular-based fractional paternity analyses. We investigate the effects of pollinator body size, plant size (as a proxy of floral display), local plant density, and local plant kinship on seed production and pollen-dispersal distance. Our results show that while large-bodied pollinators set more seeds per visit, small-bodied bees visited flowers more frequently and were responsible for more than 49% of all long-distance (beyond 1 km) pollen-dispersal events. Thus, despite their size, small-bodied bees play a critical role in facilitating long-distance pollen-mediated gene flow. We also found that both plant size and local plant kinship negatively impact pollen dispersal and seed production. By incorporating genetic and trait-based data into the quantification of pollination services, we highlight the diversity in ecological function mediated by pollinators, the influential role that plant and population attributes play in driving service provision, and the unexpected importance of small-bodied pollinators in the recruitment of plant genetic diversity.

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population-specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.

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

Across the globe, wild bees are threatened by ongoing natural habitat loss, risking the maintenance of plant biodiversity and agricultural production. Despite the ecological and economic importance of wild bees and the fact that several species are now managed for pollination services worldwide, little is known about how land use and beekeeping practices jointly influence gene flow. Using stingless bees as a model system, containing wild and managed species that are presumed to be particularly susceptible to habitat degradation, here we examine the main drivers of tropical bee gene flow. We employ a novel landscape genetic approach to analyse data from 135 populations of 17 stingless bee species distributed across diverse tropical biomes within the Americas. Our work has important methodological implications, as we illustrate how a maximum-likelihood approach can be applied in a meta-analysis framework to account for multiple factors, and weight estimates by sample size. In contrast to previously held beliefs, gene flow was not related to body size or deforestation, and isolation by geographic distance (IBD) was significantly affected by management, with managed species exhibiting a weaker IBD than wild ones. Our study thus reveals the critical importance of beekeeping practices in shaping the patterns of genetic differentiation across bee species. Additionally, our results show that many stingless bee species maintain high gene flow across heterogeneous landscapes. We suggest that future efforts to preserve wild tropical bees should focus on regulating beekeeping practices to maintain natural gene flow and enhancing pollinator-friendly habitats, prioritizing species showing a limited dispersal ability.