BGPT: Paper Review: High royal jelly production does not impact the gut microbiome of honey bees

Fuel Your Discoveries

Quick Answer Copied

Bottom line (skeptical): royal jelly producing phenotype vs unselected stock shows little gut-microbiome divergence; urban–rural location effects are clearer.

In this 16S (V4) study, RJB vs ITB showed no strong differences in alpha- or beta-diversity; the only reported stock-level differential feature was Fructobacillus higher in ITBs.
Urban vs rural differences were statistically supported (e.g., Bray–Curtis and weighted UniFrac in beta-diversity), with non-core taxa largely driving composition shifts.

Key caveat: this is 16S amplicon data from a cross-sectional, geographically confounded design (urban vs rural also differ in sampling day/weather and are missing ITBs from rural). Those constraints limit inference about stock causality vs environment.

Long Answer

Paper review: “High royal jelly production does not impact the gut microbiome of honey bees”

DOI: 10.1186/s42523-021-00124-1 • Journal year: 2021 • Topic: honey bee gut microbiome vs high royal jelly producing phenotype

Study claim (as written): RJB vs ITB show no unique bacterial community differences; environment (urban vs rural) is more important in shaping gut microbial composition.

1) What they compared (design snapshot)

Population

RJB: royal jelly producing stock (6 hives urban, plus 6 hives rural).
ITB: Italian bees only at the urban site (6 hives).

Sample units

8 bees per hive from a brood frame (likely nurses), sampled individually for gut bacterial profiling.
Total collected: 144 bees; usable DNA extracted from 101 samples after QC/processing.

Microbiome assay

16S rRNA V4 amplicon (515F/806R) sequenced on Illumina iSeq100; analysis via DADA2/ASVs and QIIME2; contamination addressed via negative control removal/subtraction approach.

2) Visual summaries of reported statistical signals

How to read: lower p-values indicate stronger evidence of community differences in the reported analyses. The paper reports no significant Bray-Curtis difference for ITB vs RJB (urban-only comparison), while urban vs rural comparisons are significant.

3) What features differed (taxonomic signal)

Stock-level differential abundance (ANCOM)

RJB vs ITB (stock effect): one taxon reported as significantly different: Fructobacillus sp. higher in ITBs than RJBs.
Urban vs rural (environment effect): multiple taxa show higher abundance in urban bees vs rural bees; reported statistical significance is stated most strongly for Lactobacillus sp. at the taxon/clustering level, with additional patterns for Fructobacillus, Bifidobacterium, and others.

Important interpretive nuance (as written): the authors argue that the Fructobacillus difference between stocks likely arises from rural RJB having very low abundance of that taxon and then collapsing urban-only vs combined group effects.

4) Core conclusion vs what the data can/can’t prove

What the evidence supports (in-scope)

Within their sampled design, the authors did not find strong community differences attributable to royal jelly production stock identity (RJB vs ITB), using diversity metrics and differential abundance tests.
Location (urban vs rural) shows clearer differences in both diversity structure and differential abundance, with core taxa relatively invariant and non-core taxa contributing more to differences.

What the evidence does not establish (key limitations)

Geographic confounding / incomplete counterfactuals: ITBs were not sampled in rural China, so “stock vs environment” separation is limited.
Temporal/day-of-sampling confounds: urban and rural were collected on different days with different weather (overcast vs sunny) and temperatures; these factors could influence what microbes the bees acquire (e.g., foraging exposure).
Resolution limits of 16S ASVs: 16S amplicon sequencing cannot reliably resolve strain-level variation or functional capabilities, which could in principle differ between RJB and ITB even if community-level composition looks similar.
Taxon presence/absence & contamination handling: the authors remove unique negative-control ASVs and subtract counts for core ASVs; however, residual contamination, index hopping, or differential DNA extraction efficiency could still influence low-abundance taxa (especially when “differences” are reported as low relative abundances).

5) Mechanistic interpretation: what could reconcile “no stock effect” with biological plausibility?

Royal jelly production may be primarily driven by host gland physiology/genetics, with any microbial contribution being modest enough to not shift broad community composition detectable by 16S diversity metrics. The paper’s motivation cites host physiological and molecular differences associated with RJ production, but the microbiome link was not supported by strong composition differences in their data.
Microbiome effects could be strain- or function-specific rather than taxon-abundance driven—precisely the gap that 16S amplicon profiling leaves open.

6) What would most efficiently disprove or revise their conclusion?

Find a stock effect under matched environments: e.g., sampling ITBs and RJBs in both urban and rural (or transplant/common-garden approaches), with identical weather/time-of-sampling windows, then re-testing diversity and differential abundance.
Use shotgun metagenomics/metatranscriptomics to test for functional divergence between RJB and ITB even when 16S composition appears similar.

7) Quick author review links (from BGPT)

Feedback:

Updated: May 02, 2026