
==== Front
Biol Open
Biol Open
BIO
Biology Open
2046-6390
The Company of Biologists Ltd

39263864
10.1242/bio.060511
BIO060511
Version of Record
Research Article
Development of germline progenitors in larval queen honeybee ovaries
Cullen Georgia
Delargy Erin
http://orcid.org/0000-0001-7790-9675
Dearden Peter K. peter.dearden@otago.ac.nz
*
Genomics Aotearoa, Bioprotection Aotearoa and Biochemistry Department, University of Otago, Dunedin 9016, Aotearoa, New Zealand
* Author for correspondence (peter.dearden@otago.ac.nz)
Competing interests

The authors declare no competing or financial interests.

15 9 2024
12 9 2024
12 9 2024
13 9 bio06051129 4 2024
8 8 2024
© 2024. Published by The Company of Biologists Ltd
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

ABSTRACT

Honeybees (Apis mellifera) are a keystone species for managed pollination and the production of hive products. Eusociality in honeybees leads to much of the reproduction in a hive driven by the queen. Queen bees have two large active ovaries that can produce large numbers of eggs if conditions are appropriate. These ovaries are also active throughout the long lives of these insects, up to 5 years in some cases. Recent studies have indicated that the germline precursors of the adult honeybee queen ovary are organized into 8-cell clusters, joined together by a polyfusome; a cytoplasmic bridge. To understand the origin of these clusters, and trace the development of the honeybee queen ovary, we examined the cell types and regionalization of the developing larval and pupal queen ovaries. We used established (nanos and castor), and novel (odd skipped) gene expression markers to determine regions of the developing ovary. Primordial germline cells develop in the honeybee embryo and are organized into ovary structures before the embryo hatches. The ovary is regionalized by larval stage 3 into terminal filaments and germaria. At this stage, clusters of germline cells in the germaria are joined by fusomes and are dividing synchronously. The origin of the 8-cell clusters in the adult germarium is therefore during larval stages. On emergence, the queen ovary has terminal filaments and germaria but has not yet developed any vitellaria, which are produced after the queen embarks on a nuptial flight. The lack of germaria, and the storing of germline progenitors as clusters, may be adaptions for queen bees to endure the metabolic demands of a nuptial flight, as well as rapidly lay large numbers of eggs to establish a hive.

Summary: Honeybee reproduction is key to ensuring the productivity on hives and the maintenance of pollination capacity. Understanding how queen ovaries develop from embryo to adult provides information about ovary stucture, and the nature and formation of the gremlin.

Honeybee (Apis mellifera)
Germline
Ovary
University of Otago http://dx.doi.org/10.13039/100008247
==== Body
pmcINTRODUCTION

Apis mellifera (Honeybees) are eusocial insects, characterised by cooperative brood care, overlapping generations and a division of reproductive labour (Wilson and Hölldobler, 2005). Individuals undergo enforced altruism (Ratnieks and Halenterä, 2009), becoming aggressive to workers with more developed ovaries (Visscher and Dukas, 1995) and killing pathogen-infected (Starks et al., 2000). The two female castes of honeybees, workers and queens, are exemplars of the polyphenism underlying the division of reproductive labour. Queen bee anatomy differs from workers, especially in the case of the abdomen, which is larger in queens, and contains two large active ovaries (Snodgrass, 1956). The queen's ovaries are divided into 116-219 ovarioles each (Jackson et al., 2011), while workers have less than 10 per ovary (Snodgrass, 1956). This polyphenism in ovary size is a result of differential feeding during larval developmental stages (Rachinsky et al., 1990). After eclosion, honeybee queens remain in the hive for a few days (Plate et al., 2019) to a week (Jackson and Robinson, 2018), before undertaking nuptial flights (Woyke, 1964). After mating, queens usually take a few days to begin laying (Kocher et al., 2008). The queen produces thousands of embryos that become either workers or drones depending on fertilization (Woyke, 1963). She is the only individual in a colony that can produce fertilized eggs, and thus female workers. Therefore, the queen's ability to lay is crucial for long-term hive success (Gill and Hammond, 2011).

Honeybees have polytrophic meroistic ovaries (Gutzeit et al., 1993; Tanaka and Hartfelder, 2004). This ovary structure is defined by the presence of nurse cells, which provide nutrients and materials to the developing oocyte, and accompany it along the ovariole (Büning, 1993). In honeybees, the ovariole is divided into three major sections: the terminal filament, germarium and vitellarium (Cullen et al., 2023; Dearden, 2006). The terminal filament consists of flattened somatic cells, that appear to provide the progenitors of the somatic follicle cells of the ovary (Hartfelder and Steinbrück, 1997; Martins et al., 2011; Tanaka and Hartfelder, 2004). The germarium contains the first germ-line progenitors, which are clusters of eight germline cells joined by a polyfusome (Cullen et al., 2023). These 8-cell clusters can divide synchronously at a rate that maintains the numbers of these clusters between laying seasons (Cullen et al., 2023).

The progenitor 8-cell clusters in the germarium will divide multiple times (Cullen et al., 2023). Upon entering the posterior half of the germarium, polyfusomes cease and ring canals, another cytoplasmic bridge, appear (Gutzeit et al., 1993). Eventually, germline clusters become oocytes and accompanying nurse cells (Cullen et al., 2023; St Johnston, 2008). This, and the further maturation of the oocyte, takes place in the vitellarium (Chapman, 1998). The oocyte must undergo meiosis, and accompanying nurse cells will divide and endoreplicate, producing polyploid nurse cells linked to each oocyte via the ring canals. The oocytes will then grow, and become provisioned with yolk until they are ready to be laid (Aamidor et al., 2022). The 8-cell cluster is a key point in reproduction as these are the primary store of germline cells in the ovary and the source of oocytes.

Despite the importance of honeybees to the economy and the environment (Khalifa et al., 2021; Tanaka and Hartfelder, 2004), honeybee populations are declining (Brosi et al., 2017; Yang et al., 2023). They experience a range of threats, such as parasites (Warner et al. 2024), viruses (Abbo et al., 2017; Kang et al., 2016), insecticide exposure (Manzoor and Pervez, 2022; Stuligross and Williams, 2021), and the reduction of natural ecosystems and diverse flora for agriculture and urbanization (Goulson et al., 2015). Due to these threats, and others, the costs of beekeeping are increasing. A better understanding of how queen honeybees consistently produce large numbers of offspring, and how to optimise conditions to support that reproduction, will help reduce the effects these threats have on the global populations.

To understand how to best support honeybee reproduction we need to understand the dynamics of germline and oocyte production. Previous studies have shown that germline cells are specified late in honeybee embryogenesis (Dearden, 2006) and that the bulk of ovary growth occurs during larval and pupal stages (Chapman, 1998; Hartfelder and Steinbrück, 1997).

Here we set out to trace ovary development in queen honeybees. Understanding when and where ovary structure appears, how the ovarioles form, and, perhaps most importantly, where/when the 8-cell clusters form, is crucial information if we are to better support reproduction in queen bees and understand the evolution of the unusual structures of the honeybee ovary.

RESULTS

Timing of development of honeybee queen ovary morphology

The precursors to honeybee ovaries form in the embryo as two lines of cells in the most dorsal regions of stage 9 embryos, between abdominal segments 3 and 6 (Dearden, 2006). These consolidate into two arches of cells underlying the dorsal surface of just-hatched (L1) larval honeybees (Dearden, 2006). Considerable growth and differentiation occur during larval and adult stages.

Stage 1 and 2 ovaries are difficult to visualize and dissect out of the larval body and do not stain well if kept as whole mounts. L2 larval ovaries lie on either side of the midline between abdominal segments 3 and 6. The ovaries are small and transparent. The cells of the ovary are linked with trachea and buried in the fat layer of the L2 larva (Fig. 1A).

Fig. 1. A developmental series of Honeybee (Apis mellifera lingustica) queen ovaries. The ovaries are oriented with the developing terminal filaments to the left, forming an oviduct to the right. Ovaries stained for DNA with DAPI (grey) and cortical actin with phalloidin (magenta). Images are organized by stage: (A) larval stage 2, white line demarks the limits of the ovary tissue, (B) larval stage 3, (C) larval stage 4, (D) early larval stage 5, (E) late larval stage 5, (F) early pupa, (G) late pupa, yellow dotted lines show crop lines due to the microscope field of view, and rotation of the specimen to match other samples (H) newly emerged queen. Scale bars: 200 μm.

L3 queen larval ovaries have contracted somewhat and are now found between abdominal segments 4 and 5 on either side of the dorsal line (data not shown). When stained with DAPI and phalloidin, the cells are organized into ‘fingers’ of tissue that resemble small future ovarioles, although they are not separated from each other (Fig. 1B).

From this stage onward, the ‘fingers’ of tissue elongate, and, by early L5, separate into individual ovarioles. The ovarioles lengthen through late larval stage 5, with connective tissue forming to connect the posterior of all the ovarioles (Fig. 1E-G). This is likely tissue that will go on to be, or connect to, calices (Kozii et al., 2022). Calices provide the transfer of eggs between the ovaries and the lateral oviduct (Kozii et al., 2022). This region connects both lateral oviducts from either ovariole with the spermatheca (median oviduct) (Kozii et al., 2022).

Growth continues in the pupal stages, with terminal filaments and germaria present. On hatching, the ovary does not appear to contain nurse cells or differentiated oocytes.

Regionalisation of the ovary during development

To understand the regions of the developing honeybee ovaries, we used markers of specific cell types to identify each region as it develops. In adults, ovary cells that express castor (cas) are somatic (Cullen et al., 2023) and so we used this gene as a potential marker of somatic cells in the ovary. Germline cells express nanos (nos) and vasa (vas) (Dearden, 2006) from late embryonic stages through to adults (Cullen et al., 2023), allowing us to use these as markers of germline fate. RNA expression from the nos gene was used in these experiments to mark germline cells

The odd skipped (odd) gene encodes a zinc finger transcription factor that in Drosophila melanogaster acts as a pair-rule gene (Coulter et al., 1990). Odd is a member of a family of similar transcription factors encoded by the brother of odd with entrails limited (bowl) and sister of odd (sob) genes (Hart et al., 1996). To determine the most likely orthologue of odd in honeybees, we identified homologues of all these genes using blastp (Altschul et al., 1990), and then carried out maximum likelihood phylogenetics. Honeybee XP_001120949 clusters with odd proteins from other holometabolous insects against sob and bowl proteins (Fig. S1). We thus designate the gene encoding this protein, honeybee odd skipped (odd).

In Drosophila, RNA from odd is expressed in a variety of cell types, but not germline or ovary cells (Coulter et al., 1990; Ward and Coulter, 2000). In experiments to examine the expression of genes involved in segmentation in honeybees (data not shown), we identified that honeybee odd is expressed during segmentation, but is also expressed in the same presumptive ovary region as nos and vas RNA in late embryos (Fig. 2). This region, located in the most dorsal regions of the germband, between segments 3 and 6, is where primordial germline cells can first be detected via expression of nos and vas RNA (Dearden, 2006).

Fig. 2. Expression of honeybee odd skipped in embryos and adult ovaries. (A) Stage 8 honeybee embryo (Cridge et al., 2017) stained for even skipped (magenta) for reference (Wilson and Dearden, 2012) and odd skipped (green) (Fig. S1). RNA from the odd skipped gene is expressed in clusters of cells interspersed with even skipped expression along the body wall of the embryos, and in a line of cells stretching from abdominal segment 3 to abdominal segment 6 in the most dorsal region of the germ band (marked by yellow line). (B and C) Expression of odd skipped (green) in terminal filament regions of an adult queen ovary. RNA expression appears in all the cells of the terminal filament (in this specimen folded over the germarium of the ovary), with expression reducing at the boundary of the germarium. DNA stained with DAPI, grey in C.

To attempt to determine what cells odd might be marking in the honeybee embryonic ovary, we investigated where odd is expressed in the adult queen ovary. RNA from odd is expressed in the cells of the terminal filament (Fig. 2). It is detected in the cells around the boundary of the terminal filament and the germarium but absent in the follicle cells surrounding the germarium (Fig. S2). These cells are thought to be the descendants of terminal filament cells (Dobens and Raftery, 2000). This encouraged us to use odd as a marker for presumptive terminal filament cells during ovary development.

At larval stage 3, odd, cas and nos RNA expression can be used to identify the cell-type structure of the developing ovary (Fig. 3). The ovary is a curved structure with the outside of the curve oriented to the midline of the larva. Many of the cells of the ovary appear to express castor RNA, especially the trachea that infiltrate the tissue. Expression of odd RNA appears in clusters of cells on the outside of the curved ovary. RNA from nos appears in groups of cells juxtaposed to odd expressing cells. Close analysis of these expression patterns (Fig. 3F, G and H) show odd RNA expression at the end of the small fingers of tissue in the ovary at this stage. RNA from nos is expressed in a band of cells abutting the odd domain further down the structure. Cells do not co-express odd and nos RNA. This relationship is very similar to the final pattern of gene expression of odd and nos, with odd RNA being expressed in the terminal filament, and its expression ceasing at the boundary of the germarium, where groups of nos RNA-expressing cells first appear in the ovary. These expression patterns imply that the finger structures are presumptive ovarioles, and these are initially separated into terminal filament and germline cells.

Fig. 3. Expression of RNA from cas (yellow), nos (cyan) and odd (magenta) in L3 larval ovaries. (A) Ovary stained with DAPI (grey), cas, nos and odd. (B) cas, nos and odd. (C) cas and odd. (D) cas and nos. (E) nos and odd. (F, G and H) Magnification of expression of (F) DAPI, cas, nos and odd. (G) cas, nos and odd and (H) nos and odd. Expression of odd RNA (magenta) can be seen at the end of the small ‘fingers’, with nos RNA expression in a cluster of cells abutting the odd domain. Scale bars indicate 100 μm.

Larval stage 4 ovaries show a very similar arrangement as those at L3 (Fig. 4). Though the ‘fingers’ are slightly longer. By early larval stage 5, however, the ‘fingers’ have separated into ovarioles, with the terminal tips of the ovarioles expressing odd and cas RNA, nos RNA being expressed by clusters of cells further down each ovariole, and cas RNA expressing cells surrounding these. No cells co-express cas and nos RNA.

Fig. 4. Expression of cas (yellow), odd (magenta) and nos (cyan) throughout larval and pupal development of queen honeybees. DNA stained with DAPI is shown in grey. (A-E) L4 ovaries, (F-J) Early L5 ovaries, (K-O) Late L5 ovaries, (P-T) Early pupal ovaries and (U-Y) Late pupal ovaries, with staining for DAPI, cas, odd and nos (A), cas, odd and nos (B), cas and odd (C), cas and nos (D) and odd and nos (E). Other stage ovaries are stained with the same genes and displayed in the same colours. Scale bars: 100 μm.

As the larval ovary extends during L5, the arrangement of cells seen in early L5 is maintained but extended. Odd and cas RNA-expressing cells make up the developing terminal filament. Moving along each ovariole, in the region in which clusters of nos RNA expressing cells appear, odd RNA-expressing cells disappear, with cells surrounding the nos RNA expressing cells expressing only cas RNA. These structures and gene expression patterns remain throughout larval and pupal development as the ovary grows. By analogy to the adult form, the somatic and germline cells of the ovary are defined in the embryo, they become organized by L3, and then divide and grow to produce the structures of the terminal filament and germarium in the adult ovary. There is no evidence of a vitellarium in larval or pupal ovaries.

These in situ hybridization data indicate that there are germline cells in groups at L3. However, it does not tell us whether these germline cells are linked via cytoplasmic bridges, and divide synchronously, as in adult germline cell clusters.

Formation of intercellular bridges

Having described the regionalisation of the ovary into terminal filaments and germaria, we next investigated when and where the germline cells form into the 8-cell clusters that comprise the germline in adult queens (Cullen et al., 2023). These clusters have two key characteristics. The first is that they are joined by a polyfusome, a cytoplasmic bridge that links the 8 cells together (Cullen et al., 2023). The second characteristic is that they divide synchronously (Cullen et al., 2023). To determine when 8-cell clusters begin to form in the ovary, we stained for polyfusomes with phalloidin (Cullen et al., 2023), and for dividing cells with antibodies directed to phosphorylated histone H3 (Hans and Dimitrov, 2001).

Fusomes and polyfusomes are cytoplasmic intercellular bridges between cells forming cystocytes and germline cell clusters (Ventelä, 2006). In Drosophila melanogaster, fusome precursors are called spectrosomes (Deng and Lin, 2001). The cystoblast spectrosome grows into a fusome, a distinct region of the cytoplasm, which bridges two germline cells (King, 1979), as germline cysts form (de Cuevas and Spradling, 1998). Polyfusomes, in Drosophila, are assembled from the further joining of fusomes creating branched structures that can link together the cytoplasm of several cells (King, 1979; Pritsch and Büning, 1989; Trauner and Büning, 2007). The appearance of polyfusomes during development implies the formation of germline cysts (Cullen et al., 2023).

Phalloidin staining should mark spectrosomes, fusomes, polyfusomes and ring canals. In the earliest larval stage ovaries we can image with confidence, larval stage 2, we cannot find accumulations of phalloidin staining consistent with any of these structures (Fig. 5A). By larval stage 3, however, phalloidin can be seen staining concentrations of material that appear to link neighbouring cells (Fig. 5B), reminiscent of Drosophila fusomes. These potential honeybee larval stage 3 fusomes are 5-7.6 μm wide comparable to Drosophila fusomes, (6-10 μm wide during cystoblast to 2-cell cyst stages) (de Cuevas and Spradling, 1998). The similarity in size, staining and shape implies the presence of fusomes at larval stage 3, and that some of the germline cells have formed 2 or 4-cell cysts.

Fig. 5. Phalloidin and DAPI stained larval and pupal ovaries focusing on the appearance of cytoplasmic bridges. (A) Stage L2 ovary. Phalloidin stains cortical actin around cells, but there is no sign of conglomerations of staining implying the formation of fusomes. (B) Stage L3, structures staining strongly with phalloidin can be seen bridging some neighbouring cells (marked with yellow arrowheads). These structures are generally to be found in cells adjacent to cells at the tips of the ‘finger’ structures in the ovary (marked where visible with yellow asterisks). (C) L4 ovariole, strongly staining polyfusome-like structures (marked with yellow arrowheads) can be seen several 10s of micrometres away from the terminal structures of the ovariole. (D) Early-stage L5 ovaries; phalloidin-stained polyfusome-like structures can be seen associated with the forming ovarioles. (E) Late-stage L5 ovaries; rows of phalloidin-stained polyfusome-like structures can be seen in each ovariole. (F) Pupal ovaries, polyfusomes structures can be seen as in late L5 ovaries (data not shown), and in more posterior regions, polyfusomes have disappeared and ring canals (indicated with arrows) are now visible (phalloidin staining rings in E). Scale bars: 20 μm.

Polyfusomes, similar in size and shape to those in adult ovaries, are present from larval stage 4 (Fig. 5C,D,E). Ring canals, which are present in the posterior germarium of adult queen honeybee ovaries, are not present in the ovary until late pupal stages (Fig. 5F).

Phalloidin staining implies that fusomes, and thus the development of cysts, appear from larval stage 3. If these are the developing germline cell clusters, we should be able to detect synchronously dividing cells at the same stages.

Synchronously dividing cells

Germline cell clusters in adult honeybee ovaries are joined by polyfusomes and divide synchronously (Cullen et al., 2023). Synchronously dividing cells, especially if associated with polyfusomes, in the developing ovary may be signs of clusters forming. As noted (Fig. 1), the ovaries grow significantly during larval and pupal stages. This growth is associated with large amounts of cell divisions in all tissues of the developing ovary (Fig. 6) making identifying synchronously dividing cells challenging.

Fig. 6. Immunohistochemistry for PHH3 in larval and pupal ovaries to detect dividing cells. Ovaries are stained with PHH3 (magenta) and phalloidin (cyan) to show where dividing cells appear in the developing ovary. As the ovary grows, many cells in the structure divide. (A) L3 ovary. (B) L4 Ovary. (C) Early L5 ovary, note the polyfusomes lined up in each ovariole stained with phalloidin (cyan) (D) Late L5 ovary. Ovaries are orientated with the terminal filament (anterior) to the left. Scale bars are 100 μm. E, F and (G) Clusters of synchronously dividing cells during larval ovary development. See also supplemental videos 1,2 and 3. (E) Larval stage 3 ovary showing a cluster of 4 synchronously dividing cells (arrowhead) (F) L4 stage ovary showing polyfusomes (asterisks) and a cluster of 4 dividing cells (arrowhead). (G) Early L5 ovary showing two ovarioles, both with polyfusomes (asterisks) and a cluster of 8 synchronously dividing cells (arrowhead). Scale bars for E,F and G: 20 μm.

Synchronously dividing cells can be detected during larval ovary development. At larval stages 3 (Fig. 6D and Movie 1) and 4 (Fig. 6E and Movie 2), there are clusters of simultaneously dividing cells in the region of nos positive cells in the ovary. At larval stages 3 and 4, these are clusters of 4 dividing cells. At larval stage 5 (Fig. 6F and Movie 3) these are clusters of eight dividing cells, having undergone duplications. The eight dividing cells, while in the same regions of polyfusomes, have very faint, or difficult-to-see polyfusomes. When the germline cell clusters in adult ovaries divide, the phalloidin staining of the polyfusomes joining them is reduced, making them hard to discern (Cullen et al., 2023). It is unclear why.

The appearance of fusomes, and then polyfusomes, as well as clusters of synchronously dividing cells, implies that 8-cell germline clusters in the adult ovary begin to form in larval stages 2 and 3, with the advent of cells joined by fusomes, which then, in larval stage 4, produce clusters of eight cells joined together by polyfusomes. This clustered organisation of the germline in larval development remains the core structure of the germline in adult queen honeybees.

Pupal growth of the ovary

As noted previously, the pupal ovary does not contain the vitellarium region of the ovary. In adults, the vitellarium contains meiosis and cysts of germ-line nurse cells and oocytes. Oocytes are provisioned and enlarge in the vitellarium, finally producing mature oocytes for laying. To determine when the vitellarium begins to form, we followed the expression of cas, nos and odd through late pupal development and into newly emerged queens as a way to understand the further differentiation of the ovary (Fig. 7).

Fig. 7. Later development of pupal and newly-emerged queen ovaries. Ovaries are stained for DAPI (grey), odd (green), cas (magenta) and nos (cyan). (A) Early pupal ovaries. (B) Late Pupal ovaries. (C) Ovaries from a newly-emerged queen. Ovaries are orientated with the terminal filament (anterior) to the left. Scale bars: 500 μm. Yellow dotted lines show crop lines due to the microscope field of view, and rotation of the specimen to match other samples.

The patterns of gene expression noted during larval development continue into late pupal and newly emerged queen ovaries. RNA expression from odd marks the terminal filaments, with all terminal filament cells expressing odd RNA, until the anterior of the germarium, where clusters of nos expressing germline cell clusters reside. Expression of cas is restricted to somatic cells of both the terminal filament and the germarium. Expression of nos RNA is limited to the germarium.

Even in newly emerged queen ovaries, there is little development of the vitellarium, supporting the idea that the ovary undergoes significant changes in structure and gene expression after mating (Kocher et al., 2008). Expression of cas RNA in newly emerged queens marks a small triangular structure at the end of each germarium (Fig. 7C), perhaps a sign of the beginnings of differentiation of the vitellarium. Cas marks somatic cells (Cullen et al., 2023), such as follicle cells. These ‘triangles’ could signify the production of follicle cells that will envelope developing oocytes (Storto and King, 1989) as they move into the vitellarium.

While the vitellarium does not appear to form until after mating (Kocher et al., 2008), we wanted to know if the division of germline cell clusters is occurring in newly emerged queens. By staining newly emerged queen ovaries with an antibody against PHH3 we can see 8-cell clusters undergoing duplication in the germarium (Fig. S2).

DISCUSSION

Honeybee queens have been reported as laying ∼500-3000 (Allen, 1955; Avni et al., 2014) eggs per day. This egg-laying occurs in bursts, with permissive environmental conditions, and space in the hive required to allow these high rates of reproduction (Fine et al., 2018; Shehata et al., 1981). While the rate of laying is subject to environmental factors, e.g. nutrition and climate (Avni et al., 2014), structures that might affect egg production, such as ovariole number, can differ between queens (Jackson et al., 2011; Tarpy et al., 2000). It is also possible that the number, and quality, of germ-cell clusters can be manipulated during development, alongside ovariole number.

Understanding the development of the larval ovaries in bees required suitable markers. While nos and cas have been used previously, we serendipitously discovered that odd is expressed from early embryonic stages, in the cells of the developing terminal filament. This expression does not occur in Drosophila, and we have, as yet no function for odd in the terminal filament. As a marker of terminal filament cells and development, odd expression demonstrates that, in bees at least, the terminal filament is an early developing component of the ovary; perhaps providing the somatic precursors to each ovariole. The initial ovary is thus a collection of terminal filament cells as somatic precursors to the ovary, and germline cells, all of which form in the embryo. Terminal filament cell number determines ovariole number in Drosophila melanogaster (Sarikaya et al., 2012), and co-varies with ovariole number in other Drosophila species (Green and Extavour, 2012). Given the unusually large number of ovarioles in Apis mellifera, differing from most other Hymenoptera (Hartfelder et al., 2018), odd may provide a key marker in understanding the genetic mechanism underlying the evolutionary increase in ovariole number. Comparing the numbers of odd RNA-expressing cells in late embryogenesis between honeybees and Hymenoptera with smaller numbers of ovarioles and determining how they are specified may help us understand how, and why, queen honeybees produce so many ovarioles.

Tracing the regionalization of the honeybee ovary during larval development indicates that the ovariole structure is determined by larval stage 3, with terminal filament and germline cells residing in different domains by this early stage (Fig. 8). In honeybee testes, the testioles are visible as finger-like structures by the first larval stage (Lago et al., 2020), so the testioles are likely defined in late embryonic stages. Due to the similarities in morphology and structure between honeybee ovarioles and testioles (Hartfelder et al., 2018), ovarioles are likely determined in the embryonic stages.

Fig. 8. Illustration of larval and pupal ovary development in honeybee queens indicating key events in the production of adult ovaries.

Later larval and pupal development involves the expansion and elongation of these domains to produce an ovary in newly emerged queens that is made up of an extensive terminal filament, and a germarium, but with no sign of the vitellarium. It is unclear if meiosis has begun in these ovaries as we have no markers to date that identify cells in meiosis in honeybees.

The key feature of the honeybee queen germline is that germline precursors are organized into 8 cell clusters, joined by a polyfusome (Cullen et al., 2023). These structures reside in the germarium, where they divide at a rate that can replace these germline cell precursors over the years of a queen's life (Cullen et al., 2023). We see evidence of the formation of these clusters beginning in larval stage 2, where we first see fusomes linking two cells in the germline regions of the ovarioles together, and stage 3 where we find synchronously dividing 4-cell clusters and polyfusomes. By early stage 5, there are numerous polyfusomes in the germarium and clusters of 8 cells dividing synchronously (Fig. 8). The germline cell clusters crucial to fertility in the adult queen honeybee thus form during early larval development and are in place in the germarium before the emergence of virgin queens.

The honeybee queen ovary is immature on emergence, as the vitellarium, where oocytes are specified and resourced, has not yet formed. Germline cell clusters are dividing at this time, but it is unclear if these clusters are being stored ready for reproduction to begin, or apoptosed. At this stage, after a few weeks in the hive (Koeniger, 1986), the virgin queens go on a mating flight, mating with multiple males on the wing, and then return to the hive. Egg laying begins 2-3 days (Koeniger, 1986) after mating. Changes in gene expression and morphology of the ovary occur after mating (Kocher et al., 2008), which produces an active vitellarium, and active egg laying.

Queen bees are constrained somewhat by their social biology. They need to be able to lay bursts of eggs when conditions are appropriate, reaching remarkable levels of egg laying (Harbo, 1986), and thus need large and active ovaries. Despite this, virgin queen bees also need to go on a metabolically challenging nuptial flight early in life. It seems likely that the immature ovaries in newly emerged queen bees reflect a trade-off in resource need for flight and reproduction. A similar trade-off occurs at the offset of swarming, as the queen decreases her egg-laying rate (Allen, 1956), thus reducing weight (Seeley and Fell, 1981), to allow for flying with the colony during swarming. By suppressing, or delaying, the development and provisioning of eggs before mating, queens have the resources for an expensive nuptial flight. Trade-offs between reproduction and flight have been identified in many insect species (for examples see Guerra, 2011). After repressing reproduction for a nuptial flight, queens need to rapidly activate their ovaries and begin laying large numbers of eggs to continue populating the hive with workers to replace the nurse workers as they age out of that role to become foragers (Seeley, 1982). The number of ovarioles in bees is unusual, especially for Hymenoptera (Church et al., 2021; Khila and Abouheif, 2010; Kugler et al., 1976; Martins and Serrão, 2004), a testament to the need to rapidly lay large numbers of eggs in ideal situations. Eight-cell clusters provide a ‘timing’ advantage to reproduction, as rather than having to develop and lay eggs from a single germ stem cell, honeybees are further along in the duplication timeline. Perhaps the large number of ovarioles in honeybee queens, and the 8-cell cluster form of the germline in the ovary (Cullen et al., 2023), are adaptations to allow the rapid activation and high rates of reproduction in young honeybee queens.

MATERIALS AND METHODS

Honeybees and staging

Honeybees were sourced from Otto's Bees, a bee-keeping supplier based in Dunedin, Otago, New Zealand. Otto's Bees keeps two stocks of bees, an ‘Italian’ strain, derived from a closed breeding programme maintained by Betta Bees Research Ltd, and a ‘Carniolan’ strain. Replicates of experiments were performed using unrelated hives.

Grafting of embryos into queen cells and queen cell raising was carried out as per Cameron et al. (2013).

Embryos and larvae were staged according to Cridge et al. (2017). We investigated larval stages 2, 3 and 4. The larvae are in stage 5 for 4 days (Cridge et al., 2017) so we investigated ‘early’ and ‘late’ periods. We also split the pupal stage, also 5 days, into early and late stages (Table S1).

The staging was performed by timing individuals from days post-laying, and based on morphology, as described below.

Dissection and fixation

Embryos

Embryos were collected, fixed and dissected using a protocol modified from Osborne and Dearden (2005), as described below. Embryos were collected from worker cells with a paintbrush into 1:1 heptane:4% formaldehyde in 1× PBS. Embryos were shaken overnight at moderate speed on a nutating mixer. After shaking, the lower phase was removed and replaced with 100% methanol at −20°C, this was then shaken vigorously for 3 min by hand. The solution was then removed, and embryos were washed three times in fresh −20°C 100% methanol. Embryos were then stored until dissections in 100% methanol at −20°C.

Embryos and early larva were transferred to a thin-walled glass test tube; 100% methanol was replaced with 2 ml PTw (1× PBS+0.1% Tween 20). Embryos were allowed to settle and then sonicated in a sonic cleaning bath for 5-10 s. Embryos were transferred to a microcentrifuge tube with 200 μg/ml of Pronase (Sigma) in PTw. After 5 min the embryos and larva were washed twice with 500 μl of PTw. The embryos were transferred to a plastic Petri dish with PTw and the outer membranes were dissected away using sharp forceps. Dissected embryos were collected into a microcentrifuge tube with PTw. Embryos were then fixed again in 10% formaldehyde in PBS for 35 min. Embryos were then rinsed six times in 500 μl of PTw.

Collection of larval and pupal stages

Larva and pupa were collected from queen cells. Larval samples were kept in the hive until necessary to continue brood care. Pupae were removed from the hive and kept in an incubator with at least 50% ROH at 34°C. Upon removal from queen cells, larvae and pupae were immediately dissected. Larval stages 2 and 3 were dissected by cutting a slit along the centre of the dorsal midline and loosening the internal organs with tweezers. The entire larva was then stained as below. For Larval stages 4 and 5, ovaries were dissected by cutting a slit along the centre of the dorsal midline and dissecting them from the larval body. For pupal stages, the head and thorax were removed from the abdomen. Then segments three and four were pulled apart and ovaries were removed from the abdomen.

Newly emerged queens were always dissected at less than 12 h old. Newly emerged queens were placed at 4°C until movement ceased, and then the heads were removed from the thorax before immediate dissection. Ovaries were dissected into 1× PBS under a dissection microscope and kept in 1× PBS on ice until fixation (less than 2 h). Once ready for fixation, ovaries were fixed in 1:1 heptane:4% formaldehyde in 1× PBS for 12 (larval stages)-15 (pupal stages and newly emerged queens) minutes, then rinsed in 1× PTx (PBS+0.1% Triton X-100) six times. All ovaries were used immediately in experiments.

In situ hybridization

Embryos

Honeybee embryos were stained via a modification of that used for Nasonia embryos in Taylor and Dearden (2022). Embryos were pre-hybridised in 200 μl of probe hybridisation buffer (2.4 M urea, 5×sodium chloride sodium citrate (SSC), 9 mM citric acid (pH 6.0), 0.1% Tween 20, 50 μg/ml heparin, 1×Denhardt's solution, 10% dextran sulfate) for at least 30 min at 37°C and probe solution was prepared by adding 6μL of each 1 μM probe (Molecular Instruments) to 200 μl probe hybridisation buffer at 37°C. The pre-hybridisation solution was removed, and the probe solution was added. These samples were then incubated for 2 days at 37°C. Unbound probes were removed by 4×15-min wash steps with 200 μl of probe wash buffer [2.4 M Urea, 5×SSC, 9 mM citric acid (pH 6.0), 0.1% Tween, 50 μg/ml heparin] at 37°C and 3×5-min washes with 200 μl 5X SSCT (5X SSC, 0.1% Tween 20) at room temperature. Embryos were pre-amplified in 200 μl of amplification buffer (5X SSC, 0.1% Tween 20, 10% dextran sulfate) for at least 30 min at room temperature. Six μL of h1 and h2 hairpins (Molecular Instruments) were prepared separately, snap-cooled by heating to 95°C for 90 s then left to cool in the dark for 30 min. The hairpin solution was then prepared by adding all the hairpins to 200 μl of amplification buffer at room temperature. The pre-amplification solution was then removed and replaced with the hairpin solution. Embryos were incubated for 2 days in the dark at room temperature. Unbound hairpins were then removed by 2×5-min washes with 200 μl of SSCT, 2×30-min washes with 200 μl of SSCT and 1×5-min wash with 200 μl of SSCT. SSCT was then removed and replaced with 500μL 70% glycerol. Embryos were then bridge-mounted for microscopy.

2-5 Stage Larva/Pupal Stages

Hybridization chain reactions of larva and pupal ovaries were performed as described in Cullen et al. (2023), with probe concentrations of 2 μl of 1 μM probes (nanos, castor, oddskipped). Hairpins were changed between replicates (−488,-546,-647), to ensure no signal bias in specific wavelengths of light.

Immunohistochemistry

Immunohistochemistry, phalloidin and DAPI staining were performed as described in Cullen et al. (2023), with primary antibody for dividing cells [1:200 α-mouse-ph3 (Anti-Histone H3; phospho S10) antibody (abcam 14955)].

Microscopy

Confocal microscopy was performed on an Olympus FV3000 confocal microscope. Ovaries were serially optically sectioned with the number of slices/depth optimized using FV3000 software. Images were produced using Icy (de Chaumont et al., 2012) or FIJI (Schindelin et al., 2012) and images were trimmed for ease of visualisation. All raw data images and media are available at Zenodo (DOI 10.5281/zenodo.10972330).

Phylogenetics

Proteins with similarities to odd skipped-like proteins were identified using blastp (Altschul et al., 1990), and aligned using Clustal (Thompson et al., 1994). Maximum likelihood phylogenetic analysis was carried out using RAxML (Stamatakis, 2014).

Supplementary Material

10.1242/biolopen.060511_sup1 Supplementary information

Acknowledgements

The authors would like to that Dr Otto Hyink for the provision of bees and grafting services and Amanda Austin for honeybee embryo staining. We thank Petra Dearden for critical review of the manuscript.

Footnotes

Author contributions

Conceptualization: G.C., P.K.D.; Methodology: G.C., E.D., P.K.D.; Validation: P.K.D.; Formal analysis: G.C., P.K.D.; Investigation: G.C., P.K.D.; Resources: P.K.D.; Data curation: G.C.; Writing - original draft: G.C.; Visualization: G.C., P.K.D.; Supervision: P.K.D.; Project administration: P.K.D.; Funding acquisition: P.K.D.

Funding

This project was supported by the New Zealand Ministry of Business, Innovation and Employment ‘Selecting Future Bees’ Programme Grant. Open Access funding provided by University of Otago. Deposited in PMC for immediate release.

Data availability

The Authors affirm that raw data files for all figures are available on Zenodo (doi:10.5281/zenodo.10972330).

First person

This article has an associated First Person interview with the first author of the paper.
==== Refs
References

Aamidor, S. E., Cardoso-Júnior, C. A. M., Harianto, J., Nowell, C. J., Cole, L., Oldroyd, B. P. and Ronai, I. (2022). Reproductive plasticity and oogenesis in the queen honey bee (Apis mellifera). J. Insect Physiol. 136 , 104347. 10.1016/j.jinsphys.2021.104347 34902433
Abbo, P. M., Kawasaki, J. K., Hamilton, M., Cook, S. C., DeGrandi-Hoffman, G., Li, W. F., Liu, J. and Chen, Y. P. (2017). Effects of Imidacloprid and Varroa destructor on survival and health of European honey bees, Apis mellifera. Insect Science 24 , 467-477. 10.1111/1744-7917.12335 26990560
Allen, M. D. (1955). Observations on honeybees attending their queen. Br. J. Anim. Behav. 3 , 66-69. 10.1016/S0950-5601(55)80015-9
Allen, M. D. (1956). The behaviour of honeybees preparing to swarm. Br. J. Anim. Behav. 4 , 14-22. 10.1016/S0950-5601(56)80011-7
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215 , 403-410. 10.1016/S0022-2836(05)80360-2 2231712
Avni, D., Hendriksma, H. P., Dag, A., Uni, Z. and Shafir, S. (2014). Nutritional aspects of honey bee-collected pollen and constraints on colony development in the eastern Mediterranean. J. Insect Physiol. 69 , 65-73. 10.1016/j.jinsphys.2014.07.001 25038311
Brosi, B. J., Delaplane, K. S., Boots, M. and de Roode, J. C. (2017). Ecological and evolutionary approaches to managing honeybee disease. Nat. Ecol. Evol. 1 , 1250-1262. 10.1038/s41559-017-0246-z 29046562
Büning, J. (1993). Germ cell cluster formation in insect ovaries. Int. J. Insect Morphol. Embryol. 22 , 237-253. 10.1016/0020-7322(93)90013-Q
Cameron, R. C., Duncan, E. J. and Dearden, P. K. (2013). Biased gene expression in early honeybee larval development. BMC Genomics 14 , 903. 10.1186/1471-2164-14-903 24350621
Chapman, R. F. (1998). The Insects: Structure and Function. Cambridge University Press.
Church, S. H., de Medeiros, B. A. S., Donoughe, S., Márquez Reyes, N. L. and Extavour, C. G. (2021). Repeated loss of variation in insect ovary morphology highlights the role of development in life-history evolution. Proc. R. Soc. B 288 , 20210150. 10.1098/rspb.2021.0150
Coulter, D. E., Swaykus, E. A., Beran-Koehn, M. A., Goldberg, D., Wieschaus, E. and Schedl, P. (1990). Molecular analysis of odd skipped, a zinc finger encoding segmentation gene with a novel pair-rule expression pattern. EMBO J. 9 , 3795-3804. 10.1002/j.1460-2075.1990.tb07593.x 2120051
Cridge, A. G., Lovegrove, M. R., Skelly, J. G., Taylor, S. E., Petersen, G. E. L., Cameron, R. C. and Dearden, P. K. (2017). The honeybee as a model insect for developmental genetics. Genesis 55 , e23019. 10.1002/dvg.23019
Cullen, G., Gilligan, J. B., Guhlin, J. G. and Dearden, P. K. (2023). Germline progenitors and oocyte production in the honeybee queen ovary. Genetics 225 , iyad138. 10.1093/genetics/iyad138 37487025
de Chaumont, F., Dallongeville, S., Chenouard, N., Hervé, N., Pop, S., Provoost, T., Meas-Yedid, V., Pankajakshan, P., Lecomte, T., Le Montegner, Y. et al. (2012). Icy: an open bioimage informatics platform for extended reproducible research. Nat. Methods 9 , 690-696. 10.1038/nmeth.2075 22743774
de Cuevas, M. and Spradling, A. C. (1998). Morphogenesisof the Drosophila fusome and its implications for oocyte specification. Development 125 , 2781-2789. 10.1242/dev.125.15.2781 9655801
Dearden, P. K. (2006). Germ cell development in the Honeybee (Apis mellifera); Vasa and Nanos expression. BMC Dev. Biol. 6 , 6. 10.1186/1471-213X-6-6 16503992
Deng, W. and Lin, H. (2001). Asymmetric germ cell division and oocyte determination during Drosophila oogenesis. Int. Rev. Cytol. 203 , 93-138. 10.1016/S0074-7696(01)03005-4 11131529
Dobens, L. L. and Raftery, L. A. (2000). Integration of epithelial patterning and morphogenesis in Drosophila ovarian follicle cells. Dev. Dyn. 218 , 80-93. 10.1002/(SICI)1097-0177(200005)218:1<80::AID-DVDY7>3.0.CO;2-8 10822261
Fine, J. D., Shpigler, H. Y., Ray, A. M., Beach, N. J., Sankey, A. L., Cash-Ahmed, A., Huang, Z. Y., Astrauskaite, I., Chao, R. and Zhao, H. (2018). Quantifying the effects of pollen nutrition on honey bee queen egg laying with a new laboratory system. PLoS ONE 13 , e0203444.30183759
Gill, R. J. and Hammond, R. L. (2011). Workers influence on royal reproduction. Proc. R. Soc. B 278 , 1524-1531. 10.1098/rspb.2010.1774
Goulson, D., Nicholls, E., Botías, C. and Rotheray, E. L. (2015). Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. Science 347 , 1233957. 10.1126/science.1255957
Green, D. A., II and Extavour, C. G. (2012). Convergent evolution of a reproductive trait through distinct developmental mechanisms in Drosophila. Dev. Biol. 372 , 120-130. 10.1016/j.ydbio.2012.09.014 23022298
Guerra, P. A. (2011). Evaluating the life-history trade-off between dispersal capability and reproduction in wing dimorphic insects: a meta-analysis. Biol. Rev. 86 , 813-835. 10.1111/j.1469-185X.2010.00172.x 21199288
Gutzeit, H. O., Zissler, D. and Fleig, R. (1993). Oogenesis in the honeybee Apis mellifera: Cytological observations on the formation and differentiation of previtellogenic ovarian follicles. Roux's Arch. Dev. Biol. 202 , 181-191. 10.1007/BF00365309 28305996
Hans, F. and Dimitrov, S. (2001). Histone H3 phosphorylation and cell division. Oncogene 20 , 3021-3027. 10.1038/sj.onc.1204326 11420717
Harbo, J. R. (1986). Effect of population size on brood production, worker survival and honey gain in colonies of Honeybees. J. Apic. Res. 25 , 22-29. 10.1080/00218839.1986.11100687
Hart, M. C., Wang, L. and Coulter, D. E. (1996). Comparison of the structure and expression of odd skipped and two related genes that encode a new family of zinc finger proteins in Drosophila. Genetics 144 , 171-182. 10.1093/genetics/144.1.171 8878683
Hartfelder, K. and Steinbrück, G. (1997). Germ cell cluster formation and cell death are alternatives in caste-specific differentiation of the larval honey bee ovary. Invertebr. Reprod. Dev. 31 , 237-250. 10.1080/07924259.1997.9672582
Hartfelder, K., Tiberio, G. J., Lago, D. C., Dallacqua, R. P. and Bitondi, M. M. G. (2018). The ovary and its genes - developmental processes underlying the establishment and function of a highly divergent reproductive system in the female castes of the honey bee, Apis mellifera. Apidologie 49 , 49-70. 10.1007/s13592-017-0548-9
Jackson, K. and Robinson, G. E. (2018). Contest experience does not increase survivorship in honey bee queen duels. Insectes Soc. 65 , 631-637. 10.1007/s00040-018-0654-3
Jackson, J. T., Tarpy, D. R. and Fahrbach, S. E. (2011). Histological estimates of ovariole number in honey bee queens, Apis mellifera, reveal lack of correlation with other queen quality measures. J. Insect Sci. 11 , 82. 10.1673/031.011.8201 21870968
Kang, Y., Blanco, K., Davis, T., Wang, Y. and DeGrandi-Hoffman, G. (2016). Disease dynamics of honeybees with Varroa destructor as parasite and virus vector. Math. Biosci. 275 , 71-92. 10.1016/j.mbs.2016.02.012 26968931
Khalifa, S. A. M., Elshafiey, E. H., Shetaia, A. A., El-Wahed, A. A. A., Algethami, A. F., Musharraf, S. G., AlAjmi, M. F., Zhao, C., Masry, S. H. D., Abdel-Daim, M. M. et al. (2021). Overview of bee pollination and its economic value for crop production. Insects 12 , 688. 10.3390/insects12080688 34442255
Khila, A. and Abouheif, E. (2010). Evaluating the role of reproductive constraints in ant social evolution. Philos. Trans. R. Soc. B Biol. Sci. 365 , 617-630. 10.1098/rstb.2009.0257
King, R. C. (1979). Aberrant fusomes in the ovarian cystocytes of the FS (1) 231 mutant of Drosophila melanogaster Meigen (Diptera: Drosophilidae). Int. J. Insect Morphol. Embryol. 8 , 297-309. 10.1016/0020-7322(79)90038-2
Kocher, S. D., Richard, F.-J., Tarpy, D. R. and Grozinger, C. M. (2008). Genomic analysis of post-mating changes in the honey bee queen (Apis mellifera). BMC Genomics 9 , 1-15. 10.1186/1471-2164-9-232 18171476
Koeniger, G. (1986). Reproduction and mating behavior. In Bee Genetics and Breeding, pp. 255-280. Louisiana, USA: Academic Press.
Kozii, I. V., Wood, S. C., Koziy, R. V. and Simko, E. (2022). Histomorphological description of the reproductive system in mated honey bee queens. J. Apic. Res. 61 , 114-126. 10.1080/00218839.2021.1900636
Kugler, J., Orion, T. and Ishay, J. (1976). The number of ovarioles in the Vespinae (Hymenoptera). Insectes Soc. 23 , 525-533. 10.1007/BF02223479
Lago, D. C., Martins, J. R., Dallacqua, R. P., Santos, D. E., Bitondi, M. M. G. and Hartfelder, K. (2020). Testis development and spermatogenesis in drones of the honey bee, Apis mellifera L. Apidologie 51 , 935-955. 10.1007/s13592-020-00773-2
Manzoor, F. and Pervez, M. (2022). Pesticide impact on honeybees declines and emerging food security crisis. In Hamadttu Abdel Farag El-Shafie (eds) Global Decline of Insects [Internet]. Rijeka: IntechOpen; 2021 p. Ch. 4. Available from: 10.5772/intechopen.98871
Martins, G. F. and Serrão, J. E. (2004). A comparative study of the ovaries in some Brazilian bees (Hymenoptera; Apoidea). Papéis Avulsos De Zoologia 44 , 45-53. 10.1590/S0031-10492004000300001
Martins, J. R., Anhezini, L., Dallacqua, R. P., Simões, Z. L. P. and Bitondi, M. M. G. (2011). A honey bee hexamerin, HEX 70a, is likely to play an intranuclear role in developing and mature ovarioles and testioles. PLoS ONE 6 , e29006. 10.1371/journal.pone.0029006 22205988
Osborne, P. and Dearden, P. K. (2005). Non-radioactive in-situ hybridisation to honeybee embryos and ovaries. Apidologie 36 , 113-118. 10.1051/apido:2004075
Plate, M., Bernstein, R., Hoppe, A. and Bienefeld, K. (2019). The importance of controlled mating in honeybee breeding. Genet. Sel. Evol. 51 , 74. 10.1186/s12711-019-0518-y 31830903
Pritsch, M. and Büning, J. (1989). Germ cell cluster in the panoistic ovary of Thysanoptera (Insecta). Zoomorphology 108 , 309-313. 10.1007/BF00312163
Rachinsky, A., Strambi, C., Strambi, A. and Hartfelder, K. (1990). Caste and metamorphosis: Hemolymph titers of juvenile hormone and ecdysteroids in last instar honeybee larvae. Gen. Comp. Endocrinol. 79 , 31-38. 10.1016/0016-6480(90)90085-Z 2354779
Ratnieks, F. L. W. and Halenterä, H. (2009). The evolution of extreme altruism and inequality in insect societies. Philos. Trans. R. Soc. B Biol. Sci. 364 , 3169-3179. 10.1098/rstb.2009.0129
Sarikaya, D. P., Belay, A. A., Ahua, A., Dorta, A., Green, D. A., II and Extavour, C. G. (2012). The roles of cell size and cell number in determining ovariole number in Drosophila. Dev. Biol. 363 , 279-289. 10.1016/j.ydbio.2011.12.017 22200592
Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., Preibisch, S., Rueden, C., Saalfeld, S., Schmid, B. et al. (2012). Fiji: an open-source platform for biological-image analysis. Nat. Methods 9 , 676-682. 10.1038/nmeth.2019 22743772
Seeley, T. D. (1982). Adaptive signficance of the age polyethism schedule in honeybee colonies. Behav. Ecol. Sociobiol. 11 , 287-293. 10.1007/BF00299306
Seeley, T. D. and Fell, R. D. (1981). Queen substance production in Honey bee (Apis mellifera) colonies preparing to swarm (Hymenoptera: Apidae). J. Kans. Entomol. Soc. 54 , 192-196.
Shehata, S. M., Townsend, G. F. and Shuel, R. W. (1981). Seasonal physiological changes in queen and worker honeybees. J. Apic. Res. 20 , 69-78. 10.1080/00218839.1981.11100475
Snodgrass, R. E. (1956). Anatomy of the Honey Bee. Cornell University Press.
Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30 , 1312-1313. 10.1093/bioinformatics/btu033 24451623
Starks, P. T., Blackie, C. A. and Seeley, T. D. (2000). Fever in honeybee colonies. Naturwissenschaften 87 , 229-231. 10.1007/s001140050709 10883439
St Johnston, D. (2008). Drosophila oogenesis. Curr. Biol. 18 , R1082-R1087. 10.1016/j.cub.2008.09.011 19081037
Storto, P. D. and King, R. C. (1989). The role of polyfusomes in generating branched chains of cystocytes during Drosophila oogenesis. Dev. Genet. 10 , 70-86. 10.1002/dvg.1020100203 2499436
Stuligross, C. and Williams, N. M. (2021). Past insecticide exposure reduces bee reproduction and population growth rate. Proc. Natl. Acad. Sci. U.S.A. 118 , e2109909118. 10.1073/pnas.2109909118 34810261
Tanaka, E. D. and Hartfelder, K. (2004). The initial stages of oogenesis and their relation to differential fertility in the honey bee (Apis mellifera) castes. Arthropod. Struct. Dev. 33 , 431-442. 10.1016/j.asd.2004.06.006 18089049
Tarpy, D. R., Hatch, S. and Fletcher, D. J. C. (2000). The influence of queen age and quality during queen replacement in honeybee colonies. Anim. Behav. 59 , 97-101. 10.1006/anbe.1999.1311 10640371
Taylor, S. E. and Dearden, P. K. (2022). The Nasonia pair-rule gene regulatory network retains its function over 300 million years of evolution. Development 149 , dev199632. 10.1242/dev.199632 35142336
Thompson, J. D., Higgins, D. G. and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22 , 4673-4680. 10.1093/nar/22.22.4673 7984417
Trauner, J. and Büning, J. (2007). Germ-cell cluster formation in the telotrophic meroistic ovary of Tribolium castaneum (Coleoptera, Polyphaga, Tenebrionidae) and its implication on insect phylogeny. Dev. Genes Evol. 217 , 13-27. 10.1007/s00427-006-0114-3 17123126
Ventelä, S. (2006). Cytoplasmic bridges as cell-cell channels of germ cells. In F. Baluska, D. Volkmann and P. W. Barlow (eds). Cell-Cell Channels. New York: Springer, pp. 208-216.
Visscher, P. K. and Dukas, R. (1995). Honey bees recognize development of nestmates’ ovaries. Anim. Behav. 49 , 542-544. 10.1006/anbe.1995.0074
Ward, E. J. and Coulter, D. E. (2000). Odd skipped is expressed in multiple tissues during Drosophila embryogenesis. Mech. Dev. 96 , 233-236. 10.1016/S0925-4773(00)00389-0 10960790
Warner, S., Pokhrel, L. R., Akula, S. M., Ubah, C. S., Richards, S. L., Jensen, H. and Kearney, G. D. (2024). A scoping review on the effects of Varroa mite (Varroa destructor) and global honey bee decline. Sci. Total Environ. 906 , 167492. 10.1016/j.scitotenv.2023.167492 37778563
Wilson, M. J. and Dearden, P. K. (2012). Pair-rule gene orthologues have unexpected maternal roles in the honeybee (Apis mellifera). PLOS ONE 7 , e46490. 10.1371/journal.pone.0046490 23029534
Wilson, E. O. and Hölldobler, B. (2005). Eusociality: origin and consequence. Proc. Natl Acad. Sci. USA 102 , 13367-13371. 10.1073/pnas.0505858102 16157878
Woyke, J. (1963). Drone larvae from fertilized eggs of the honeybee. J. Apic. Res. 2 , 19-24. 10.1080/00218839.1963.11100052
Woyke, J. (1964). Causes of repeated mating flights by queen honeybees. J. Apic. Res. 3 , 17-23. 10.1080/00218839.1964.11100077
Yang, Y., Wu, Y., Long, H., Ma, X., Shariati, K., Webb, J., Guo, L., Pan, Y., Ma, M., Deng, C. et al. (2023). Global honeybee health decline factors and potential conservation techniques. Food Sec. 15 , 855-875. 10.1007/s12571-023-01346-8
