Juvenile hormone is considered to be a master regulator of polyphenism in social insects. In the ant Cardiocondyla obscurior, whether a female egg develops into a queen or a worker is determined maternally and caste-specific differentiation occurs in embryos, so that queens and workers can be distinguished in a non-invasive manner from late embryogenesis onwards. This ant also exhibits two male morphs – winged and wingless males. Here, we used topical treatment with juvenile hormone III and its synthetic analogue methoprene, a method that influences caste determination and differentiation in some ant species, to investigate whether hormone manipulation affects the development and growth of male, queen- and worker-destined embryos and larvae. We found no effect of hormone treatment on female caste ratios or body sizes in any of the treated stages, even though individuals reacted to heightened hormone availability with increased expression of krüppel-homolog 1, a conserved JH first-response gene. In contrast, hormone treatment resulted in the emergence of significantly larger males, although male morph fate was not affected. These results show that in C. obscurior, maternal caste determination leads to irreversible and highly canalized caste-specific development and growth.

Queen-worker caste polyphenism in social Hymenoptera (ants, bees, wasps) is a prime example of developmental plasticity. In ants, various factors can be associated with caste-specific development, including genotype (Helms-Cahan et al., 2002; Julian et al., 2002; Pearcy et al., 2004), maternal effects (Bier, 1952; Passera, 1980; Schwander et al., 2008), and the abiotic and biotic environment (e.g. Wesson, 1940; Brian, 1954, 1955, 1973; Vargo and Passera, 1992; Penick and Liebig, 2012). To explain how intrinsic and extrinsic factors are translated into differential caste development, Wheeler proposed juvenile hormone (JH) as a master regulator (Wheeler, 1986), inspired by work on other insects showing that JH, together with ecdysone, controls the duration of development and growth via the timing of moults (Nijhout, 2003), thus making it an important factor in regulating size polymorphism (Nijhout and Wheeler, 1982). For ants, this insight was based on two studies of Pheidole spp. soldier caste determination (Wheeler and Nijhout, 1981, 1983) and five studies in three species (Pheidole pallidula, Myrmica rubra, Solenopsis invicta) demonstrating that treatment with juvenile hormone or its synthetic analogues could result in increased queen production (Brian, 1974; Vinson and Robeau, 1974; Robeau and Vinson, 1976; Banks et al., 1978; Passera and Suzzoni, 1979).

Since Wheeler's seminal work, to the best of our knowledge, a mere five additional studies have investigated the relationship between hormone signalling and queen/worker development in ants (Schrempf and Heinze, 2006; Helms-Cahan et al., 2011; Penick et al., 2012; Libbrecht et al., 2013; Kuhn et al., 2018), and only one has followed up on previous reports (de Menten et al., 2005), so that empirical evidence for the generality of JH-mediated regulation of caste and caste-specific body size is still surprisingly limited. One reason for this is the difficulty of studying development in ants: many do not mate in the lab or only produce sexuals when colonies are sufficiently large and, unlike honey bees, ants do not rear castes in distinct brood cells, making it challenging to follow individual development and to identify castes during early development (Schultner and Pulliainen, 2020).

The genus Cardiocondyla is one of 10 ant genera in which complete worker sterility has evolved, meaning that adult workers lack reproductive organs. We recently showed that in Cardiocondyla obscurior, queen- and worker-destined embryos and larvae can be distinguished in a non-invasive manner and with high accuracy using caste-specific crystalline deposits (Schultner et al., 2023). Along with differences in the presence of crystalline deposits, embryos of the two castes are distinct in their miRNA and gene expression (Oettler et al., 2023 preprint). Final (third) instar larvae also exhibit caste-specific gene expression (Schrader et al., 2015; Klein et al., 2016) and wing disc morphology (Oettler et al., 2019), and castes differ in size from the second larval instar onwards (Schultner et al., 2023). Single-queen colonies produce more queen-destined offspring as queens age (Jaimes-Niño et al., 2022) and workers do not engage in differential rearing or culling according to caste (Schultner et al., 2023). Together, this indicates that female caste in this species is determined very early in development, likely via maternal effects. In addition to the two female castes, C. obscurior produces two discrete male morphs – large, winged males and small, wingless males (Kinomura and Yamauchi, 1987; Heinze and Hölldobler, 1993; Oettler et al., 2010). Like female castes, the two male morphs differ in wing disc development (Oettler et al., 2019) and gene expression in the final larval instar (Schrader et al., 2015; Klein et al., 2016). However, male morphs cannot be distinguished visually before this stage and when male morph is determined is unknown, although it is presumed to occur later than in females (Schrempf and Heinze, 2006).

Using this system, we tested (1) how a maternal mode of female caste determination, associated with very early morphological and molecular differentiation, affects the action of JH and the degree of flexibility in caste-specific development and (2) whether male and female polyphenism are regulated in a similar manner. We did this by subjecting embryos and larvae to (synthetic) juvenile hormone treatment, which has previously been linked to increased queen and winged male production in C. obscurior (Schrempf and Heinze, 2006). We documented effects on adult caste and morph ratios and on larva and adult body sizes, as well as on the expression of the conserved JH-response gene krüppel-homolog 1.

Study species

Cardiocondyla obscurior Wheeler 1929 is a myrmicine tramp ant found in the tropics and subtropics. These small ants live in colonies with a few dozen workers and several queens (Heinze, 2017; Oettler, 2020). Queens, workers and males develop via an egg stage and three larval instars (first instar: L1; second instar: L2, third instar: L3), which can be identified by body shape and melanization of the mandibles (Schrempf and Heinze, 2006). Queen and worker-destined eggs and larvae can be distinguished using caste-specific crystalline deposit patterns (Schultner et al., 2023). All individuals used in this study derived from laboratory stock colonies originally collected in Japan in 2010 and in Brazil in 2009 (Schrader et al., 2014; Errbii et al., 2021; Ün et al., 2021b) and were treated according to 3R principles (https://www.nc3rs.org.uk/). Stock and experimental colonies were kept in 9.6×9.6×3 cm plastic boxes with plastered bottoms. In each colony, a plastic insert nest with three chambers and a microscope slide covered with black foil as a lid was placed on a notch in the plaster to prevent the insert nest from slipping and consequently damaging the ants and brood. To keep workers from escaping, the upper third of the nest was paraffined weekly. Three metal plates were placed inside each colony holding sponges for water supply, honey and insect prey. All colonies were kept in a climate chamber under a 12 h:12 h and 22°C:26°C light:dark cycle at 75% relative humidity and fed three times a week with chopped cockroaches (once per week) or fruit flies (twice per week), honey and water.

Hormone treatment of larvae of unknown caste and sex

In order to replicate the experimental set-up of a previous study (Schrempf and Heinze, 2006), we treated first, second and third instar larvae of unknown sex and caste with 2 µl of either methoprene (Sigma-Aldrich; 1 mg ml−1 diluted in 70% ethanol) (L1: n=157, L2: n=95, L3: n=40), juvenile hormone III (JH III) (Sigma-Aldrich; 1 mg ml−1 diluted in 70% ethanol) (L1: n=129, L2: n=36, L3: n=24), 70% ethanol (L1: n=224, L2: n=102, L3: n=21), or acetone (L1: n=66, L2: n=107, L3: n=23) using a Hamilton syringe. We also included a handling control (L1: n=468, L2: n=223, L3: n=30). Individuals were collected from several stock colonies and pooled prior to treatment to remove any colony-level effects. Individuals were treated in groups of five and then transferred on filter paper to rearing nests containing workers from stock colonies. We chose to treat brood in groups because the small size of C. obscurior makes individual treatment of eggs and larvae difficult. For practical reasons, each group received the same amount (2 µg) of JH III or methoprene, which represents doses of ∼20,000 µg g−1 for L1 larvae, 10,000 µg g−1 for L2 larvae and 1000 µg g−1 for L3 larvae; these doses are far higher than the effective doses used by Wheeler and Nijhout (1981) to induce soldier development in Pheidole (50–250 µg g−1 body weight), but about four times lower than those used by Schrempf and Heinze (2006). Rearing colonies contained between 10 and 30 workers (∼1:2 worker:brood ratio) and worker numbers were standardized weekly until no more brood remained. The sex and caste/morph of each emerging individual was documented and survival and caste/morph ratios were calculated. Caste/morph ratios were compared between treatments using Fisher's exact tests in R v.4.2.2 (https://www.r-project.org/).

Hormone treatment of queen- and worker-destined late-stage embryos and larvae

To examine whether juvenile hormone can influence caste-specific growth and development after caste has been determined, queen- and worker-destined eggs and L1, L2, L3 larvae were treated with either methoprene, JH III, ethanol (solvent control) or left untreated (handling control). For better identification of caste-specific crystalline patterns in eggs, these were selected in PBT solution (0.3%). For each developmental stage and caste, individuals were collected from several stock colonies, pooled and then randomly divided into groups of five on filter paper to remove any colony-level effects. Each group of five individuals was topically treated with either 1 μl methoprene (1 mg ml−1 diluted in 70% ethanol), 1 μl JH III (1 mg ml−1 diluted in 70% ethanol), 1 μl of 70% ethanol, or none of the above, using a Hamilton syringe. We reduced doses compared with the previous experiment to minimize mortality and because we found no effect even with higher doses (see next paragraph). This process was conducted twice per treatment day, resulting in 10 individuals being treated per day, treatment, developmental stage and caste. Following treatment, the 10 individuals were transferred on filter paper to a nest containing 15 adult workers collected from several stock colonies the previous day. Two hours after treatment, two individuals were randomly selected from each experimental colony, individually flash-frozen in liquid nitrogen and stored in −80°C for gene expression analyses (see below). The remaining individuals remained in the experimental nests to be reared to adulthood. A total of 2240 individuals were treated, of which a subset were sampled for gene expression analyses and the rest were allowed to develop until adulthood (n=56 per treatment, developmental stage and caste).

To assess the effect of hormone dose, in a separate experiment we treated L2 queen- and worker-destined larvae with 2 µl of either methoprene (1 mg ml−1 diluted in 70% ethanol) (queen-destined: n=10, worker-destined: n=33), JH III (1 mg ml−1 diluted in 70% ethanol) (queen-destined: n=15, worker-destined: n=15), 70% ethanol (queen-destined: n=11, worker-destined: n=15) or acetone (queen-destined: n=15, worker-destined: n=20) using a Hamilton syringe. Individuals were collected from several stock colonies, pooled prior to treatment to remove any colony-level effects and treated in groups of five larvae.

Caste and morph ratios

Experimental colonies were observed once a week to record brood survival and developmental stage. Additionally, the number of workers was counted and any dead or missing workers were replaced from stock colonies. Once the first brood of a colony completed pupation, colonies were checked every other day and when melanization set in, colonies were observed daily to ensure that newly emerged, light-coloured adult workers were isolated before they darkened, so as to remain distinguishable from rearing workers. All freshly emerged adults from one experimental colony were transferred together to a smaller round nest with a plastered bottom containing sponges for water supply, a metal plate with honey, and a folded black foil as a shelter, where they were kept for 5–7 days to allow cuticle hardening. The caste and sex of each fully melanized adult was documented and all adults were frozen in Eppendorf tubes at −20°C for morphometric measurements. From these data, we calculated survival and caste prediction accuracy for each treatment, developmental stage and caste, and compared these using Fisher's exact tests in R (Table 1). Because worker- and male-destined eggs and larvae exhibit similar crystalline deposit patterns (E.S., pers. obs.), in each treatment a small proportion of individuals developed into males (31 males in total, Table 1). These were excluded prior to survival and caste prediction accuracy analyses. Note that the survival and caste prediction accuracy results from the handling control have been published as supplemental data in a previous study (Schultner et al., 2023); we nevertheless include them here for comparison.

Table 1.

Survival and caste ratios of queen- and worker-destined embryos and larvae in the ant Cardiocondyla obscurior after hormone treatment

Survival and caste ratios of queen- and worker-destined embryos and larvae in the ant Cardiocondyla obscurior after hormone treatment
Survival and caste ratios of queen- and worker-destined embryos and larvae in the ant Cardiocondyla obscurior after hormone treatment

Development time

Development time was defined as the time in days between treatment and emergence of the first adult from each experimental colony. Development time was calculated separately for each treatment, developmental stage and caste. In cases where more than one caste emerged from an experimental colony, development time was calculated for the first hatched individual of each caste. In late-stage embryos, L1 and L2 larvae, development time was compared between castes and treatments using a linear regression with caste and treatment as explanatory variables, as well as an interaction term to account for caste-specific responses to treatment: development time∼caste×treatment. Model fit was assessed using residual tests implemented in the DHarma package in R (https://cran.r-project.org/package=DHARMa). Tukey-corrected pairwise P-values were calculated using the emmeans package in R (https://CRAN.R-project.org/package=emmeans). In L3 larvae, Kruskal–Wallis tests (function kruskal.test in R v. 4.2.2) were used to test the effect of caste and treatment on development time separately, as the data did not follow a normal distribution.

Body size measurements

To assess the effect of hormone treatment on adult body size, all adults which emerged in experimental colonies were dried under a stereomicroscope light and pinned for subsequent measurements (Table 1). Pinned individuals were photographed at 200× using a Keyence stereomicroscope connected to a camera (Keyence VHX). Head width (HW), head length (HL), mesosoma width (MW), mesosoma length (ML) and petiole width (PW) were measured from photographs after appropriate scaling using ImageJ. In the cases in which adult female caste did not align with presumed caste based on crystalline deposit patterns [queens: 12.8% (51/397) of individuals; workers: 3% (10/331) of individuals; Table 1], the emerging adults were nonetheless measured and included in subsequent analyses. Similarly, all emerged males were measured and included in analyses. The final data set contained 759 individuals (queens – late-stage embryos: handling control, n=9; solvent control, n=13; JH III, n=8; methoprene, n=13; L1 larvae: handling control, n=14; solvent control, n=23; JH III, n=18; methoprene, n=15; L2 larvae: handling control, n=26; solvent control, n=25; JH III, n=28; methoprene, n=25; L3 larvae: handling control, n=43; solvent control, n=38; JH III, n=27; methoprene, n=31; workers – late-stage embryos: handling control, n=27; solvent control, n=16; JH III, n=20; methoprene, n=20; L1 larvae: handling control, n=24; solvent control, n=19; JH III, n=19; methoprene, n=18; L2 larvae: handling control, n=29; solvent control, n=27; JH III, n=27; methoprene, n=29; L3 larvae: handling control, n=23; solvent control, n=24; JH III, n=29; methoprene, n=21; ergatoid males – late-stage embryos: handling control, n=1; solvent control, n=1; JH III, n=4; L1 larvae: handling control, n=1; solvent control, n=2; JH III, n=1; methoprene, n=3; L2 larvae: handling control, n=1; solvent control, n=2; JH III, n=1; methoprene, n=1; L3 larvae: handling control, n=2; solvent control, n=5; JH III, n=1; methoprene: n=5). A principal component analysis (PCA) of the five measured traits was done using the NIPALS algorithm implemented in the pcaMethods package in R v4.2.2 (Stacklies et al., 2007) to estimate missing values. Principal component 1 scores were compared between the three castes (queens, workers, males), the four developmental stages (late-stage embryos, L1, L2, L3 larvae) and between treatments using a linear regression, with caste, treatment and developmental stage as explanatory variables and an interaction term to account for caste-specific responses to treatment: PC score∼caste×treatment+developmental stage. Residual tests implemented in the DHarma package in R v4.2.2 (Hartig, 2020) were used to assess model fit. There was no significant effect of developmental stage on adult body size (Table 2), so this term was removed from the model before pairwise P-values were calculated and Tukey-corrected for multiple comparisons using the emmeans package in R (Lenth et al., 2020).

Table 2.

Effects of hormone treatment on body size and development time in the ant Cardiocondyla obscurior

Effects of hormone treatment on body size and development time in the ant Cardiocondyla obscurior
Effects of hormone treatment on body size and development time in the ant Cardiocondyla obscurior

Since the treatment was applied in a single topical application in each developmental stage, it is possible that an effect on body size was diluted over time. To assess the short-term effect of hormonal manipulation, we repeated hormone treatments using only L2 larvae, which were measured once they reached the L3 stage. As in the previous experiment, L2 larvae of each caste were chosen based on crystalline deposit patterns and treated in groups of five with 1 µl of either JH III, methoprene or ethanol as a solvent control, or subjected to a handling control. The five larvae were transferred to a nest with 12 adult workers originating from several stock colonies which had been set up one day earlier. This procedure was repeated for a total of 35 L2 larvae per treatment and caste. After treatment, nests were checked every day for moulted larvae and all larvae that had reached the L3 stage were anaesthetized under CO2 for 4 h. Larvae were then placed on their backs and photographed at 200x using a Keyence stereomicroscope connected to a camera (Keyence VHX). The order in which the pictures were taken was randomized to buffer the effect of time as larvae were more likely to die, lose water (change of body volume) or wake up after a long waiting period. Larval head capsule width, body width and body length were measured after appropriate scaling using ImageJ (v. 1.53 t). As some larvae died or were damaged in the pre-imaging process, the final data set consisted of 29–35 individuals per caste and treatment (Queens - handling control: n= 32, solvent control: n=29, JH III: n=35, methoprene: n=32; Workers - handling control: n=32, solvent control: n=32, JH III: n=30, methoprene: n=28). Like in adults, a principal component analysis including all traits was run using the NIPALS algorithm implemented in the pcaMethods package in R v4.2.2 (Stacklies et al., 2007) to estimate missing values. Principal component 1 scores were then compared between castes and between treatments using a linear regression, with caste and treatment as explanatory variables and an interaction term to account for caste-specific responses to treatment: PC 1 score∼caste×treatment. Residual tests implemented in the DHarma package in R v4.2.2 (Hartig, 2020) were used to test model fit. Pairwise P-values were calculated and Tukey-corrected for multiple comparisons using the emmeans package in R v4.2.2 (Lenth et al., 2020).

Hormone treatment of early-stage embryos

To test whether caste/morph determination and differentiation can be manipulated early in development before caste-specific crystalline deposits are visible, we treated female and male eggs which were at most 24 h old. To obtain female eggs, experimental colonies containing three mated queens of unknown age, 10 adult workers and five worker pupae were set up from stock colonies (n=25 colonies). To obtain male eggs, experimental colonies containing three queen pupae prior to mating, 15 adult workers and five worker pupae were set up from stock colonies (n=40 colonies). Ants exhibit haplodiploid sex determination, so unmated queens only produce male-destined eggs. Experimental colonies were monitored once per week until the first eggs were observed. Once experimental colonies began producing sufficient numbers of eggs, eggs were removed daily to ensure that only eggs which were at most 24 h old were included in experiments. Eggs from several experimental colonies were pooled to remove any colony-level effects and then divided into groups of five on filter paper, after which they were topically treated with either 1 μl of methoprene (1 mg ml−1 diluted in 70% ethanol) or 1 μl of 70% ethanol using a Hamilton syringe. This was done in parallel for both sexes to avoid bias. Treated eggs were transferred on filter paper to rearing colonies containing 15 workers randomly collected from stock colonies; rearing colonies were always set up one day before egg transfer to allow workers to acclimate. Each rearing colony contained 5-25 eggs, as the number of treated eggs depended on the fecundity of queens, which is much higher in mated compared to unmated queens and varies over time (E.S., pers. obs.). In total, we treated 650 early-stage male embryos (solvent control: n=325, methoprene: n=325) and 570 early-stage female embryos (solvent control: n=285, methoprene: n=285). Rearing colonies were monitored twice per week and fed in the same manner as stock colonies. Freshly emerged adults were transferred to a smaller nest with a plastered bottom and sponges for water supply, a metal plate with honey, and a folded black foil as a shelter, and kept for 3-5 days to allow cuticle hardening. The caste/morph of each fully melanized adult was documented and all adults were frozen in Eppendorf tubes at −20°C for morphometric measurements. From these data, survival and caste/morph ratios were calculated and survival and treatment effects were analysed using Fisher's exact tests in R v4.2.2.

Body size measurements of adults emerging from treatment of early-stage embryos were obtained as described above for later developmental stages. The final data set consisted of 74 individuals (queens - solvent control: n=1, methoprene: n=4; workers - solvent control: n=24, methoprene: n=24; ergatoid males - solvent control: n=10, methoprene: n=11). A principal component analysis using the NIPALS algorithm to estimate missing values implemented in the pcaMethods package in R (Stacklies et al., 2007) was run on all traits. Principal component 1 and 2 scores were compared between castes and treatments using a linear regression, with caste and treatment as explanatory variables, as well as an interaction term to account for caste-specific responses to treatment: PC 1/2 score∼ caste×treatment. Residual tests implemented in the DHarma package in R (https://cran.r-project.org/package=DHARMa) were used to test model fit and pairwise P-values were calculated and Tukey-corrected for multiple comparisons using the emmeans package in R (https://CRAN.R-project.org/package=emmeans).

Gene expression analysis

qPCR was used to assess the effect of sex, caste and treatment on the expression of the JH-response gene Krüppel-homolog 1 (kr-h1). C. obscurior sequences for kr-h1 were retrieved by running a BLAST search of kr-h1 protein sequences obtained for Drosophila melanogaster (CG45074, www.flybase.org) and Apis mellifera (GB45427, www.hymenopteragenome.org) against the C. obscurior 1.4 genome (Schrader et al., 2014). Intron-spanning primers for the extracted C. obscurior ortholog (Cobs_15554) were designed in Geneious Prime vR10 (https://www.geneious.com). Specific primers (kr-h1 forward primer sequence: CTTGGTGTGCAGCCCGGACC; kr-h1 reverse primer sequence: ACCGGTACGGATCCTCGCCC) were tested on pooled samples of queen, worker, ergatoid male and winged male DNA in a temperature gradient PCR (60°C, 63°C, 66°C) and PCR products subsequently sequenced to confirm amplification of the desired transcript. Primer efficiency was tested in a five-step dilution series and was confirmed to be within the acceptable range.

To establish a timeline of kr-h1 expression and to investigate sex-specific responses to JH, second instar male and female larvae (of unknown caste) were collected from several stock colonies, pooled by sex and then treated in groups of five with 2 µl of methoprene (1 mg ml−1 diluted in 70% ethanol, female larvae only), JH III (1 mg ml−1 diluted in 70% ethanol) or 70% ethanol using a Hamilton syringe. We further included a handling control, which consisted of larvae which were removed from colonies and immediately flash-frozen in liquid nitrogen and stored at −80°C. Female larvae were collected from stock colonies containing mated queens whereas male larvae were collected from colonies containing only unmated queens. After treatment, larvae were moved to rearing nests containing workers from stock colonies, which had been set up the day before. One, two, six and 24 h after treatment, treated larvae were removed from rearing colonies, individually flash-frozen in liquid nitrogen and stored in −80°C. We ensured that only larvae which had been accepted by workers, i.e. which had been moved inside the nest after treatment, were included. The final sample consisted of 177 females (handling control: 0 h, n=19; ethanol: 1 h, n=14; 2 h, n=11; 6 h, n=13; 24 h, n=10; JH III: 1 h, n=12; 2 h, n=12; 6 h, n=15; 24 h, n=15; methoprene: 1 h, n=16; 2 h, n=9; 6 h, n=13; 24 h, n=18) and 111 males (handling control: 0 h, n=21; ethanol: 1 h, n=13; 2 h, n=13; 6 h, n=8; 24 h, n=6; JH III: 1 h, n=13; 2 h, n=12; 6 h, n=6; 24 h, n=19).

To assess caste-specific responses to hormone treatment, kr-h1 expression was measured in queen- and worker-destined late-stage embryos, L1, L2 and L3 larvae. All samples were collected from the experiment described in ‘Hormone manipulations of queen- and worker-destined late-stage embryos and larvae’, 2 h after treatment. This timepoint was chosen based on the timeline of kr-h1 expression in female embryos of unknown caste described above (Fig. 3). The final data set consisted of 5–14 individuals per treatment, developmental stage and caste (Queens – late-stage embryos: handling control, n=14; solvent control, n=9; JH, n=12; methoprene, n=10; L1 larvae: handling control, n=5; solvent control, n=8; JH, n=10; methoprene, n=10; L2 larvae: handling control, n=5; solvent control, n=8; JH, n=8; methoprene, n=8; L3 larvae: handling control, n=7; solvent control, n=9; JH, n=9; methoprene, n=9; Workers – late-stage embryos: handling control, n=12; solvent control, n=10; JH, n=9; methoprene, n=9; L1 larvae: handling control, n=6; solvent control, n=7; JH, n=7; methoprene, n=8; L2 larvae: handling control, n=5; solvent control, n=6; JH, n=8; methoprene, n=6; L3 larvae: handling control, n=6; solvent control, n=9; JH, n=7; methoprene, n=5).

Each sample was homogenized with ceramic beads in a shaker, total RNA was extracted using a modified protocol of the ReliaPrep™ RNA Cell Miniprep System (Promega, USA), and cDNA was synthesized with the iScript™ gDNA Clear cDNA Synthesis Kit (Bio-Rad, USA), following the manufacturer's instructions. cDNA concentration was measured using the Qubit dsDNA HS assay kit (Thermo Fisher Scientific, USA) and standardized to 1.5–3 ng µl−1, depending on developmental stage. Each qPCR reaction was run in 10 µl (5 µl SYBR green, 3 µl purified H2O, 0.5 µl forward primer, 0.5 ml reverse primer, 1 µl cDNA) in a qPCR machine (CFX, Bio-Rad, USA). Three technical replicates were run for each sample and gene. In addition to kr-h1 and the housekeeping gene y45, in all samples we analysed expression of the two sex-specific doublesex isoforms dsxf and dsxm, which vary in their expression between males and females (Klein et al., 2016). This was done to identify males in samples collected from stock colonies containing mated queens, as males cannot be distinguished morphologically from workers before the pupal stage. All male samples thus identified were excluded before statistical analyses. For amplification of dsxF, dsxM and y45, primer sequences designed for a previous study were used (Klein et al., 2016). The housekeeping gene y45 was used for normalization according to the deltaCq method (Schmittgen and Livak, 2008). For timeline data, differences in kr-h1 expression between treatments were compared separately for each sex and timepoint with linear regressions, using log-transformed expression as a response and treatment as explanatory variable: kr-h1 expressionsex,timepoint∼treatment. To assess caste-specific responses to treatment, kr-h1 expression was compared separately for each developmental stage, using log-transformed expression as a response and caste and treatment as explanatory variables: kr-h1 expressiondevelopmental stage∼caste+treatment. For all analyses, model fit was assessed using residual tests implemented in the DHarma package in R and pairwise P-values were calculated and Tukey-corrected for multiple comparisons using the R emmeans package.

Hormone treatment does not affect caste or morph ratios in C. obscurior

A previous study reported that treatment of C. obscurior larvae of unknown sex and caste with the synthetic JH analogue methoprene always resulted in the production of winged morphs – queens and winged males (Schrempf and Heinze, 2006). We first attempted to replicate this result by treating larvae collected from stock colonies without regard to caste-specific crystalline deposits. Compared with a solvent control, methoprene treatment (2 µl of a 1 mg ml−1 solution in 70% ethanol) did not have a significant influence on winged morph production in any of the three larval stages (Table S1, Fisher's exact test, ethanol solvent control versus methoprene: L1, P=0.241; L2, P=0.573; L3, P=0.86).

Similarly, experimental hormone treatment had no effect on the caste of emerging females when queen and worker-destined late-stage embryos and larvae were treated separately (Table 1). Across all developmental stages and treatments, presumed queen-destined individuals developed into adult queens in 87% (346/397) of cases and presumed worker-destined individuals developed into workers in 97% (321/331) of cases. In the three larval stages, the accuracy of caste fate prediction based on crystalline deposit patterns was >90% in individuals subjected to a handling control, and this remained similar following treatment with a solvent control as well as after treatment with JH III or methoprene (Table 1). We confirmed the irreversibility of caste determination by treating queen- and worker destined second instar larvae with higher doses of hormone; in all cases, treated individuals developed into the predicted caste (Table S1). Compared with larvae, queen caste was more difficult to predict in eggs, with only around half of untreated, presumed queen-destined eggs developing into queens (Schultner et al., 2023); this proportion did not change after treatment (Table 1). In contrast, all untreated, presumed worker-destined eggs developed into workers, indicating that worker-destined eggs are more easily identified than queen-destined eggs and/or more likely to end up in experiments because they outnumber queen-destined eggs in stock colonies. The accuracy of caste prediction in early developmental stages furthermore appears to be influenced by sampling precision, as a previous study showed higher prediction accuracy rates for queen-destined eggs, but lower rates for queen-destined first instar larvae (Schultner et al., 2023). Caste prediction accuracy was lowest for methoprene-treated worker-destined eggs, but still relatively high at 80% (12/15), suggesting minor, if any, effects of methoprene on caste once it can be identified by caste-specific crystalline deposits. Of the presumed queen- and worker-destined individuals which emerged as adults, 4.1% (31/759) developed into males, the majority in replicates with presumed worker-destined individuals. Male morph fate was not affected by treatment, as only wingless males emerged (Table 1).

Caste ratios were also not affected by treatment of female early-stage embryos (1–24 h old), with 2.9% (1/35) of solvent control-treated eggs developing into queens, compared with 7% (4/57) of methoprene-treated eggs (Fisher's exact test, odds ratio=0.393, P=0.646). Methoprene-treated eggs showed slightly higher survival [methoprene, 20% (57/285); solvent control, 12.3% (35/285); Fisher's exact test: odds ratio=1.784, P=0.016]. When early-stage male embryos produced by unmated queens were treated, only wingless males emerged in both treatments [solvent control: 100% (10/10), methoprene: 100% (11/11), Fisher's exact test: odds ratio=infinite, P=1], and treatment did not affect male survival, which was extremely low [solvent control: 3.1% (10/325), methoprene: 3.4% (11/325), Fisher's exact test: odds ratio=1.103, P=1].

Juvenile hormone analogues such as methoprene can be used to kill insects because high doses of these chemicals perturb development (Troisi and Riddiford, 1974; Tay and Lee, 2014; Nur Aliah et al., 2021). Similarly, in C. obscurior, survival of treated brood dropped considerably when high doses (5 mg ml−1) were administered (Schrempf and Heinze, 2006). The moderate dose used here did not have adverse effects on survival in queen- and worker-destined individuals, which was similar between treatments in each of the developmental stages (Table 1). Across all developmental stages and treatments, a higher proportion of queen-destined compared with worker-destined individuals survived until pupation (queen-destined: 44.5% (346/893), worker-destined: 38.1% (321/868); Fisher's exact test, odds ratio=0.770, P=0.008); this was mostly driven by differential survival of castes in the third larval instar (Table 1). Adult workers do not discriminate between developing queens and workers, so differential treatment is unlikely to explain this difference (Schultner et al., 2023). Instead, increased survival of queen-destined larvae may stem from size differences between the two castes (Oettler et al., 2019; Schultner et al., 2023).

Male but not female body size is affected by hormone treatment

Ant body size is strongly associated with caste polyphenism, with queens typically larger than workers [see Trible and Kronauer (2017, 2021a,b) and Abouheif (2021) for a discussion about the role of body size in ant caste development]. Diverse factors have been shown to be associated with body size variation in queens and workers, including genotype, maternal effects, nutrition, as well as social and abiotic environment (e.g. Heinze et al., 2003; Hughes et al., 2003; Bargum et al., 2004; Fjerdingstad, 2005; Schwander et al., 2005; Meunier and Chapuisat, 2009; Linksvayer et al., 2011). How caste-specific queen and worker body sizes are attained is not well understood, but studies focusing on size variation within the worker caste have suggested that methylation of epidermal growth factor receptor, which links to JH via the insulin signalling pathway, plays a key role in determining worker size, both in species with distinct worker castes (Alvarado et al., 2015) and in species with monomorphic workers (Renard et al., 2022). JH itself has been implicated in increased body size of Pheidole soldiers (Wheeler and Nijhout, 1981, 1983; Rajakumar et al., 2012) and Pogonomyrmex and Camponotus workers (Helms-Cahan et al., 2011; LeBoeuf et al., 2016).

In C. obscurior, hormone manipulations did not result in larger queens or workers. After treatment in the second larval stage, queen-destined third instar larvae were always larger than worker-destined third instar larvae (Fig. 1, Table 2, Table S2), confirming previous results (Schultner et al., 2023). Treatment had no effect on the size of worker larvae. Methoprene-treated queen-destined larvae were smaller than those subjected to a handling control, but there was no difference between any of the other treatments (Fig. 1, Table 2, Table S2).

Fig. 1.

Body size of queen- and worker-destined larvae in the ant Cardiocondyla obscurior after hormone treatment. (A) Principal component analysis of larval head width, body width and body length separates larvae by overall size on PC1. (B) Queen larvae are larger than worker larvae across treatments. Different letters indicate significant differences between groups. See Table S2 for all Tukey-corrected pairwise P-values. JH III, juvenile hormone III.

Fig. 1.

Body size of queen- and worker-destined larvae in the ant Cardiocondyla obscurior after hormone treatment. (A) Principal component analysis of larval head width, body width and body length separates larvae by overall size on PC1. (B) Queen larvae are larger than worker larvae across treatments. Different letters indicate significant differences between groups. See Table S2 for all Tukey-corrected pairwise P-values. JH III, juvenile hormone III.

Adult queens emerging from treatments of late-stage embryos and larvae were also larger than adult workers (Fig. 2, Table 2, Table S3). This was independent of the timing of hormone treatment, i.e. the developmental stage at which individuals were treated, so that data were pooled across developmental stages for subsequent analyses (Table 2). As in larvae, treatment did not result in larger queens or workers (Fig. 2, Table 2, Table S3). Similarly, treatment had no clear effect on development time until pupation in any of the four treated stages (Table 2). Caste affected development time of late-stage embryos and second instar larvae, with queens exhibiting longer development times than workers; there was no effect of caste on development time of first and third instar larvae (Table 2).

Fig. 2.

Body size of adult queens, workers and males emerging from hormone treatment of late-stage embryos and larvae. (A) Principal component analysis separates individuals by size on PC1. (B) PC 1 scores plotted by caste and treatment show that queens are the largest caste, and workers are larger than males. Female size is not affected by hormone treatment, whereas larger males emerge from methoprene treatments. Different letters indicate significant differences between groups. See Table S3 for all Tukey-corrected pairwise P-values.

Fig. 2.

Body size of adult queens, workers and males emerging from hormone treatment of late-stage embryos and larvae. (A) Principal component analysis separates individuals by size on PC1. (B) PC 1 scores plotted by caste and treatment show that queens are the largest caste, and workers are larger than males. Female size is not affected by hormone treatment, whereas larger males emerge from methoprene treatments. Different letters indicate significant differences between groups. See Table S3 for all Tukey-corrected pairwise P-values.

In contrast to females, males treated with methoprene emerged as significantly larger adults than males emerging from the other treatments (Table 2; Fig S1, Table S3). Wingless males are typically smaller than females (Oettler et al., 2019), but males emerging from methoprene treatments exhibited body sizes similar to those of workers (Fig. 2, Table 2; Table S3).

Adult queens emerging from treatment of female early-stage embryos of unknown caste were larger than workers and males, and again, methoprene treatment did not result in larger body size of females (Table 2; Fig S1, Table S4). Male body size did respond to treatment of early-stage embryos, but in an unexpected manner: males emerging from methoprene treatments tended to exhibit decreased thorax lengths and petiole widths compared with solvent-control treated males (Table 2; Fig. S1, Table S4).

Conserved expression of a JH first-response gene in queens, workers and males

Molecular pathways associated with juvenile hormone have been extensively studied in a few model organisms such as the tobacco hornworm Manduca sexta and the flour beetle Tribolium castaneum (Riddiford, 2012). These have revealed a suite of conserved genes acting downstream of JH, including methoprene-tolerant and krüppel-homolog 1 (Kr-h1) (Konopova and Jindra, 2007; Minakuchi et al., 2008, 2009; Jindra et al., 2013). Kr-h1 expression in social Hymenoptera has mainly been studied in adults, and varies with task in bee workers (Whitfield et al., 2003; Grozinger and Robinson, 2007; Shpigler et al., 2010; Johnson and Jasper, 2016), and with caste (Glastad et al., 2021; Zhang et al., 2021) and reproductive status in ants (Araki et al., 2020; Gospocic et al., 2021). During worker development, pharmacologically induced hypomethylation resulting in increased expression of the growth-promoting egfr gene was associated with decreased kr-h1 expression, and an increase in body size (Renard et al., 2022). Similarly, knockout of an oestrogen-related receptor resulted in decreased kr-h1 expression (Zhang et al., 2021). While these studies have helped confirm the molecular links between insulin and hormone signalling, growth and reproduction, it remains unknown whether the function of kr-h1 as a JH first-response gene is conserved across sexes and castes in developing ants. We used quantitative PCR to measure kr-h1 expression following hormone treatment in C. obscurior. In female larvae, kr-h1 expression increased significantly 6 and 24 h after hormone treatment compared with a solvent control (Fig. 3, linear regression, factor: treatment, 1 h: F2,39=0.661, P=0.522; 2 h: F2,29=2.044, P=0.144; 6 h: F2,38=3.997, P=0.027; 24 h: F2,40=10.016, P<0.001). In males, hormone treatment led to increased kr-h1 expression after 1, 2 and 24 h (Fig. 3, linear regression, factor: treatment, 1 h: F1,24=6.767, P=0.016; 2 h: F1,23=15.993, P<0.001; 6 h: F2,12=1.27, P=0.282; 24 h: F1,23=4.665, P=0.041).

Fig. 3.

Timeline of krüppel-homolog 1 expression in female and male larvae after hormone treatment. Expression in larvae from the three larval instars was compared between treatments, separately for each sex and timepoint, using linear regressions followed by Tukey correction of P-values (*P<0.05, **P<0.01).

Fig. 3.

Timeline of krüppel-homolog 1 expression in female and male larvae after hormone treatment. Expression in larvae from the three larval instars was compared between treatments, separately for each sex and timepoint, using linear regressions followed by Tukey correction of P-values (*P<0.05, **P<0.01).

Krüppel-homolog 1 expression was not affected by the caste of late-stage embryos or larvae (linear regression, factor: caste, eggs: F1,77=0.045, P=0.833; L1 larvae: F1,54=3.785, P=0.057; L3 larvae: F1,42=0.035, P=0.853; L2 larvae: as model residuals did not fit assumptions of normality, linear regressions comparing caste-specific expression were run separately by treatment: handling control: F1,8=0.033, P=0.86; solvent control: F1,12=0.188, P=0.672; JH III: F1,13=4.139, P=0.063; methoprene: F1,10=0.032, P=0.862). There was furthermore no effect of treatment in late-stage embryos, first instar and second instar larvae (linear regression, factor: treatment, eggs: F3,77=0.522, P=0.668; L1 larvae: F3,54=1.934, P=0.135; L2 larvae: as model residuals did not fit assumptions of normality, linear regressions comparing expression by treatment were run separately in each caste: queens: F3,23=1.392, P=0.86; workers: F1,20=3.016, P=0.054). Third instar larvae of both castes exhibited significantly higher kr-h1 expression after JH III treatment compared with handling control and solvent control-treated larvae (linear regression, factor: treatment, L3 larvae, F3,42=5.574, P=0.002; pairwise-corrected P-values: worker-destined L3 larvae: handling control versus JH III, P=0.037; solvent control versus JH III, P=0.027; queen-destined L3 larvae: handling control versus JH III, P=0.037; solvent control versus JH III, P=0.027).

Juvenile hormone studies in ants have employed a number of different methods, from topical application with fine capillaries or brushes to supplementation of liquid and solid food, to physical contact with hormone-soaked objects. In addition, a wide range of doses have been used and sometimes not fully reported, making it difficult to compare between studies. Moreover, mortality was not always reported accurately, or controls were missing. This is further complicated by the range of factors associated with caste determination in ants, including genotype, maternal effects, nutrition, social environment, temperature and combinations thereof (e.g. Gösswald and Bier, 1953; Brian, 1975; Passera, 1980; Helms-Cahan et al., 2002; Anderson et al., 2008; Smith et al., 2008; Penick and Liebig, 2012), and the consequences for the timing of determination. Accordingly, the interpretation of experimental studies on the role of JH in ant queen–worker caste development is not straightforward.

In Cataglyphis mauritanica and Pogonomyrmex barbatus× rugosus, two species with a genetic component to caste determination (Helms-Cahan et al., 2002; Julian et al., 2002), methoprene treatment of eggs (C. mauritanica, Kuhn et al., 2018) but not queens (P. barbatus× rugosus, Helms-Cahan et al., 2011) increased production of queens with worker genotypes. Colonies with treated P. barbatus× rugosus queens produced larger workers, but this presumed hormone effect cannot be disentangled from treatment-induced changes to colony size (Helms-Cahan et al., 2011). In P. pallidula, only eggs produced by overwintered queens shortly after hibernation can give rise to new queens (Passera, 1980); treating queens and eggs with JH I outside this time period resulted in queen production, whereas treatment of larvae had no effect (Passera and Suzzoni, 1979). Analogous to results on JH, differences in ecdysteroid titers have also been detected in queen- and worker-destined eggs and larvae (Suzzoni et al., 1980, 1983). In another Pheidole species, methoprene treatment of queens, eggs and larvae did not result in queen production (Ono, 1982). In Pogonomyrmex rugosus, as in P. pallidula, only hibernated queens are capable of producing queen-potent eggs (Schwander et al., 2008); here, feeding hibernated colonies with methoprene strongly increased the proportion of queens produced (Libbrecht et al., 2013), although colony productivity and brood mortality were not monitored. In Solenopsis invicta, JH-analogue treatment of larvae but not eggs or queens resulted in queen production (Vinson and Robeau, 1974; Robeau and Vinson, 1976; Banks et al., 1978). Finally, in three species with caste determination during larval development, topical treatment with JH analogues resulted in development of more and larger queens in Myrmica rubra (Brian, 1974) and increased queen production in Harpegnathos saltator (Penick et al., 2012), but feeding with JH III had no effect on queen production in Camponotus floridanus, although worker size increased (LeBoeuf et al., 2016).

How does C. obscurior fit into this admittedly confusing picture? A previous study found that methoprene treatment of larvae resulted in increased production of winged individuals (Schrempf and Heinze, 2006) but we could not replicate this result. Although we did not systematically test different doses of methoprene, it is unlikely that this discrepancy is caused by a dose effect as the doses used here (∼400–40,000 µg g−1 body weight depending on the experiment) were within the general range reported in other studies (∼60–300,000 µg g−1 body weight) and much higher than natural titers of JH reported for insects. For example, whole body extracts of adult female mosquitos contained ∼0.008 µg g−1 JH (Zhao et al., 2016) and JH titers in queens of the fire ant Solenopsis invicta ranged from 0.008 pmol to a maximum of ∼0.02 pmol (Brent and Vargo, 2003). The efficacy of doses used in the present study is further supported by the result that JH treatment led to increased expression of the conserved gene krüppel-homolog 1. There are several plausible explanations for the conflicting results obtained by Schrempf and Heinze (2006). First, in the present study queen-destined larvae showed a trend toward increased survival compared with workers, so that previous results may have been confounded by differential survival of the two castes. Second, castes differ in size by the second larval instar (Schultner et al., 2023), which may have led to a sampling bias in the previous experiment, so that more larger, queen destined, larvae were included. Third, caste proportions obtained from solvent control treatments (acetone) were not reported in the previous study. In our experiments, treatment of second instar larvae with acetone resulted in a worker-biased caste ratio (Table S1), indicating side-effects of this solvent. In contrast to Schrempf and Heinze (2006), we therefore used ethanol as a solvent for the majority of experiments, also because we found no general negative effects of ethanol on survival (Table 1). Finally, the studies used two different populations of C. obscurior, which are known to differ in key traits such as body size, behavior, genome and microbiome (Schrader et al., 2014; Errbii et al., 2021; Ün et al., 2021a,b), as well as the tendency to produce winged males (E.S., unpublished data).

Together with the direct lines of evidence for extreme morphological and molecular differentiation between embryos of the two castes (Wallner, 2021; Oettler et al., 2023 preprint; Schultner et al., 2023), the results of the present study indicate that queen and worker growth and development are highly canalized in C. obscurior once caste has been determined. It remains to be tested whether JH treatment of queens leads to the increased production of queen-destined eggs, as in other species with maternal caste determination. In contrast to females, male development appears to have retained a substantial level of plasticity, perhaps because male morph is determined during larval development. The different degrees of canalization in males and females may also be explained by the evolutionary age of male polymorphism, which although basal to Cardiocondyla (Oettler et al., 2010), evolved later than queen–worker polyphenism.

Juvenile hormone has long been considered the holy grail in ant queen–worker caste polyphenism, even though no studies have followed up on the original reports in S. invicta, P. pallidula and M. rubra. Wheeler herself acknowledged that after this first intense phase of study, investigations into JH have not been as fruitful as initially hoped (Wheeler, 2003). From the studies summarized above as well as our own results, we nevertheless believe that some conclusions can be drawn. First, the mode and the timing of caste determination are clearly critical factors which can influence JH-responsiveness. Second, JH treatment may override even seemingly ‘hard-wired’ caste determination modes such as the hybridogenetic system in Cataglyphis mauritanica, but the strength of this effect appears to depend on additional individual- and/or colony-level traits. More studies on species with genetic modes of caste determination are needed to validate these findings. Third, any influence experimental JH manipulation may have on caste and body size can be mediated by mortality rates, as these affect brood:worker ratios and colony size, two important social factors in ant development. To avoid this caveat, experiments should ideally be done with standardized colony sizes. Finally, replication of previous experiments coupled with modern molecular methods is needed to evaluate general assumptions about the role of JH in ant caste development, and to help to move this traditional research field into a new era.

The authors thank H. Lowack, B. Grüneberg, S. Stewart, J. Wallner, L. Drechsel, H. Sammereier and G. Gebhard for help with experiments.

Author contributions

Conceptualization: J.B., J.O., J.H., E.S.; Methodology: J.B., E.S.; Validation: E.S.; Formal analysis: J.B., A.S., E.S.; Investigation: J.B., A.S., E.S.; Resources: J.H., E.S.; Data curation: J.B., A.S., E.S.; Writing - original draft: J.B., E.S.; Writing - review & editing: J.B., J.O., J.H., E.S.; Visualization: E.S.; Supervision: E.S.; Project administration: E.S.; Funding acquisition: J.B., J.O., J.H., E.S.

Funding

This work was funded by a research stipend from the Deutscher Akademischer Austauschdienst to J.B., a research grant from the Deutsche Forschungsgemeinschaft (DFG) to J.O. (Oe549/3-1,2) and a research stipend from the Bayerisches Staatsministerium für Wissenschaft und Kunst to E.S. Funders had no role in the design of the study, data collection or interpretation of results. Open Access funding provided by University of Regensburg. Deposited in PMC for immediate release.

Data availability

All relevant data can be found within the article and its supplementary information.

Abouheif
,
E.
(
2021
).
Ant caste evo-devo: it's not all about size
.
Trends Ecol. Evol.
36
,
668
-
670
.
Alvarado
,
S.
,
Rajakumar
,
R.
,
Abouheif
,
E.
and
Szyf
,
M.
(
2015
).
Epigenetic variation in the Egfr gene generates quantitative variation in a complex trait in ants
.
Nat. Commun.
6
,
1
-
9
.
Anderson
,
K.
,
Linksvayer
,
T. A.
and
Smith
,
C. R.
(
2008
).
The causes and consequences of genetic caste determination in ants (Hymenoptera: Formicidae)
.
Myrmecological News
11
,
119
-
132
.
Araki
,
M.
,
Miyakawa
,
M. O.
,
Suzuki
,
T.
and
Miyakawa
,
H.
(
2020
).
Two insulin–like peptides may regulate egg production in opposite directions via juvenile hormone signaling in the queenless ant Pristomyrmex punctatus
.
J. Exp. Zool. B Mol. Dev. Evol.
334
,
225
-
234
.
Banks
,
W. A.
,
Lofgren
,
C. S.
and
Plumley
,
J. K.
(
1978
).
Red imported fire ants: effects of insect growth regulators on caste formation and colony growth and survival
.
J. Econ. Entomol.
71
,
75
-
78
.
Bargum
,
K.
,
Boomsma
,
J.
and
Sundström
,
L.
(
2004
).
A genetic component to size in queens of the ant, Formica truncorum
.
Behav. Ecol. Sociobiol.
57
,
9
-
16
.
Bier
,
K.
(
1952
).
Beziehungen zwischen Nährzellkerngröße und Ausbildung ribonukleinsäurehaltiger Strukturen in den Oocyten von Formica rufa rufo-pratensis minor Gößwald
.
Verh. Dtsch. Zool. Ges. Freiburg
40
,
369
-
374
.
Brent
,
C. S.
and
Vargo
,
E. L.
(
2003
).
Changes in juvenile hormone biosynthetic rate and whole body content in maturing virgin queens of Solenopsis invicta
.
J. Insect Physiol.
49
,
967
-
974
.
Brian
,
M. V.
(
1954
).
Studies of caste differentiation in Myrmica rubra L. 1. The growth of queens and males
.
Insectes Soc.
1
,
101
-
122
.
Brian
,
M. V.
(
1955
).
Studies of caste differentiation in Myrmica rubra L. 3. Larval dormancy, winter size and vernalisation
.
Insectes Soc.
2
,
85
-
114
.
Brian
,
M. V.
(
1973
).
Caste control through worker attack in the ant Myrmica
.
Insectes Soc.
20
,
87
-
102
.
Brian
,
M. V.
(
1974
).
Caste differentiation in Myrmica rubra: the role of hormones
.
J. Insect Physiol.
20
,
1351
-
1365
.
Brian
,
M. V.
(
1975
).
Caste determinaton through a queen influence on diapause in larvae of the ant Myrmica rubra
.
Entomol. Exp. Appl.
18
,
429
-
442
.
de Menten
,
L.
,
Fournier
,
D.
,
Brent
,
C. S.
,
Passera
,
L.
,
Vargo
,
E. L.
and
Aron
,
S.
(
2005
).
Dual mechanism of queen influence over sex ratio in the ant Pheidole pallidula
.
Behav. Ecol. Sociobiol.
58
,
527
-
533
.
Errbii
,
M.
,
Keilwagen
,
J.
,
Hoff
,
K. J.
,
Steffen
,
R.
,
Altmüller
,
J.
,
Oettler
,
J.
and
Schrader
,
L.
(
2021
).
Transposable elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior
.
Mol. Ecol.
30
,
6211
-
6228
.
Fjerdingstad
,
E. J.
(
2005
).
Control of body size of Lasius niger ant sexuals - worker interests, genes and environment
.
Mol. Ecol.
14
,
3123
-
3132
.
Glastad
,
K. M.
,
Ju
,
L.
and
Berger
,
S. L.
(
2021
).
Tramtrack acts during late pupal development to direct ant caste identity
.
PLoS Genet.
17
,
1
-
22
.
Gospocic
,
J.
,
Glastad
,
K. M.
,
Sheng
,
L.
,
Shields
,
E. J.
,
Berger
,
S. L.
and
Bonasio
,
R.
(
2021
).
Kr-h1 maintains distinct caste-specific neurotranscriptomes in response to socially regulated hormones
.
Cell
184
,
5807
-
5823.e14
.
Gösswald
,
K.
and
Bier
,
K.
(
1953
).
Untersuchungen zur Kastendetermination in der Gattung Formica
.
Naturwissenschaften
40
,
38
-
39
.
Grozinger
,
C. M.
and
Robinson
,
G. E.
(
2007
).
Endocrine modulation of a pheromone-responsive gene in the honey bee brain
.
J. Comp. Physiol. A Neuroethol. Sens. Neural Behav. Physiol.
193
,
461
-
470
.
Heinze
,
J.
(
2017
).
Life-history evolution in ants: the case of Cardiocondyla
.
Proc. R. Soc. B
284
,
20161406
.
Heinze
,
J.
,
Foitzik
,
S.
,
Fischer
,
B.
,
Wanke
,
T.
and
Kipyatkov
,
V. E.
(
2003
).
The significance of latitudinal variation in body size in a holarctic ant, Leptothorax acervorum
.
Ecography
26
,
349
-
355
.
Heinze
,
J.
and
Hölldobler
,
B.
(
1993
).
Fighting for a harem of queens: physiology of reproduction in Cardiocondyla male ants
.
Proc. Natl. Acad. Sci. USA
90
,
8412
-
8114
.
Helms-Cahan
,
S.
,
Parker
,
J. D.
,
Rissing
,
S. W.
,
Johnson
,
R. A.
,
Polony
,
T. S.
,
Weiser
,
M. D.
and
Smith
,
D. R.
(
2002
).
Extreme genetic differences between queens and workers in hybridizing Pogonomyrmex harvester ants
.
Proc. R. Soc. B
269
,
1871
-
1877
.
Helms-Cahan
,
S.
,
Graves
,
C. J.
and
Brent
,
C. S.
(
2011
).
Intergenerational effect of juvenile hormone on offspring in Pogonomyrmex harvester ants
.
J. Comp. Physiol. B Biochem. Syst. Envir. Physiol.
181
,
991
-
999
.
Hughes
,
W. O. H.
,
Sumner
,
S.
,
Van Borm
,
S.
, and
Boomsma
,
J.
(
2003
).
Worker caste polymorphism has a genetic basis in Acromyrmex leaf-cutting ants
.
Proc. Natl. Acad. Sci. USA
100
,
9394
-
9397
.
Jaimes-Nino
,
L. M.
,
Heinze
,
J.
and
Oettler
,
J.
(
2022
).
Late-life fitness gains and reproductive death in Cardiocondyla obscurior ants
.
eLife
11
,
e74695
.
Jindra
,
M.
,
Palli
,
S. R.
and
Riddiford
,
L. M.
(
2013
).
The juvenile hormone signaling pathway in insect development
.
Annu. Rev. Entomol.
58
,
181
-
204
.
Johnson
,
B. R.
and
Jasper
,
W. C.
(
2016
).
Complex patterns of differential expression in candidate master regulatory genes for social behavior in honey bees
.
Behav. Ecol. Sociobiol.
70
,
1033
-
1043
.
Julian
,
G. E.
,
Fewell
,
J. H.
,
Gadau
,
J.
,
Johnson
,
R. A.
and
Larrabee
,
D.
(
2002
).
Genetic determination of the queen caste in an ant hybrid zone
.
Proc. Natl. Acad. Sci. USA
99
,
8157
-
8160
.
Kinomura
,
K.
and
Yamauchi
,
K.
(
1987
).
Fighting and mating behaviors of dimorphic males in the ant Cardiocondyla wroughtoni
.
J. Ethol.
5
,
75
-
81
.
Klein
,
A.
,
Schultner
,
E.
,
Lowak
,
H.
,
Schrader
,
L.
,
Heinze
,
J.
,
Holman
,
L.
and
Oettler
,
J.
(
2016
).
Evolution of social insect polyphenism facilitated by the sex differentiation cascade
.
PLoS Genet.
12
,
e1005952
.
Konopova
,
B.
and
Jindra
,
M.
(
2007
).
Juvenile hormone resistance gene Methoprene-tolerant controls entry into metamorphosis in the beetle Tribolium castaneum
.
Proc. Natl. Acad. Sci. USA
104
,
10488
-
10493
.
Kuhn
,
A.
,
Darras
,
H.
and
Aron
,
S.
(
2018
).
Phenotypic plasticity in an ant with strong caste-genotype association
.
Biol. Lett.
14
,
20170705
.
LeBoeuf
,
A.
,
Waridel
,
P.
,
Brent
,
C. S.
,
Goncalves
,
A.
,
Menin
,
L.
,
Ortiz
,
D.
,
Riba-Grognuz
,
O.
,
Koto
,
A.
,
Soares
,
Z. G.
,
Privman
,
E.
et al. 
(
2016
).
Oral transfer of chemical cues, growth proteins and hormones in social insects
.
eLife
5
,
e20375
.
Libbrecht
,
R.
,
Corona
,
M.
,
Wende
,
F.
,
Azevedo
,
D. O.
,
Serrao
,
J. E.
and
Keller
,
L.
(
2013
).
Interplay between insulin signaling, juvenile hormone, and vitellogenin regulates maternal effects on polyphenism in ants
.
Proc. Natl. Acad. Sci. USA
110
,
11050
-
11055
.
Linksvayer
,
T. A.
,
Kaftanoglu
,
O.
,
Akyol
,
E.
,
Blatch
,
S.
,
Amdam
,
G. V.
and
Page
,
R. E.
Jr
(
2011
).
Larval and nurse worker control of developmental plasticity and the evolution of honey bee queen-worker dimorphism
.
J. Evol. Biol.
24
,
1939
-
1948
.
Meunier
,
J. M.
and
Chapuisat
,
M.
(
2009
).
The determinants of queen size in a socially polymorphic ant
.
J. Evol. Biol.
22
,
1906
-
1913
.
Minakuchi
,
C.
,
Zhou
,
X.
and
Riddiford
,
L. M.
(
2008
).
Krüppel homolog 1 (Kr-h1) mediates juvenile hormone action during metamorphosis of Drosophila melanogaster
.
Mech. Dev.
125
,
91
-
105
.
Minakuchi
,
C.
,
Namiki
,
T.
and
Shinoda
,
T.
(
2009
).
Krüppel homolog 1, an early juvenile hormone-response gene downstream of Methoprene-tolerant, mediates its anti-metamorphic action in the red flour beetle Tribolium castaneum
.
Dev. Biol.
325
,
341
-
350
.
Nijhout
,
H. F.
(
2003
).
The control of body size in insects
.
Dev. Biol.
261
,
1
-
9
.
Nijhout
,
H.
and
Wheeler
,
D. E.
(
1982
).
Juvenile hormone and the physiological basis of insect polymorphisms
.
Q. Rev. Biol.
57
,
109
-
133
.
Nur Aliah
,
N. A.
,
Ab-Rahim
,
S.
,
Moore
,
H. E.
and
Heo
,
C. C.
(
2021
).
Juvenile hormone: Production, regulation, current application in vector control and its future applications
.
Trop. Biomed.
38
,
254
-
264
.
Oettler
,
J.
(
2020
).
Cardiocondyla: Heart node ants
. In
Encyclopedia of Social Insects
(
ed.
C.
Starr
), pp.
1
-
3
.
Cham
:
Springer Nature Switzerland
.
Oettler
,
J.
,
Suefuji
,
M.
and
Heinze
,
J.
(
2010
).
The evolution of alternative reproductive tactics in male Cardiocondyla ants
.
Evolution
64
,
3310
-
3317
.
Oettler
,
J.
,
Platschek
,
T.
,
Schmidt
,
C.
,
Rajakumar
,
R.
,
Favé
,
M.-J.
,
Khila
,
A.
,
Heinze
,
J.
and
Abouheif
,
E.
(
2019
).
Interruption points in the wing gene regulatory network underlying wing polyphenism evolved independently in male and female morphs in Cardiocondyla ants
.
J. Exp. Zool. B Mol. Dev. Evol.
332
,
7
-
16
.
Oettler
,
J.
,
Wallner
,
T.
,
Dofka
,
B.
,
Heinze
,
J.
,
Eichner
,
N.
,
Meister
,
G.
,
Errbii
,
M.
,
Rehli
,
M.
,
Gebhard
,
C.
and
Schultner
,
E.
(
2023
).
Darwin's “neuters” and the evolution of the sex continuum in a superorganism
.
bioRxiv
.
Ono
,
S.
(
1982
).
Effect of juvenile hormone on the caste determination in the ant, Pheidole fervida
.
Appl. Entomol. Zool.
17
,
1
-
7
.
Passera
,
L.
(
1980
).
La ponte d’œufs préorientés chez la fourmi Pheidole pallidula (Nyl.) (Hymenoptera-Formicidae)
.
Insectes Soc.
27
,
79
-
95
.
Passera
,
L.
and
Suzzoni
,
J. P.
(
1979
).
Le role de la reine de Pheidole pallidula (Nyl.) (Hymenoptera, Formicidae) dans la sexualisation du couvain après traitement par l'hormone juvénile
.
Insectes Soc.
26
,
343
-
353
.
Pearcy
,
M.
,
Aron
,
S.
,
Doums
,
C.
and
Keller
,
L.
(
2004
).
Conditional use of sex and parthenogenesis for worker and queen production in ants
.
Science
306
,
1780
-
1783
.
Penick
,
C. A.
and
Liebig
,
J.
(
2012
).
Regulation of queen development through worker aggression in a predatory ant
.
Behav. Ecol.
23
,
992
-
998
.
Penick
,
C. A.
,
Prager
,
S. S.
and
Liebig
,
J.
(
2012
).
Juvenile hormone induces queen development in late-stage larvae of the ant Harpegnathos saltator
.
J. Insect Physiol.
58
,
1643
-
1649
.
Rajakumar
,
R.
,
San Mauro
,
D.
,
Dijkstra
,
M. B.
,
Huang
,
M. H.
,
Wheeler
,
D. E.
,
Hiou-Tim
,
F.
,
Khila
,
A.
,
Cournoyea
,
M.
and
Abouheif
,
E.
(
2012
).
Ancestral developmental potential facilitates parallel evolution in ants
.
Science
335
,
79
-
82
.
Renard
,
T.
,
Gueydan
,
C.
and
Aron
,
S.
(
2022
).
DNA methylation and expression of the egfr gene are associated with worker size in monomorphic ants
.
Sci. Rep.
12
,
1
-
9
.
Riddiford
,
L. M.
(
2012
).
How does juvenile hormone control insect metamorphosis and reproduction?
Gen. Comp. Endocrinol.
179
,
477
-
484
.
Robeau
,
R.
and
Vinson
,
S.
(
1976
).
Effects of juvenile hormone analogues on caste differentiation in the imported fire ant, Solenopsis invicta
.
J. Georgia Entomol. Soc.
11
,
198
-
203
.
Schmittgen
,
T. D.
and
Livak
,
K. J.
(
2008
).
Analyzing real-time PCR data by the comparative CT method
.
Nat. Protoc.
3
,
1101
-
1108
.
Schrader
,
L.
,
Kim
,
J. W.
,
Ence
,
D.
,
Zimin
,
A.
,
Klein
,
A.
,
Wyschetzki
,
K.
,
Weichselgartner
,
T.
,
Kemena
,
C.
,
Stökl
,
J.
,
Schultner
,
E.
et al. 
(
2014
).
Transposable element islands facilitate adaptation to novel environments in an invasive species
.
Nat. Commun.
5
,
5495
.
Schrader
,
L.
,
Simola
,
D.
,
Heinze
,
J.
and
Oettler
,
J.
(
2015
).
Sphingolipids, transcription factors and conserved tool kit genes: developmental plasticity in the ant Cardiocondyla obscurior
.
Mol. Biol. Evol.
32
,
1474
-
1486
.
Schrempf
,
A.
and
Heinze
,
J.
(
2006
).
Proximate mechanisms of male morph determination in the ant Cardiocondyla obscurior
.
Evol. Dev.
8
,
266
-
272
.
Schultner
,
E.
and
Pulliainen
,
U.
(
2020
).
Brood recognition and discrimination in ants
.
Insectes Soc.
67
,
11
-
34
.
Schultner
,
E.
,
Wallner
,
T.
,
Dofka
,
B.
,
Brülhart
,
J.
,
Heinze
,
J.
,
Freitak
,
D.
,
Pokorny
,
T.
and
Oettler
,
J.
(
2023
).
Queens control caste allocation in the ant Cardiocondyla obscurior
.
Proc. R. Soc. B
290
,
20221784
.
Schwander
,
T.
,
Rosset
,
H.
and
Chapuisat
,
M.
(
2005
).
Division of labour and worker size polymorphism in ant colonies: the impact of social and genetic factors
.
Behav. Ecol. Sociobiol.
59
,
215
-
221
.
Schwander
,
T.
,
Humbert
,
J.-Y.
,
Brent
,
C. S.
,
Helms-Cahan
,
S.
,
Chapuis
,
L.
,
Renai
,
E.
and
Keller
,
L.
(
2008
).
Maternal effect on female caste determination in a social insect
.
Curr. Biol.
18
,
265
-
269
.
Shpigler
,
H.
,
Patch
,
H. M.
,
Cohen
,
M.
,
Fan
,
Y.
,
Grozinger
,
C. M.
and
Bloch
,
G.
(
2010
).
The transcription factor Krüppel homolog 1 is linked to hormone mediated social organization in bees
.
BMC Evol. Biol.
10
,
120
.
Smith
,
C. R.
,
Anderson
,
K.
,
Tillberg
,
C. V.
,
Gadau
,
J.
and
Suarez
,
A. V.
(
2008
).
Caste determination in a polymorphic social insect: nutritional, social, and genetic factors
.
Am. Nat.
172
,
497
-
507
.
Stacklies
,
W.
,
Redestig
,
H.
,
Scholz
,
M.
,
Walther
,
D.
and
Selbig
,
J.
(
2007
).
pcaMethods -- a Bioconductor package providing PCA methods for incomplete data
.
Bioinformatics
23
,
1164
-
1167
.
Suzzoni
,
J. P.
,
Passera
,
L.
and
Strambi
,
A.
(
1980
).
Ecdysteroid tire and caste determination in the ant, Pheidole pallidula (Nyl.) (Hymenoptera: Formicidae)
.
Experientia
36
,
1228
-
1229
.
Suzzoni
,
J.-P.
,
Passera
,
L.
and
Strambi
,
A.
(
1983
).
Ecdysteroid production during caste differentiation in larvae of the ant, Plagiolepis pygmaea
.
Physiol. Entomol.
8
,
93
-
96
.
Tay
,
J. W.
and
Lee
,
C. Y.
(
2014
).
Influences of pyriproxyfen on fecundity and reproduction of the pharaoh ant (Hymenoptera: Formicidae)
.
J. Econ. Entomol.
107
,
1216
-
1223
.
Trible
,
W.
and
Kronauer
,
D. J. C.
(
2017
).
Caste development and evolution in ants: it's all about size
.
J. Exp. Biol.
220
,
53
-
62
.
Trible
,
W.
and
Kronauer
,
D. J. C.
(
2021a
).
Hourglass model for developmental evolution of ant castes
.
Trends Ecol. Evol.
36
,
100
-
103
.
Trible
,
W.
and
Kronauer
,
D. J. C.
(
2021b
).
Ant caste evo-devo: size predicts caste (almost) perfectly
.
Trends Ecol. Evol.
36
,
671
-
673
.
Troisi
,
S. J.
and
Riddiford
,
L. M.
(
1974
).
Juvenile hormone effects on metamorphosis and reproduction of the fire ant, Solenopsis invicta
.
Environ. Entomol.
3
,
112
-
116
.
Ün
,
Ç.
,
Hacker
,
J.
,
Heinze
,
J.
,
Oettler
,
J.
and
Schultner
,
E.
(
2021a
).
Endosymbionts mediate the effects of antibiotic exposure in the tramp ant Cardiocondyla obscurior
.
Ecol. Entomol.
47
,
170
-
180
.
Ün
,
Ç.
,
Schultner
,
E.
,
Manzano–Marín
,
A.
,
Flórez
,
L. V.
,
Seifert
,
B.
,
Heinze
,
J.
and
Oettler
,
J.
(
2021b
).
Cytoplasmic incompatibility between Old and New World populations of a tramp ant
.
Evolution
75
,
1775
-
1791
.
Vargo
,
E. L.
and
Passera
,
L.
(
1992
).
Gyne development in the Argentine ant Iridomyrmex humilis: role of overwintering and queen control
.
Physiol. Entomol.
17
,
193
-
201
.
Vinson
,
S. B.
and
Robeau
,
R.
(
1974
).
Insect growth regulator effects on colonies of the imported fire ant
.
J. Econ. Entomol.
67
,
584
-
587
.
Wallner
,
T.
(
2021
).
Caste differentiation in Cardiocondyla obscurior
. PhD thesis, University of Regensburg. doi:10.5283/epub.50832
Wesson
,
L. G.
(
1940
).
An experimental study on caste determination in ants
.
Psyche
47
,
105
-
111
.
Wheeler
,
D. E.
(
1986
).
Developmental and physiological determinants of caste in social Hymenoptera: evolutionary implications
.
Am. Nat.
128
,
13
-
34
.
Wheeler
,
D. E.
(
2003
).
One hundred years of caste determination in Hymenoptera
. In
Genes, Behaviors and Evolution of Social Insects
(
ed.
Kikuchi
,
T.
,
Azuma
,
N.
,
Higashi
,
S.
), pp.
35
-
53
.
Sapporo, Japan
:
Hokkaido University Press
.
Wheeler
,
D. E.
and
Nijhout
,
H. F.
(
1981
).
Soldier determination in ants: new role for juvenile hormone
.
Science
213
,
361
-
363
.
Wheeler
,
D. E.
and
Nijhout
,
H. F.
(
1983
).
Soldier determination in Pheidole bicarinata: effect of methoprene on caste and size within castes
.
J. Insect Physiol.
29
,
847
-
854
.
Whitfield
,
C.
,
Cziko
,
A.-M.
and
Robinson
,
G. E.
(
2003
).
Gene expression profiles in the brain predict behavior in individual honeybees
.
Science
302
,
296
-
299
.
Zhang
,
J. J.
,
Xi
,
G. S.
and
Zhao
,
J.
(
2021
).
Vitellogenin regulates estrogen-related receptor expression by crosstalk with the JH and IIS-TOR signaling pathway in Polyrhachis vicina Roger (Hymenoptera, Formicidae)
.
Gen. Comp. Endocrinol.
310
,
113836
.
Zhao
,
B.
,
Hou
,
Y.
,
Wang
,
J.
,
Kokoza
,
V. A.
,
Saha
,
T. T.
,
Wang
,
X. L.
,
Lin
,
L.
,
Zou
,
Z.
and
Raikhel
,
A. S.
(
2016
).
Determination of juvenile hormone titers by means of LC-MS/MS/MS and a juvenile hormone-responsive Gal4/UAS system in Aedes aegypti mosquitoes
.
Insect Biochem. Mol. Biol.
77
,
69
-
77
.

Competing interests

The authors declare no competing or financial interests.

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.

Supplementary information