Kwantificeren van de dansvloer van de honingbij: een datagedreven methode voor het definiëren en vergelijken van regio's voor de kwispeldans

Samenvatting

Foerageerders van de honingbij (Apis mellifera) voeren in de korf kwispeldansen uit om de locatie van rendabele voedselbronnen aan hun nestgenoten te communiceren. Deze wervingsdansen vinden plaats binnen een specifiek gebied van de raat, bekend als de 'dansvloer', maar de locatie en structuur hiervan zijn historisch gezien alleen kwalitatief beschreven.

In dit onderzoek introduceren we een datagedreven geometrische methode om de dansvloer te definiëren en te kwantificeren op basis van coördinaten van de kwispeldans. De aanpak combineert convex hulls (convexe omhullingen) en betrouwbaarheidsellipsen om een gesloten regio te creëren die het gebied met de hoogste dansdichtheid vertegenwoordigt. Dit levert interpreteerbare ruimtelijke statistieken op, waaronder de locatie van het zwaartepunt (centroid), oppervlakte, omtrek, de hoofd- en kleine as, en de oriëntatie.

Om de prestaties van de methode te demonstreren, hebben we deze toegepast op 155 waarnemingen van acht kolonies in observatiekasten die verschilden in grootte en datum. Met behulp van complementaire univariate en multivariate analyses bleek het raamwerk consistent ongeveer 91% van de dansen te vatten, wat overeenkomt met historische schattingen gebaseerd op de afstand tot de ingang. Bovendien detecteerde de methode systematische verschillen tussen waarnemingen gerelateerd aan kastgrootte, dag en tijd (bijv. grootte-afhankelijke verschuivingen in de positie van het zwaartepunt en de breedte, en effecten van het tijdstip van de dag op de oriëntatie). Dit illustreert de gevoeligheid van de methode voor experimentele en temporele context. Dit werk biedt een expliciete kwantitatieve definitie van de dansvloer van de honingbij en een reproduceerbaar analytisch raamwerk voor het vergelijken van ruimtelijke wervingspatronen tussen kolonies, omgevingen en toekomstige experimentele ontwerpen.

---

Inleiding

Een van de meest intrigerende vormen van dierlijke communicatie is de kwispeldans van honingbijen (Apis mellifera). Succesvolle foerageerders adverteren aan hun nestgenoten een mathematische vector (richting en afstand) naar een waardevolle bron via een reeks complexe manoeuvres op de raat binnen de donkere korf. Deze wervingsdansen vinden niet willekeurig in het nest plaats, maar zijn doorgaans geconcentreerd in een specifiek gebied, bekend als de dansvloer. De ruimtelijke lokalisatie van dit gedrag verhoogt waarschijnlijk de efficiëntie van de communicatie door zowel werkloze foerageerders als ervaren, terughoudende foerageerders in dezelfde fysieke regio te concentreren waar actieve rekruteerders en succesvolle verkenners hun dansen uitvoeren. Toch is, ondanks decennia van onderzoek, nog steeds onduidelijk welke factoren de locatie en structuur van de dansvloer bepalen.

In zowel natuurlijke als beheerde nesten vinden kwispeldansen meestal plaats nabij de ingang op vrijhangende, lege raat, wat de overdracht van trillingen kan verbeteren. Echter, vrijhangende, lege raat is geen vereiste; dansen komen ook voor in commerciële en observatiekasten op raten die volledig zijn geïntegreerd in houten frames, inclusief broedraten. In vrijwel alle tot nu toe uitgevoerde studies wordt de dansvloer simpelweg beschreven als zijnde "nabij de ingang van de korf". Zo merkten Seeley en Towne op dat 94% van de dansen in een observatiekast met twee frames werd uitgevoerd binnen 24 cm van de ingang, met de hoogste dichtheid tussen 4 en 18 cm. Andere voorgestelde factoren zijn onder meer geurmarkering door vroege terugkerende foerageerders, het vrijgeven van chemische signalen die werving stimuleren, en voorkeuren voor specifieke substraattypes zoals broedraten.

Ondanks deze interesse is de dansvloer zelden kwantitatief gedefinieerd. De meeste studies beschrijven deze alleen kwalitatief, wat de mogelijkheid beperkt om ruimtelijk wervingsgedrag systematisch te vergelijken tussen kolonies, contexten of experimentele condities. Onze geometrische aanpak genereert een expliciete, datagedreven grens die direct vergeleken kan worden. De methode combineert twee geometrische instrumenten: een convex hull, die de buitenste grens van alle danspunten traceert, en een betrouwbaarheidsellipse, die de statistische spreiding samenvat. Samen vormen zij een gesloten regio die de kern van de dansactiviteit omvat.

In dit onderzoek passen we deze methode toe op data van acht honingbijkolonies in glazen observatiekasten van twee verschillende groottes. We onderzoeken of de structuur van de dansvloer:

1. Varieert gedurende de dag;
2. Verandert van dag tot dag;
3. Verschilt tussen kolonies;
4. Verschuift gedurende het seizoen;
5. Afhangt van de grootte van de kast.

---

Materialen en methoden

Dieren, experimentele opzet en gegevensverzameling

Kolonies van de westerse honingbij (Apis mellifera L.) werden gehuisvest in observatiekasten met glazen wanden van 8 frames (100 × 100 cm) of 4 frames (100 × 50 cm) in een gematigd gebied. De 8-frame kasten stonden in een houten schuur in een veld met wilde bloemen en bomen, waarbij natuurlijk licht werd geblokkeerd en diffuus kunstlicht werd gebruikt. De 4-frame kasten stonden in een laboratorium met gecontroleerde licht- en temperatuuromstandigheden op een academische campus.

Elke kolonie werd één keer gebruikt in een proef die 3 tot 5 dagen duurde. De ingangen van de kasten werden aangepast zodat binnenkomende bijen alleen toegang hadden tot één zijde van het frame dichtst bij de ingang, waardoor alle kwispeldansen werden beperkt tot één raatoppervlak. Een raster van 5 × 5 cm op de glazen zijden hielp waarnemers bij het schatten van de locaties. In totaal werden 7.444 dansen geanalyseerd.

Tijdens elke observatieperiode werden de locaties van alle kwispeldansende bijen geregistreerd door het raster in een horizontale zigzagbeweging te scannen, beginnend bij het kwadrant het verst van de ingang. De scanprocedure werd ongeveer 20 keer herhaald per observatieperiode van 30 minuten. Waarnemingen vonden elke 2 uur plaats. De kaarten werden gedigitaliseerd om de horizontale en verticale afstand van elke dans tot de ingang te berekenen.

Definitie van de dansvloer

Alle stappen werden uitgevoerd in R versie 4.2.1. De dansvloeren werden gedefinieerd in drie stappen:

1. Convexe omhulling (Convex Hull): De kleinste convexe polygoon die alle danslocaties van een waarneming omsluit (vergelijkbaar met een elastiek dat rond de buitenste punten wordt gespannen). Dit definieert de buitenste grens van de distributie.
2. Betrouwbaarheidsellipse (Confidence Ellipse): Om de invloed van uitschieters te beperken, werd een ellipse berekend op basis van de ruimtelijke distributie. De ellipse is gecentreerd op de gemiddelde locatie, met assen geschaald naar twee standaarddeviaties (2-SD). Dit vat de kern van de ruimtelijke spreiding samen.
3. Ruimtelijke intersectie: De dansvloer werd gedefinieerd als de ruimtelijke doorsnede (intersectie) van de convex hull en de 2-SD betrouwbaarheidsellipse. Deze aanpak minimaliseert overbodige oppervlakte wanneer een klein aantal dansen ver van de hoofdbundel plaatsvindt.

Statistieken van de dansvloer

Uit elke dansvloerpolygoon werden zeven statistieken geëxtraheerd:

• Oppervlakte (cm²)
• Omtrek (cm)
• X en Y (horizontale en verticale coördinaten van het zwaartepunt ten opzichte van de ingang, in cm)
• Lengte en Breedte (hoofdas en kleine as van de betrouwbaarheidsellipse, in cm)
• Hoek (oriëntatie van de hoofdas in graden ten opzichte van de horizontale lijn)

Statistische analyses

De analyses werden uitgevoerd op het niveau van individuele observatieperioden van 30 minuten. Om consistentie te bevestigen, werd een lineair model gebruikt om de proportie gevangen dansen tussen 4-frame en 8-frame kasten te vergelijken.

Univariate analyses: Er werden lineaire mixed-effects modellen gebouwd voor elke afhankelijke variabele. De vaste effecten waren Tijd (uren sinds zonsopgang), Proefdag (dag 1–5), en Dag van het Jaar (DOY), inclusief hun interacties met de Grootte van de kast (4-frame of 8-frame).

Multivariate analyse: Er werd een MANOVA-model geconstrueerd om de algehele configuratie van de dansvloer te beoordelen. Hierbij werd Grootte als primaire verklarende variabele behouden.

---

Resultaten

De meeste dansvloeren waren min of meer elliptisch en bevonden zich nabij de ingang van de korf. Oppervlakte en omtrek waren sterk gecorreleerd, waardoor alleen Oppervlakte werd behouden voor verdere analyses.

Vorm van de dansvloer

De vorm van de dansvloer werd beïnvloed door meerdere voorspellers:

• Oppervlakte: Varieerde significant per proefdag, maar zonder consistent patroon.
• Lengte: Werd beïnvloed door zowel de proefdag als de kastgrootte; dansvloeren in 4-frame kasten waren over het algemeen langer dan in 8-frame kasten.
• Breedte: Toonde een significante interactie tussen kastgrootte en dag van het jaar (DOY). In 8-frame kasten piekte de breedte midden in het seizoen, terwijl deze in 4-frame kasten juist midden in het seizoen het laagst was en later toenam.
• Hoek (Oriëntatie): De enige vormparameter die werd beïnvloed door de tijd van de dag. In 4-frame kasten bleef de oriëntatie stabiel, terwijl deze in 8-frame kasten sterk varieerde gedurende de dag.

Positie van de dansvloer

De positie van het zwaartepunt (centroid X en Y) varieerde significant op basis van DOY, kastgrootte en hun interactie:

• X-coördinaat (horizontaal): In 4-frame kasten verschoof de positie later in het seizoen naar rechts, terwijl deze in 8-frame kasten stabiel bleef.
• Y-coördinaat (verticaal): Steeg gedurende het seizoen in 4-frame kasten, maar bleef consistenter in 8-frame kasten.

Dit suggereert dat kleinere kolonies de plaatsing van hun dansactiviteit dynamischer aanpasten aan de hand van het seizoen, waarbij de dansvloer later in het seizoen verder van de ingang kwam te liggen.

Multivariate resultaten en robuustheid

De MANOVA bevestigde dat kastgrootte, DOY en hun interactie sterke effecten hadden op de algehele structuur van de dansvloer. De interactie tussen kastgrootte en tijdstip van de dag was ook significant, hoewel zwakker.

De definitie van de dansvloer ving consistent de meerderheid van de kwispeldansen op, met een gemiddelde inclusieratio van 90,8% ± 4,9%. Er was geen significant verschil in deze ratio tussen de verschillende kastgroottes.

---

Discussie

Onze geometrische definitie biedt een nieuw raamwerk voor het kwantificeren van de ruimtelijke organisatie van wervingsdansen. De resultaten tonen aan dat de structuur van de dansvloer systematisch varieert op basis van tijd en kastcondities.

Analyse van ruimtelijke dynamiek

De dynamiek van de verticale positie van het zwaartepunt suggereert dat grote en kleine kasten de danslocaties op verschillende tijdschalen aanpassen. Dit is consistent met eerdere waarnemingen dat dansgedrag flexibel is en reageert op veranderingen in de koloniestructuur, drukte en omgevingscontext.

De divergentie in breedte tussen kastgroottes suggereert dat kleinere en grotere kolonies de ruimte anders reorganiseren in reactie op seizoensgebonden druk. Opvallend is dat de oriëntatie (hoek) alleen in grotere kasten gedurende de dag verschoof. Dit kan erop wijzen dat beperkte raatoppervlakte in kleine kasten de bijen dwingt tot een stabiele oriëntatie, terwijl grotere kasten meer flexibiliteit bieden in reactie op veranderende verkeerspatronen of drukte.

Interpretatie van de resultaten

De MANOVA-resultaten suggereren dat seizoensgebonden en contextuele factoren de dansvloer als een geïntegreerd systeem vormen, en niet slechts als een verzameling onafhankelijke statistieken. De hittekaarten bevestigen dat dansactiviteit sterk geconcentreerd blijft in een relatief beperkt gebied nabij het centrum van het frame dichtst bij de ingang.

Kastgrootte bleek de meest consistente voorspeller. Hoewel we niet met zekerheid kunnen zeggen of de verschillen enkel door de grootte komen (omdat 8-frame kasten in 2012 in een schuur stonden en 4-frame kasten in 2013 in een laboratorium), bevestigt dit de bruikbaarheid van de methode om structurele variatie te detecteren.

Theoretische context en beperkingen

Hoewel de exacte mechanismen onduidelijk blijven, is het waarschijnlijk dat bijen informatie integreren uit chemische signalen (geurmarkeringen) en een voorkeur hebben voor dansen op broedcellen. Daarnaast veroorzaakt de 'drift' (het verschuiven van de positie tijdens het dansen) dat dansen voor verre bronnen gemiddeld verder van de ingang plaatsvinden dan dansen voor nabije bronnen.

Beperkingen van dit onderzoek zijn onder meer het gebruik van niet-gemarkeerde bijen, waardoor individuele variabiliteit niet kon worden vastgesteld, en het feit dat observatiekasten kleiner zijn dan commerciële Langstroth-kasten.

Conclusies op de centrale onderzoeksvragen

1. Verandert de dansvloer gedurende de dag? Ja, de hoek varieerde in grote kasten, maar bleef stabiel in kleine kasten.
2. Verandert de dansvloer van dag tot dag? Ja, proefdag had een significant effect op oppervlakte, lengte en verticale positie.
3. Verschillen dansvloeren tussen kolonies? Ja, de identiteit van de kolonie beïnvloedde alle statistieken.
4. Beïnvloedt het seizoen de dansvloer? Ja, de dag van het jaar voorspelde de horizontale en verticale positie en de breedte (afhankelijk van kastgrootte).
5. Beïnvloedt de kastgrootte de dansvloer? Ja, dit was het meest consistente effect, met significante invloed op lengte, breedte en positie.

Onze datagedreven aanpak biedt een schaalbaar instrument voor toekomstig onderzoek naar de lokalisatie van de kwispeldans, waardoor consistente vergelijkingen mogelijk worden tussen kolonies, omgevingen en experimentele ontwerpen.