Spatio-temporal shifts driven by climate change threaten persistence and resilience of honey bee populations

This study utilizes genomic data and climate modeling to demonstrate that climate change will disrupt the spatial genetic structure of honey bee populations in Anatolia and Thrace, leading to reduced persistence and resilience for most lineages while expanding those adapted to warmer, drier conditions, thereby rendering current conservation areas ineffective.

The Gist

Imagine the honey bee (Apis mellifera) not just as a single insect, but as a massive, global family with five distinct "cousin" groups living in Turkey. Each cousin has developed a unique personality and set of survival skills over thousands of years to match their specific neighborhood: some are experts at surviving cold, wet mountains, while others are masters of hot, dry deserts.

This paper is like a genetic weather forecast. The author, Mert Kükrer, asks a scary question: If the climate keeps changing, will these unique cousins survive, or will they all get mixed up into one generic, confused soup?

Here is the story of the research, broken down into simple concepts:

1. The Map of the Family Tree

First, the researcher took DNA samples from hundreds of bees across Turkey (a place where five different bee subspecies meet). Think of this as taking a family photo of every cousin in the neighborhood.

• The Finding: He confirmed there are indeed five distinct groups. However, they aren't living in neat, separate boxes. They are mixing together in the middle, like colors blending on a painter's palette.
• The Analogy: Imagine a neighborhood where the "Mountain People" live in the north, the "Desert People" live in the south, and in the middle, you have families who are a mix of both. The study mapped exactly where these mixtures happen.

2. The Climate "Recipe" for Survival

The study used powerful computer models (called Gradient Forests and Generalized Dissimilarity Modeling) to figure out why the bees look different in different places.

• The Finding: It's not just about how far apart the bees live; it's about the weather.
  • The Caucasian bees (from the mountains) need cool, moist air. They are like a delicate houseplant that needs shade and water.
  • The Levantine bees (from the hot south) love heat and dryness. They are like a cactus that thrives in the sun.
• The Analogy: Think of each bee group as a specific recipe. The "Mountain Recipe" requires cool temperatures and rain. The "Desert Recipe" requires heat and dryness. If you try to bake the Mountain Cake in a desert oven, it will burn.

3. The Future Forecast: A Storm is Coming

The researcher then asked the computer: "What happens if we run these recipes under the climate scenarios predicted for the year 2100?"

• The Bad News: The "Mountain Recipes" are in trouble. As the world gets hotter and drier, the cool, moist zones where the mountain bees live are shrinking. These bees are losing their homes.
• The "Invasive" Winner: The "Desert Recipes" (the hot-climate bees) are expanding. They are moving into new areas, taking over territory that used to belong to the cooler-climate bees.
• The Analogy: Imagine a game of musical chairs where the chairs (habitats) are disappearing. The "Mountain Bees" are being forced out of their seats. The "Desert Bees" are running faster, taking over the empty chairs, and pushing the mountain bees into smaller and smaller corners.

4. The "Genetic Melting Pot" Problem

As the climate changes, the distinct boundaries between these bee groups are dissolving.

• The Risk: When the cool-climate bees are forced into the hot-climate zones, they have to mate with the local hot-climate bees to survive. This creates a "genetic soup."
• The Consequence: We might lose the unique "Mountain Bee" entirely. It won't just move; its unique DNA will get diluted until it's gone. This is like losing a rare spice in a soup because you kept adding water and common salt. The soup still tastes like soup, but it's no longer the special dish it used to be.

5. Are Our Safety Nets Working?

Turkey has set up "Genetic Conservation Areas" (GCAs)—protected zones where beekeepers are forbidden from moving bees around, to keep the local breeds pure.

• The Finding: The study checked these protected zones and found they are not enough.
  • Some unique bee groups (like the Zagrosian bee) aren't protected at all.
  • Even the protected zones are in areas that will become too hot or dry for the bees living there.
• The Analogy: It's like building a lifeboat for a ship that is sinking, but the lifeboat is made of wood that will rot in the rising water. The current safety zones are in the wrong places for the future climate.

The Bottom Line

Climate change isn't just about bees dying; it's about bees changing. The unique, locally adapted bees that have survived for thousands of years are being pushed out by a warming world.

The study suggests we need to:

1. Move the Safety Nets: Create new protected areas in places that will remain cool and suitable in the future (like higher mountains).
2. Stop the Mixing: Be careful about moving bees around, as it speeds up the loss of unique genetic traits.
3. Watch the Weather: Understand that the bees we have today might not be the bees we have in 50 years.

In short: The climate is changing the menu, and unless we protect the specific ingredients (the unique bee breeds) in the right kitchens (the right habitats), we risk losing the special flavors of nature forever.

Technical Summary

1. Problem Statement

The study addresses a critical gap in understanding how climate change impacts intraspecific genetic diversity, specifically in the western honey bee (Apis mellifera), a species vital for global pollination and agriculture. While the effects of climate change on species distribution ranges are well-documented, there is limited research on how ancestry composition and genetic turnover within a species respond to environmental gradients, particularly when confounded by human-mediated dispersal (e.g., migratory beekeeping and trade).

The specific challenges identified include:

• The lack of predictive models for how climate change will alter the spatial genetic structure of honey bee subspecies in historical refugia like Anatolia and Thrace.
• The uncertainty regarding whether existing Genetic Conservation Areas (GCAs) in Turkey will remain effective under future climate scenarios.
• The need to distinguish between natural adaptation and anthropogenic gene flow in shaping population dynamics.

2. Methodology

The research employs a multi-step framework integrating population genetics, machine learning, and ecological forecasting:

• Data Collection & Genotyping:

  • Samples: 672 individuals genotyped at 30 microsatellite loci.
  • Scope: 392 localities across 75 provinces in Turkey (covering five native subspecies: A. m. syriaca, A. m. caucasica, A. m. anatoliaca, A. m. meda, and the Thracian ecotype), plus reference samples from Austria, Germany, and Georgia.
  • Filtering: Samples were filtered to remove obvious mismatches and siblings, resulting in 592 samples for downstream analysis.

• Population Structure Analysis:

  • Used Structure 2.3.4 to estimate individual membership coefficients, identifying five distinct ancestral groups.
  • Applied UPGMA phylogenies, PCA, and AMOVA to characterize genetic differentiation.

• Modeling Intra-specific Turnover:

  • Gradient Forests (GFs): A machine learning approach used to model compositional turnover as a function of 19 bioclimatic variables (WorldClim 2.1) and 18 environmental variables (Envirem), plus spatial predictors (Moran's Eigenvector Maps - MEMs). This identified key environmental drivers of genetic variation.
  • Generalized Dissimilarity Modelling (GDM): Used to relate biological distance (ancestry composition) to environmental and geographic distances. This modeled the rate of turnover across the landscape.

• Spatio-Temporal Forecasting:

  • Climate Scenarios: Projections were made using six CMIP6 Global Circulation Models under four Shared Socioeconomic Pathways (SSPs: 126, 245, 370, 585) for the years 2030, 2050, 2070, and 2090.
  • Vulnerability Indices: Four novel indices were calculated to assess climate impact:
    1. Persistence: Inverse of the offset between current and future ecological similarity to a site.
    2. Resilience: Ratio of future-to-current similarity to neighbors (identifying refugia).
    3. Disappearance: Loss of current ancestry composition.
    4. Emergence: Appearance of novel ancestry compositions.
  • Conservation Assessment: Evaluated the "resemblance" of the landscape to existing GCAs and simulated the addition of new sites (Hakkari, Çankırı, Muş) to improve representativeness.

3. Key Contributions

• Novel Modeling Framework: This is the first study (to the authors' knowledge) to use ancestry composition as a multidimensional proxy for genetic variation in GDM and GF models to forecast intraspecific biodiversity dynamics under climate change.
• Integration of Human and Natural Factors: The study successfully disentangles the effects of natural environmental gradients from human-mediated gene flow in a complex landscape.
• Quantitative Conservation Planning: It provides a scalable, data-driven method to evaluate the future efficacy of protected areas (GCAs) and identify gaps in genetic representation.
• High-Resolution Forecasting: The study moves beyond species-level range shifts to predict specific changes in subspecies distribution and admixture zones at a fine spatial scale (3km buffer).

4. Key Results

• Population Structure: Five major ancestral groups were identified (Thracian, Anatolian, Caucasian, Levantine, Zagrosian). While geographic distance explains significant variation, climatic variables account for a large fraction of ancestry turnover.
• Environmental Drivers:
  • Spatial processes (MEMs) were the strongest predictors, followed by climatic variables.
  • Key drivers included precipitation during the driest period ($P_{driest}$), minimum temperature of the warmest period ($minT_{warm}$), and potential evapotranspiration.
  • Caucasian ancestry was highly sensitive to moisture-related variables, while Levantine ancestry responded strongly to temperature variables.
• Climate Change Impacts (21st Century Projections):
  • Shift in Dominance: Ancestral groups associated with warmer/drier conditions (specifically the Levantine group) are projected to expand significantly, nearly doubling their range by 2090 under high-emission scenarios (SSP585).
  • Decline of Native Groups: The Thracian, Caucasian, and Anatolian groups face substantial range contractions. The Thracian group is particularly threatened, potentially disappearing from within Turkey's borders in some scenarios.
  • Loss of Uniqueness: Regions of high genetic uniqueness and refugia are contracting. The area classified as "highly unique" decreased from ~100,000 km² to ~75,000 km².
  • Increased Turnover: The proportion of the landscape experiencing rapid genetic turnover (fast turnover sites) increased from ~100,000 km² to ~150,000 km², indicating increased hybridization and genetic swamping.
• Conservation Efficacy:
  • Existing GCAs provide incomplete representation, particularly for the Zagrosian group.
  • Adding hypothetical GCAs in Hakkari, Çankırı, and Muş significantly improved the resemblance index and coverage, nearly covering the entire surface area of Turkey with high-resemblance zones.
  • Maladaptation Risk: Many current GCAs are located in areas with low persistence indices, suggesting that locally adapted stocks within these protected areas may become maladapted to future climates.

5. Significance

• Conservation Urgency: The findings highlight that current conservation strategies are insufficient to protect the unique genetic diversity of Anatolian honey bees against climate-driven homogenization. The study calls for the immediate establishment of new GCAs in high-persistence and high-resilience zones (e.g., mountainous refugia).
• Biodiversity Monitoring: The proposed framework serves as a blueprint for monitoring intraspecific genetic diversity in other managed and wild species facing climate change.
• Economic and Ecological Impact: By predicting the loss of locally adapted traits (e.g., disease resistance, honey yield, swarming behavior), the study underscores the economic risks to beekeeping and the ecological risks to pollination services.
• Policy Relevance: The results provide a scientific basis for updating national genetic monitoring programs and breeding strategies to ensure the resilience of honey bee populations in the face of global change.

In conclusion, the paper demonstrates that climate change will likely disrupt the spatial genetic structure of Apis mellifera, promoting admixture and threatening the persistence of distinct subspecies. It advocates for a shift from static conservation to dynamic, climate-informed management strategies.