QUALIS: The journal ranking system undermining the impact of Brazilian science.

A system called QUALIS was implemented in Brazil in 2009, intended to rank graduate programs from different subject areas and promote selected national journals. Since this system uses a complicated suit of criteria (differing among subject areas) to group journals into discrete categories, it could potentially create incentives to publish in low-impact journals ranked highly by QUALIS. Here I assess the influence of the QUALIS journal ranking system on the global impact of Brazilian science. Brazil shows a steeper decrease in the number of citations per document since the implementation of this QUALIS system, compared to the top Latin American countries publishing more scientific articles. All subject areas showed some degree of bias, with social sciences being usually more biased than natural sciences. Lastly, the decrease in the number of citations over time proved steeper in a more biased subject area, suggesting a faster shift towards low-impact journals. Overall, the findings documented here suggest that the QUALIS system has undermined the global impact of Brazilian science, and reinforce a recent recommendation from an official committee evaluating graduate programs to eliminate QUALIS. A system based on impact metrics could avoid introducing distorted incentives, and thereby boost the global impact of Brazilian science.

Integrating environmental variables by multivariate ordination enables the reliable estimation of mineland rehabilitation status.

Despite the wide variety of variables commonly employed to measure the success of rehabilitation, the assessment and subsequent definition of indicators of environmental rehabilitation status are not simple tasks. The main challenges are comparing rehabilitated sites with target ecosystems as well as integrating individual environmental and eventually collinear variables into a single tractable measure for the state of a system before effective indicators that track rehabilitation may be modeled. Furthermore, a consensus is lacking regarding which and how many variables need to be surveyed for a reliable estimation of rehabilitation status. Here, we propose a multivariate ordination to integrate variables related to ecological processes, vegetation structure, and community diversity into a single estimation of rehabilitation status. As a case, we employed a curated set of 32 environmental variables retrieved from nonrevegetated, rehabilitating and reference sites associated with iron ore mines from the Urucum Massif, Mato Grosso do Sul, Brazil. By integrating this set of environmental variables into a single estimation of rehabilitation status, the proposed multivariate approach is straightforward and able to adequately address collinearity among variables. The proposed methodology allows for the identification of biases towards single variables, surveys or analyses, which is necessary to rank environmental variables regarding their importance to the assessment. Furthermore, we show that bootstrapping permitted the detection of the minimum number of environmental variables necessary to achieve reliable estimations of the rehabilitation status. Finally, we show that the proposed variable integration enables the definition of case-specific environmental indicators for more rapid assessments of mineland rehabilitation. Thus, the proposed multivariate ordination represents a powerful tool to facilitate the diagnosis of rehabilitating sites worldwide provided that sufficient environmental variables related to ecological processes, diversity and vegetation structure are gathered from nonrehabilitated, rehabilitating and reference study sites. By identifying deviations from predicted rehabilitation trajectories and providing assessments for environmental agencies, this proposed multivariate ordination increases the effectiveness of (mineland) rehabilitation.

Competitive males have higher quality sperm in a monogamous social bee.

BACKGROUND: Reproductive success is determined by the interplay between mating and fertilization success. In social insect species with male-biased sex ratios and queen monogamy, males face particularly strong pre-copulatory sexual selection since they must compete with thousands of other males for a unique mating opportunity. Ejaculate quality is also expected to be under selection, because queens are long-lived and store sperm for life, so males with higher quality ejaculates are expected to provide queens with larger and longer-lived colonies, which in turn may produce more daughter queens (the only direct fitness gains of haplodiploid males). Considering the action of pre and post-copulatory sexual selection on male traits, three scenarios might thus be expected: positive, negative or no association between male mating ability and fertilization success. Here we explored these scenarios in the stingless bee Scaptotrigona aff. depilis, where males gather in large aggregations and queens mate with a single male. Male mating ability was assessed through the capacity of a male to reach an aggregation and persist on it; while sperm viability, sperm number, and sperm morphology were used as proxies for sperm quality. RESULTS: Sperm viability was associated with persistence time in the aggregation, and males that persisted longer presented shorter spermatozoa and higher variation in sperm length than recently arrived males. However, sperm traits of males that reached aggregations did not differ from those of males collected inside their colonies. In addition, males that persisted longer in aggregations were smaller than other males. Male size and sperm viability were not correlated, suggesting that the observed patterns were not due to trade-offs in male resource allocation. CONCLUSIONS: Persistence in male aggregations thus seems to select more competitive males with higher quality sperm. Our work is the first one to reveal an association between male competitive ability and fertilization success in a monogamous social insect. This finding sheds important light on the evolution of male traits in social insects and the general mechanisms of sexual selection.

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.

Genetic variability in captive populations of the stingless bee Tetragonisca angustula.

Low genetic variability has normally been considered a consequence of animal husbandry and a major contributing factor to declining bee populations. Here, we performed a molecular analysis of captive and wild populations of the stingless bee Tetragonisca angustula, one of the most commonly kept species across South America. Microsatellite analyses showed similar genetic variability between wild and captive populations However, captive populations showed lower mitochondrial genetic variability. Male-mediated gene flow, transport and division of nests are suggested as the most probable explanations for the observed patterns of genetic structure. We conclude that increasing the number of colonies kept through nest divisions does not negatively affect nuclear genetic variability, which seems to be maintained by small-scale male dispersal and human-mediated nest transport. However, the transport of nests from distant localities should be practiced with caution given the high genetic differentiation observed between samples from western and eastern areas. The high genetic structure verified is the result of a long-term evolutionary process, and bees from distant localities may represent unique evolutionary lineages.

Temporal Variation in Honey Production by the Stingless Bee Melipona subnitida (Hymenoptera: Apidae): Long-Term Management Reveals its Potential as a Commercial Species in Northeastern Brazil.

Even though stingless beekeeping has a great potential as a sustainable development tool, the activity remains essentially informal, technical knowledge is scarce, and management practices lack the sophistication and standardization found in apiculture. Here, we contributed to the further development of stingless beekeeping by investigating the long-term impact of management and climate on honey production and colony survival in the stingless bee Melipona subnitida Ducke (1910). We analyzed a 10-yr record of 155 M. subnitida colonies kept by a commercial honey producer of northeastern Brazil. This constitutes the longest and most accurate record available for a stingless bee. We modeled honey production in relation to time (years), age, management practices (colony division and food supplementation), and climatic factors (temperature and precipitation), and used a model selection approach to identify which factors best explained honey production. We also modeled colony mortality in relation to climatic factors. Although the amount of honey produced by each colony decreased over time, we found that the probability of producing honey increased over the years. Colony divisions decreased honey production, but did not affect honey presence, while supplementary feeding positively affected honey production. In warmer years, the probability of producing honey decreased and the amount of honey produced was lower. In years with lower precipitation, fewer colonies produced honey. In contrast, colony mortality was not affected by climatic factors, and some colonies lived up to nine years, enduring extreme climatic conditions. Our findings provide useful guidelines to improve management and honey production in stingless bees.

Monogamy in large bee societies: a stingless paradox.

High genetic diversity is important for the functioning of large insect societies. Across the social Hymenoptera (ants, bees, and wasps), species with the largest colonies tend to have a high colony-level genetic diversity resulting from multiple queens (polygyny) or queens that mate with multiple males (polyandry). Here we studied the genetic structure of Trigona spinipes, a stingless bee species with colonies an order of magnitude larger than those of polyandrous honeybees. Genotypes of adult workers and pupae from 43 nests distributed across three Brazilian biomes showed that T. spinipes colonies are usually headed by one singly mated queen. Apart from revealing a notable exception from the general incidence of high genetic diversity in large insect societies, our results reinforce previous findings suggesting the absence of polyandry in stingless bees and provide evidence against the sperm limitation hypothesis for the evolution of polyandry. Stingless bee species with large colonies, such as T. spinipes, thus seem promising study models to unravel alternative mechanisms to increase genetic diversity within colonies or understand the adaptive value of low genetic diversity in large insect societies.

Conserving genetic diversity in the honeybee: comments on Harpur et al. (2012).

The article by Harpur et al. (2012) 'Management increases genetic diversity of honey bees via admixture' concludes that '…honey bees do not suffer from reduced genetic diversity caused by management and, consequently, that reduced genetic diversity is probably not contributing to declines of managed Apis mellifera populations'. In the light of current honeybee and beekeeping declines and their consequences for honeybee conservation and the pollination services they provide, we would like to express our concern about the conclusions drawn from the results of Harpur et al. (2012). While many honeybee management practices do not imply admixture, we are convinced that the large-scale genetic homogenization of admixed populations could drive the loss of valuable local adaptations. We also point out that the authors did not account for the extensive gene flow that occurs between managed and wild/feral honeybee populations and raise concerns about the data set used. Finally, we caution against underestimating the importance of genetic diversity for honeybee colonies and highlight the importance of promoting the use of endemic honeybee subspecies in apiculture.

Global genetic diversity status and trends: towards a suite of Essential Biodiversity Variables (EBVs) for genetic composition.

Biodiversity underlies ecosystem resilience, ecosystem function, sustainable economies, and human well-being. Understanding how biodiversity sustains ecosystems under anthropogenic stressors and global environmental change will require new ways of deriving and applying biodiversity data. A major challenge is that biodiversity data and knowledge are scattered, biased, collected with numerous methods, and stored in inconsistent ways. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has developed the Essential Biodiversity Variables (EBVs) as fundamental metrics to help aggregate, harmonize, and interpret biodiversity observation data from diverse sources. Mapping and analyzing EBVs can help to evaluate how aspects of biodiversity are distributed geographically and how they change over time. EBVs are also intended to serve as inputs and validation to forecast the status and trends of biodiversity, and to support policy and decision making. Here, we assess the feasibility of implementing Genetic Composition EBVs (Genetic EBVs), which are metrics of within-species genetic variation. We review and bring together numerous areas of the field of genetics and evaluate how each contributes to global and regional genetic biodiversity monitoring with respect to theory, sampling logistics, metadata, archiving, data aggregation, modeling, and technological advances. We propose four Genetic EBVs: (i) Genetic Diversity; (ii) Genetic Differentiation; (iii) Inbreeding; and (iv) Effective Population Size (Ne ). We rank Genetic EBVs according to their relevance, sensitivity to change, generalizability, scalability, feasibility and data availability. We outline the workflow for generating genetic data underlying the Genetic EBVs, and review advances and needs in archiving genetic composition data and metadata. We discuss how Genetic EBVs can be operationalized by visualizing EBVs in space and time across species and by forecasting Genetic EBVs beyond current observations using various modeling approaches. Our review then explores challenges of aggregation, standardization, and costs of operationalizing the Genetic EBVs, as well as future directions and opportunities to maximize their uptake globally in research and policy. The collection, annotation, and availability of genetic data has made major advances in the past decade, each of which contributes to the practical and standardized framework for large-scale genetic observation reporting. Rapid advances in DNA sequencing technology present new opportunities, but also challenges for operationalizing Genetic EBVs for biodiversity monitoring regionally and globally. With these advances, genetic composition monitoring is starting to be integrated into global conservation policy, which can help support the foundation of all biodiversity and species' long-term persistence in the face of environmental change. We conclude with a summary of concrete steps for researchers and policy makers for advancing operationalization of Genetic EBVs. The technical and analytical foundations of Genetic EBVs are well developed, and conservation practitioners should anticipate their increasing application as efforts emerge to scale up genetic biodiversity monitoring regionally and globally.

Optimizing speleological monitoring efforts: insights from long-term data for tropical iron caves.

Understanding the factors underpinning species abundance patterns in space and time is essential to implement effective cave conservation actions. Yet, the methods employed to monitor cave biodiversity still lack standardization, and no quantitative assessment has yet tried to optimize the amount and type of information required to efficiently identify disturbances in cave ecosystems. Using a comprehensive monitoring dataset for tropical iron caves, comprising abundance measurements for 33 target taxa surveyed across 95 caves along four years, here we provide the first evidence-based recommendations to optimize monitoring programs seeking to follow target species abundance through time. We found that seasonality did not influence the ability to detect temporal abundance trends. However, in most species, abundance estimates assessed during the dry season resulted in a more accurate detection of temporal abundance trends, and at least three surveys were required to identify global temporal abundance trends. Finally, we identified a subset of species that could potentially serve as short-term disturbance indicators. Results suggest that iron cave monitoring programs implemented in our study region could focus sampling efforts in the dry season, where detectability of target species is higher, while assuring data collection for at least three years. More generally, our study reveals the importance of long-term cave monitoring programs for detecting possible disturbances in subterranean ecosystems, and for using the generated information to optimize future monitoring efforts.

Combining genotype, phenotype, and environmental data to delineate site-adjusted provenance strategies for ecological restoration.

Despite the importance of climate-adjusted provenancing to mitigate the effects of environmental change, climatic considerations alone are insufficient when restoring highly degraded sites. Here we propose a comprehensive landscape genomic approach to assist the restoration of moderately disturbed and highly degraded sites. To illustrate it we employ genomic data sets comprising thousands of single nucleotide polymorphisms from two plant species suitable for the restoration of iron-rich Amazonian Savannas. We first use a subset of neutral loci to assess genetic structure and determine the genetic neighbourhood size. We then identify genotype-phenotype-environment associations, map adaptive genetic variation, and predict adaptive genotypes for restoration sites. Whereas local provenances were found optimal to restore a moderately disturbed site, a mixture of genotypes seemed the most promising strategy to recover a highly degraded mining site. We discuss how our results can help define site-adjusted provenancing strategies, and argue that our methods can be more broadly applied to assist other restoration initiatives.

Conservation implications of genetic structure in the narrowest endemic quillwort from the Eastern Amazon.

The quillwort Isoëtes cangae is a critically endangered species occurring in a single lake in Serra dos Carajás, Eastern Amazon. Low genetic diversity and small effective population sizes (N e) are expected for narrow endemic species (NES). Conservation biology studies centered in a single species show some limitations, but they are still useful considering the limited time and resources available for protection of species at risk of extinction. Here, we evaluated the genetic diversity, population structure, N e, and minimum viable population (MVP) of I. cangae to provide information for effective conservation programs. Our analyses were based on 55 individuals collected from the Amendoim Lake and 35,638 neutral SNPs. Our results indicated a single panmictic population, moderate levels of genetic diversity, and N e in the order of thousands, contrasting the expected for NES. Negative FIS values were also found, suggesting that I. cangae is not under risk of inbreeding depression. Our findings imply that I. cangae contains enough genetic diversity to ensure evolutionary potential and that all individuals should be treated as one demographic unit. These results provide essential information to optimize ex situ conservation efforts and genetic diversity monitoring, which are currently applied to guide I. cangae conservation plans.

Mating flights select for symmetry in honeybee drones (Apis mellifera).

Males of the honeybee (Apis mellifera) fly to specific drone congregation areas (DCAs), which virgin queens visit in order to mate. From the thousands of drones that are reared in a single colony, only very few succeed in copulating with a queen, and therefore, a strong selection is expected to act on adult drones during their mating flights. In consequence, the gathering of drones at DCAs may serve as an indirect mate selection mechanism, assuring that queens only mate with those individuals having a better flight ability and a higher responsiveness to the queen's visual and chemical cues. Here, we tested this idea relying on wing fluctuating asymmetry (FA) as a measure of phenotypic quality. By recapturing marked drones at a natural DCA and comparing their size and FA with a control sample of drones collected at their maternal hives, we were able to detect any selection on wing size and wing FA occurring during the mating flights. Although we found no solid evidence for selection on wing size, wing FA was found to be significantly lower in the drones collected at the DCA than in those collected at the hives. Our results demonstrate the action of selection during drone mating flights for the first time, showing that developmental stability can influence the mating ability of honeybee drones. We therefore conclude that selection during honeybee drone mating flights may confer some fitness advantages to the queens.

Estimating the density of honeybee colonies across their natural range to fill the gap in pollinator decline censuses.

Although pollinator declines are a global biodiversity threat, the demography of the western honeybee (Apis mellifera) has not been considered by conservationists because it is biased by the activity of beekeepers. To fill this gap in pollinator decline censuses and to provide a broad picture of the current status of honeybees across their natural range, we used microsatellite genetic markers to estimate colony densities and genetic diversity at different locations in Europe, Africa, and central Asia that had different patterns of land use. Genetic diversity and colony densities were highest in South Africa and lowest in Northern Europe and were correlated with mean annual temperature. Confounding factors not related to climate, however, are also likely to influence genetic diversity and colony densities in honeybee populations. Land use showed a significantly negative influence over genetic diversity and the density of honeybee colonies over all sampling locations. In Europe honeybees sampled in nature reserves had genetic diversity and colony densities similar to those sampled in agricultural landscapes, which suggests that the former are not wild but may have come from managed hives. Other results also support this idea: putative wild bees were rare in our European samples, and the mean estimated density of honeybee colonies on the continent closely resembled the reported mean number of managed hives. Current densities of European honeybee populations are in the same range as those found in the adverse climatic conditions of the Kalahari and Saharan deserts, which suggests that beekeeping activities do not compensate for the loss of wild colonies. Our findings highlight the importance of reconsidering the conservation status of honeybees in Europe and of regarding beekeeping not only as a profitable business for producing honey, but also as an essential component of biodiversity conservation.

Range-wide neutral and adaptive genetic structure of an endemic herb from Amazonian Savannas.

Conserving genetic diversity in rare and narrowly distributed endemic species is essential to maintain their evolutionary potential and minimize extinction risk under future environmental change. In this study we assess neutral and adaptive genetic structure and genetic diversity in Brasilianthus carajensis (Melastomataceae), an endemic herb from Amazonian Savannas. Using RAD sequencing we identified a total of 9365 SNPs in 150 individuals collected across the species' entire distribution range. Relying on assumption-free genetic clustering methods and environmental association tests we then compared neutral with adaptive genetic structure. We found three neutral and six adaptive genetic clusters, which could be considered management units (MU) and adaptive units (AU), respectively. Pairwise genetic differentiation (F ST) ranged between 0.024 and 0.048, and even though effective population sizes were below 100, no significant inbreeding was found in any inferred cluster. Nearly 10 % of all analysed sequences contained loci associated with temperature and precipitation, from which only 25 sequences contained annotated proteins, with some of them being very relevant for physiological processes in plants. Our findings provide a detailed insight into genetic diversity, neutral and adaptive genetic structure in a rare endemic herb, which can help guide conservation and management actions to avoid the loss of unique genetic variation.

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.

Mapping and quantification of ferruginous outcrop savannas in the Brazilian Amazon: A challenge for biodiversity conservation.

The eastern Brazilian Amazon contains many isolated ferruginous savanna ecosystem patches (locally known as 'canga vegetation') located on ironstone rocky outcrops on the top of plateaus and ridges, surrounded by tropical rainforests. In the Carajás Mineral Province (CMP), these outcrops contain large iron ore reserves that have been exploited by opencast mining since the 1980s. The canga vegetation is particularly impacted by mining, since the iron ores that occur are associated with this type of vegetation and currently, little is known regarding the extent of canga vegetation patches before mining activities began. This information is important for quantifying the impact of mining, in addition to helping plan conservation programmes. Here, land cover changes of the Canga area in the CMP are evaluated by estimating the pre-mining area of canga patches and comparing it to the actual extent of canga patches. We mapped canga vegetation using geographic object-based image analysis (GEOBIA) from 1973 Landsat-1 MSS, 1984 and 2001 Landsat-5 TM, and 2016 Landsat-8 OLI images, and found that canga vegetation originally occupied an area of 144.2 km2 before mining exploitation. By 2016, 19.6% of the canga area was lost in the CMP due to conversion to other land-use types (mining areas, pasturelands). In the Carajás National Forest (CNF), located within the CMP, the original canga vegetation covered 105.2 km2 (2.55% of the CNF total area), and in 2016, canga vegetation occupied an area of 77.2 km2 (1.87%). Therefore, after more than three decades of mineral exploitation, less than 20% of the total canga area was lost. Currently, 21% of the canga area in the CMP is protected by the Campos Ferruginosos National Park. By documenting the initial extent of canga vegetation in the eastern Amazon and the extent to which it has been lost due to mining operations, the results of this work are the first step towards conserving this ecosystem.

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.

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.