The research, led by first author Ahmet Kılınç of the General Surgery Department at Elmalı State Hospital in Antalya, with corresponding author Mehmet Eşref Ulutaş and colleagues including Çiğdem Damla Deniz, Haydar Öztaş, İsmail Hasırcı, and Adil Kartal, was designed to test honey in two distinct modes: prophylaxis and therapy. Rather than studying established, late-stage tumors, the investigators focused on the very beginning of the disease process, when pancreatic tissue first begins to drift toward malignancy. This choice matters because pancreatic cancer is usually diagnosed late, and any intervention that could act at the precursor stage would be of enormous clinical interest if it could ever be translated to humans.

To model the disease, the team used the azaserine-rat model, a well-established experimental system in pancreatic cancer research. Azaserine is a chemical carcinogen that, when administered to rats, reliably induces lesions in the pancreas that mirror the early steps of human pancreatic tumor formation. The key structures the researchers tracked are called atypical acinar cell foci, or AACF, which are microscopic clusters of abnormal acinar cells regarded as precursor lesions of pancreatic cancer. Because these foci appear early, are quantifiable under the microscope, and correlate with later tumor burden, they serve as sensitive surrogate endpoints for chemoprevention studies, allowing researchers to measure an intervention’s effect long before overt tumors would develop.

The experimental design divided the animals into five groups: a control group, a honey-only group, an azaserine group that received the carcinogen alone, a therapeutic group, and a prophylactic group. The distinction between the last two is central to the study’s message. In the prophylactic group, honey was administered at the same time as azaserine, testing whether the substance could protect pancreatic tissue while it was being assaulted by the carcinogen. In the therapeutic group, honey was started three weeks after azaserine administration, a delay designed to simulate a more clinically realistic scenario in which the disease process is already underway before treatment begins. After eleven weeks of follow-up, the animals were sacrificed and their pancreatic tissues removed for detailed histopathological examination.

The quantitative readout of the study was rigorous. The researchers did not simply report whether lesions were present or absent; they measured three separate morphometric parameters of the atypical acinar cell foci. The first was the number of focus intersections per unit area, essentially a density measure of how many precursor lesions were packed into a given region of pancreatic tissue. The second was the ratio of the total AACF size to the size of the whole pancreas, capturing how much of the organ had been overtaken by abnormal cells. The third was the estimated number of foci per unit volume, a stereological estimate that corrects for the fact that a two-dimensional tissue section underestimates three-dimensional lesion burden. Together, these metrics provide a multidimensional picture of how extensively the carcinogen had transformed the pancreas in each group.

The results were consistent across both intervention arms. In the prophylactic group, all three quantitative measures of AACF were significantly lower than in the azaserine-only group, with statistical significance reported at p < 0.05. Critically, the therapeutic group, in which honey was introduced only after the carcinogenic process had been initiated, showed the same pattern of significant reduction. This means the honey did not merely shield healthy tissue from damage; it also appeared to slow or partially reverse the progression of lesions that had already begun to form. In parallel, the researchers compared serum levels of cancer antigen 19-9, the most widely used clinical biomarker for pancreatic cancer, and C-reactive protein, a general marker of inflammation, across the groups, adding a biochemical dimension to the histological findings.

What might explain such an effect? Honey is not a single compound but a chemically intricate matrix containing hundreds of substances, including polyphenols, flavonoids, hydrogen peroxide generated by the enzyme glucose oxidase, methylglyoxal, organic acids, and trace elements. Many of these constituents have documented antioxidant and anti-inflammatory activities in laboratory systems, and oxidative stress and chronic inflammation are recognized drivers of pancreatic carcinogenesis. Quercus pyrenaica honey, derived from an oak species native to parts of Europe and Anatolia, may possess a particularly favorable composition of these bioactive molecules, although the study did not perform a full chemical characterization of the honey used. The authors are careful to frame their work as a foundation for further investigation rather than a demonstration of mechanism.

The study’s limitations are those inherent to any animal experiment. The azaserine model, while valuable, does not recapitulate every feature of human pancreatic ductal adenocarcinoma, the most common and deadliest form of the disease in patients. Doses, timing, and species differences all complicate extrapolation to human beings, and the eleven-week follow-up window captures precursor lesions rather than invasive cancer or survival outcomes. The authors themselves acknowledge that their conclusions will gain real significance only if supported by clinical studies in this field. There is also the practical question of standardization: honey varies enormously by floral source, season, and geography, so any future therapeutic application would require rigorous quality control of the active constituents.

Nevertheless, the dual prophylactic and therapeutic signal is what makes the study notable. Many candidate chemopreventive agents show protection only when given before carcinogen exposure, a scenario with little clinical relevance, since people at risk of pancreatic cancer cannot usually be identified early enough. Demonstrating an effect when treatment begins after tumor initiation suggests that the bioactive components of Quercus pyrenaica honey may influence processes such as cell proliferation, apoptosis, or inflammation within already-transformed tissue, not just antioxidant shielding of healthy cells. If subsequent work identifies the responsible compounds and their molecular targets, it could open a path toward standardized extracts rather than raw honey as an intervention.

The research, which received ethics approval from the KONÜDAM Animal Ethics Committee of Necmettin Erbakan University and was conducted in accordance with the Declaration of Helsinki, was published open access with no external funding and no declared competing interests. For now, the message for the public is one of cautious fascination rather than immediate action: no one should interpret these findings as a reason to self-treat with honey, and patients with pancreatic disease should continue to rely on evidence-based oncological care. But the study adds a compelling data point to the argument that natural products deserve systematic scientific scrutiny, and it demonstrates that even a substance as ancient and familiar as honey can still surprise researchers when put to a properly controlled experimental test.

Subject of Research: Anti-tumoral effects of Quercus pyrenaica honey in an experimental rat model of pancreatic cancer

Article Title: Investigation of the anti-tumoral effect of honey in experimental pancreatic cancer

Article References: Kılınç, A., Ulutaş, M. E., Deniz, Ç. D., Öztaş, H., Hasırcı, İ., & Kartal, A. (2026). Investigation of the anti-tumoral effect of honey in experimental pancreatic cancer. BMC Complementary Medicine and Therapies. https://doi.org/10.1186/s12906-026-05565-2

Image Credits: AI Generated

DOI: 10.1186/s12906-026-05565-2

Keywords: honey, pancreatic cancer, Quercus pyrenaica, azaserine rat model, atypical acinar cell foci, CA 19-9, C-reactive protein, chemoprevention, experimental oncology, natural products, complementary medicine, preclinical research

The team focused on three oak species that dominate the climax forests of acidic substrates across the peninsula: cork oak (Quercus suber), Pyrenean oak (Quercus pyrenaica) and pedunculate oak (Quercus robur). These species sit at different points along a moisture gradient. Pedunculate oak, a temperate and humidity-loving species, thrives in the cool Atlantic north. Pyrenean oak, a marcescent species that holds its withered leaves through winter, occupies a transitional submediterranean zone. Cork oak, the iconic evergreen of Mediterranean silviculture, tolerates pronounced summer drought. Because the three species often coexist in the same landscapes, distributed along gradients from rocky ridge tops to moist valley bottoms, the researchers reasoned that climate change could trigger a domino effect: as each species retreats from its drier margins, a more drought-adapted congener may move in to take its place.

To test this idea, the researchers built high-resolution species distribution models at a one-square-kilometre resolution, integrating nearly a quarter of a million cleaned occurrence records from GBIF, national biodiversity databases and expert-curated sources. For pedunculate oak, the team applied a spatial thinning procedure that reduced more than 205,000 records to about 27,000, mitigating oversampling while preserving the species’ environmental coverage. The models combined nineteen bioclimatic indices from the CHELSA dataset with soil pH, soil texture, a Topographic Ruggedness Index and a Topographic Wetness Index, all chosen through hypothesis-driven screening that excluded collinear predictors.

The modelling framework was an ensemble approach implemented in the biomod2 package in R, combining eight techniques ranging from generalized linear models and generalized additive models to random forests, artificial neural networks and maximum entropy. Each model was validated over ten rounds of cross-validation, and only algorithms exceeding a True Skill Statistic threshold of 0.7 were averaged into the final ensembles. Performance was strong across the board: TSS values exceeded 0.83 and AUC scores surpassed 0.97 for all three species, with cork oak models performing best. Future projections drew on five general circulation models selected for their strong performance over Europe, run under two Shared Socioeconomic Pathways: an intermediate SSP3-7.0 scenario and a pessimistic SSP5-8.5 scenario, for the mid-century period 2041 to 2070 and the late-century period 2071 to 2100.

The environmental variables that best explained each species’ current distribution revealed their distinct ecological thresholds. For pedunculate oak, the minimum temperature of the coldest month was the dominant predictor, reflecting its tolerance of subzero winters, followed by topographic ruggedness and precipitation of the warmest quarter. Pyrenean oak’s distribution hinged on annual precipitation and warm-season rainfall together with soil pH, underscoring its dependence on water availability during the hot, dry summer. Cork oak, by contrast, was governed primarily by soil pH and mean annual temperature, favouring sites where winter minima stay above freezing. All three species are strictly silicicole, confined to acidic substrates, which means that climatically suitable areas on limestone soils will remain off-limits no matter how the climate shifts.

The projected trajectories diverge sharply. Pedunculate oak’s current range is expected to remain largely stable, with between 73 and 84 percent of its present suitable area persisting across scenarios, though its total suitable area declines slightly as losses outweigh gains. Cork oak emerges as the clear winner: its suitable range expands in every scenario and time frame, gaining between 32 and 44 percent of new territory, pushing northward and inland, with potential altitudinal advances reaching as far as the Pyrenees, while losing almost nothing at its southern edge. Pyrenean oak fares worst. Its total suitable area shrinks in all scenarios, with losses of up to nearly 60 percent of its current range concentrated in the southern and inland parts of the peninsula, and only limited compensatory gains at higher elevations.

When the researchers overlaid the range dynamics of different species, the turnover patterns became strikingly clear. Wherever one species loses suitability while another gains or holds steady in the same grid cells, a potential replacement is flagged. The dominant signal is a turnover from Pyrenean oak to cork oak, projected across roughly 20 to 47 percent of Pyrenean oak’s currently suitable range depending on the scenario, intensifying toward the end of the century and under harsher emissions. This replacement is expected to begin at the species’ southern fringes and spread into the interior. A smaller turnover, from pedunculate oak to Pyrenean oak, is projected across about 5 to 8 percent of the former’s range, mostly in low-altitude inland areas at the southern edge of its distribution. In some coastal zones of western and northeastern Iberia, cork oak may ultimately replace both of the other species where they both become unsuitable.

A crucial dimension of the study is its analysis of protected areas. Between a quarter and nearly half of the projected turnover zones fall inside classified conservation sites, including Natura 2000 areas designated under the European Union’s Habitats Directive, national parks and nature reserves. The turnover from Pyrenean oak to cork oak carries a direct policy consequence: it implies a shift from habitat 9230, the Galicio-Portuguese oak woods, to habitat 9330, cork oak forests, within sites whose management plans were written for the habitats of the past. The researchers argue that conservation targets in these sites will need periodic revision to maintain habitat continuity under changing climates, and that protected areas can serve as demonstration cases for adaptive management on the privately owned lands that dominate the Portuguese landscape.

The authors are careful to frame these projections as shifts in potential climatic and edaphic suitability rather than deterministic predictions of stand-level replacement. The models do not explicitly simulate dispersal, disturbance regimes or biotic interactions, and oak species are notoriously slow to migrate relative to the pace of climate change. Yet the convergent trends across multiple climate models and scenarios lend weight to the central inference: that a gradual Mediterraneanization of the Iberian landscape is underway, marked by longer summer dry periods and reduced annual precipitation, in which drought-adapted oaks progressively occupy the space vacated by their moisture-dependent relatives.

The practical implications reach deep into forestry. The researchers recommend protecting natural regeneration of the species likely to dominate in the future, since cork oak individuals already occur near Pyrenean oak stands on rocky outcrops and within forests. Where natural regeneration is insufficient, they point to seed sowing and planting techniques, including autumn sowing to maximize root development before summer drought, selection of larger acorns, protection from seed predators such as wild boar, and the use of large containers to avoid root deformation. They also advocate multi-provenance strategies that draw reproductive material from southern or low-altitude populations already adapted to drier conditions, combined with mixed-species stands to hedge risk. If guided by such adaptive management, the authors conclude, species turnover need not spell the loss of Iberian forests; it may instead be the mechanism by which they retain canopy cover, ecological function and resilience through a century of rapid environmental change.

Subject of Research: Projected climate-driven species turnover among Iberian oak forests

Article Title: Exploring turnover dynamics in Iberian Oak forests under climate change scenarios

Article References: Passos, I., Figueiredo, A., Gonçalves, J., Ribeiro, M. M., & Vila-Viçosa, C. (2026). Exploring turnover dynamics in Iberian Oak forests under climate change scenarios. Discover Conservation, 3(1), Article 14. https://doi.org/10.1007/s44353-026-00084-0

Image Credits: AI Generated

DOI: 10.1007/s44353-026-00084-0

Keywords: climate change, oak forests, Quercus suber, Quercus pyrenaica, Quercus robur, species distribution models, species turnover, Iberian Peninsula, Mediterranean forests, forest management, conservation, habitat suitability