Oral transfer of chemical cues, growth proteins and hormones in social insects
Abstract
Social insects frequently engage in oral fluid exchange – trophallaxis – between adults, and between adults and larvae. Although trophallaxis is widely considered a food-sharing mechanism, we hypothesized that endogenous components of this fluid might underlie a novel means of chemical communication between colony members. Through protein and small-molecule mass spectrometry and RNA sequencing, we found that trophallactic fluid in the ant Camponotus floridanus contains a set of specific digestion- and non-digestion related proteins, as well as hydrocarbons, microRNAs, and a key developmental regulator, juvenile hormone. When C. floridanus workers’ food was supplemented with this hormone, the larvae they reared via trophallaxis were twice as likely to complete metamorphosis and became larger workers. Comparison of trophallactic fluid proteins across social insect species revealed that many are regulators of growth, development and behavioral maturation. These results suggest that trophallaxis plays previously unsuspected roles in communication and enables communal control of colony phenotypes.
https://doi.org/10.7554/eLife.20375.001eLife digest
Ants, bees and other social insects live in large colonies where all the individuals work together to gather food, rear young and defend the colony. This level of cooperation requires the insects in the colony to be able to communicate with each other.
Most social insects share fluid mouth-to-mouth with other individuals in their colony. This behavior, called trophallaxis, allows these species to pass around food, both between adults, and between adults and larvae. Trophallaxis therefore creates a network of interactions linking every member of the colony. With this in mind, LeBoeuf et al. investigated whether trophallaxis may also be used by ants to share information relevant to the colony as a form of chemical communication.
The experiments show that in addition to food, carpenter ants also pass small ribonucleic acid (RNA) molecules, chemical signals that help them recognize nestmates, and many proteins that appear to be involved in regulating the growth of ants. LeBoeuf et al. also found that trophallactic fluid contains juvenile hormone, an important regulator of insect growth and development. Adding juvenile hormone to the food that adult ants pass to the larvae made it twice as likely that the larvae would survive to reach adulthood. This indicates that proteins and other molecules transferred mouth-to-mouth over this social network could be used by the ants to regulate how the colony develops.
The next steps following on from this work will be to investigate the roles of the other components of trophallactic fluid, and to examine how individual ants adapt the contents of the fluid in different social and environmental conditions. Another challenge will be to determine how specific components passed to larvae in this way can control their growth and development.
https://doi.org/10.7554/eLife.20375.002Introduction
Many fluids shared between individuals of the same species, such as milk or semen, can exert significant physiological effects on recipients (Poiani, 2006; Liu and Kubli, 2003; Liu et al., 2014; Bernt and Walker, 1999). While the functions of these fluids are well known in some cases (Liu and Kubli, 2003; Liu et al., 2014; Bernt and Walker, 1999; Perry et al., 2013), the role(s) of other socially exchanged fluids (e.g., saliva) are less clear. The context-specific transmission and interindividual confidentiality of socially exchanged fluids raise the possibility that this type of chemical exchange mediates a private means of chemical communication.
Social insects are an interesting group of animals to investigate the potential role of socially exchanged fluids. Colonies of ants, bees and termites are self-organized systems that rely on a set of simple signals to coordinate the development and behavior of individual members (Bonabeau et al., 1997). While colony-level phenotypes may arise simply from the independent behavior of individuals with similar response thresholds, many group decisions require communication between members of a colony (LeBoeuf and Grozinger, 2014). In ants, three principal means of communication have been described: pheromonal, acoustic, and tactile (Hölldobler and Wilson, 1990). Pheromones, produced by a variety of glands, impart diverse information, including nestmate identity and environmental dangers. Acoustic communication, through substrate vibration or the rubbing of specific body parts against one another, often conveys alarm signals. Tactile communication encompasses many behaviors, from allo-grooming and antennation to the grabbing and pulling of another ant’s mandibles for recruitment to a new nest site or resource.
Ants, like many social insects and some vertebrates (Boulay et al., 2000a; Greenwald et al., 2015; Malcolm and Marten, 1982; Wilkinson, 1984), also exhibit an important behavior called trophallaxis, during which liquid is passed mouth-to-mouth between adults or between adults and juveniles. The primary function of trophallaxis is considered to be the exchange of food, as exemplified by the transfer of nutrients from foragers to nurses and from nurses to larvae (Wilson and Eisner, 1957; Buffin et al., 2009; Cassill and Tschinkel, 1995, 1996; Wheeler, 1918, Wheeler, 1986). The eusocial Hymenopteran forgeut has evolved a specialized distensible crop and a restrictive proventriculus (the separation between foregut and midgut) enabling frequent fluid exchange and regulation of resource consumption (Terra, 1990; Eisner and Brown, 1958; Lanan et al., 2016; Hunt, 1991). In addition to simple nourishment, trophallaxis can provide information for outgoing foragers about available food sources (Grüter et al., 2006; Farina et al., 2007).
Trophallaxis also occurs in a number of non-food related contexts, such as reunion with a nestmate after solitary isolation (Boulay et al., 2000b), upon microbial infection (Hamilton et al., 2011), and in aggression/appeasement interactions (Liebig et al., 1997). Furthermore, adult ants have been suggested to use a combination of trophallaxis and allo-grooming to share cuticular hydrocarbons (CHCs) which are important in providing a specific ‘colony odor’ (Boulay et al., 2000a; Soroker et al., 1995), although the presence of CHCs in trophallactic fluid has not been directly demonstrated. Given these food-independent trophallaxis events, and the potential of this behavior to permit both ‘private’ inter-individual chemical exchange as well as rapid distribution of fluids over the social network of a colony, we tested the hypothesis that trophallaxis serves as an additional means of chemical communication and/or manipulation. To identify the endogenous molecules exchanged during this behavior, we used mass-spectrometry and RNA sequencing to characterize the contents of trophallactic fluid, and identified many growth-related proteins, CHCs, small RNAs, and the insect developmental regulator, juvenile hormone. We also obtained evidence that some components of trophallactic fluid can be modulated by social environment, and may influence larval development.
Results
Collection and proteomic analysis of trophallactic fluid
Analysis of the molecules exchanged during trophallaxis necessitated development of a robust method for acquiring trophallactic fluid (TF). We focused on the Florida carpenter ant, Camponotus floridanus, which is a large species whose genome has been sequenced (Bonasio et al., 2010). We first attempted to collect fluid from unmanipulated pairs of workers engaged in trophallaxis, but it was impossible to predict when trophallaxis would occur, and events were usually too brief to collect the fluid being exchanged. We found that after workers were starved and isolated from their colony, then fed a 25% sucrose solution and promptly reunited with a similarly conditioned nestmate, approximately half of such pairs displayed trophallaxis within the first minute of reunion (as observed previously [Boulay et al., 2000a; Dahbi et al., 1999]), and were more likely to remain engaged in this behavior for many seconds or even minutes. Under these conditions, it was sometimes possible to collect small quantities of fluid from the visible droplet between their mouthparts (referred to as ‘voluntary’ samples). However, even under these conditions, trophallaxis was easily interrupted, making this mode of collection extremely low-yield.
To obtain larger amounts of TF and avoid the confounding factors of social isolation and feeding status, we developed a non-lethal method to collect the contents of the crop by lightly-squeezing the abdomen of CO2-anesthetised ants (referred to as ‘forced’ samples, similar to [Hamilton et al., 2011]). This approach yielded a volume of 0.34 ± 0.27 µL (mean ± SD) of fluid per ant. To determine whether the ‘forced’ fluid collected under anesthesia was similar to the fluid collected from ants voluntarily engaged in trophallaxis, and ensure that it was not contaminated with hemolymph or midgut contents, we also collected samples of these fluid sources. To compare the identities and quantities of the different proteins found in each fluid, all samples were analyzed by nanoscale liquid chromatography coupled to tandem mass spectrometry (nano-LC-MS/MS) (Figure 1A).
Hierarchical clustering of normalized spectral counts of the proteins identified across our samples revealed high similarity between voluntary and forced TF, but a clear distinction of these fluids from midgut contents or hemolymph (Figure 1A, for protein names and IDs see Figure 1—figure supplement 1). While the voluntary TF samples contained fewer identified proteins than the forced TF samples – likely due to lower total collected TF volume per analyzed sample (< 1 µL voluntary vs. > 10 µL forced) – the most abundant proteins were present across all samples in both methods of collection (Figure 1A). To investigate whether the few differences observed between the voluntary and forced TF in the less abundant proteins might be due to the starvation and/or social isolation conditions used to collect voluntary TF, we isolated groups of 25–30 ants from their respective queens and home colonies for 14 days, with constant access to food and water, and collected TF both directly before and after the period of isolation. Social isolation affected the ratios of proteins in TF, with five of the top 40 proteins in TF becoming significantly less abundant and one more abundant when ants were socially isolated (Figure 1B; see Figure 1—figure supplement 2 for names and IDs). Three of the proteins down-regulated in social isolation were also significantly less abundant in voluntary TF samples (from socially isolated ants) relative to forced TF samples (from within-colony ants, Figure 1—figure supplement 2). Together these results provide initial evidence that the composition of this fluid is influenced by social and/or environmental experience of an ant, and support the validity of our methodology to force collect TF.
Of the 50 most abundant proteins found in TF (Figure 1C), many are likely to be digestion-related (e.g. three maltases, one amylase, various proteases, two glucose dehydrogenases, one DNase) and three Cathepsin D homologs might have immune functions (Hamilton et al., 2011). However, at least 10 of the other proteins have putative roles in the regulation of growth and development, including two hexamerins (nutrient storage proteins [Zhou et al., 2007]), a yellow/major royal jelly protein (MRJP) homolog (likely nutrient storage [Drapeau et al., 2006]) and an apolipophorin (vitellogenin-domain containing lipid-transport protein [Kutty et al., 1996]). Most notably, TF contained several proteins that are homologous to insect juvenile hormone esterases (JHEs) and Est-6 in Drosophila melanogaster (Figure 1D). JHEs are a class of carboxylesterases that degrade the developmental regulator, juvenile hormone (JH) (Kamita and Hammock, 2010; Nijhout et al., 2014), and Est-6 is an abundant esterase in D. melanogaster seminal fluid (Mane et al., 1983; Richmond et al., 1980; Younus et al., 2014). While most insects have only one or two JHE-like proteins (Nijhout et al., 2014; Qu, 2015), C. floridanus has an expansion of more than 10 related proteins, eight of which were detected in TF (Figure 1D, Figure 1—figure supplements 1 and 2) and two in hemolymph (Figure 1—figure supplement 1). Three of the abundant JHE-like proteins are significantly down-regulated in the TF of ants that have undergone social isolation (Figure 1—figure supplement 2).
Thirty-three of the 50 most abundant TF proteins had predicted N-terminal signal peptides, suggesting that they can be secreted directly into this fluid by cells lining the lumen or glands connected to the alimentary canal (Figure 5—figure supplement 1). Moreover, half of the proteins without such a secretion signal had gene ontology terms indicating extracellular or lipid-particle localization, which suggests that they may gain access to the lumen of the foregut through other transport pathways.
Trophallaxis microRNAs
Many small RNAs have been found in externally secreted fluids across taxa, such as seminal fluid, saliva, milk and royal jelly (Sarkies and Miska, 2013; Weber et al., 2010; Guo et al., 2013). Although the functions of extracellular RNAs remain unclear (Turchinovich et al., 2012), we investigated if TF also contains such molecules by isolating and sequencing small RNAs from C. floridanus TF. After filtering out RNA corresponding to potential commensal microorganisms (including the known symbiont, Blochmannia floridanus; Supplementary file 1) and other organic food components, we detected 64 miRNAs. Forty-six of these were identified based upon their homology to miRNAs of the honey bee Apis mellifera (Guo et al., 2013; Søvik et al., 2015; Greenberg et al., 2012), while 18 sequences (bearing the structural stem-loop hallmarks of miRNA transcripts) were specific to C. floridanus (Figure 2, Figure 2—source data 1). The most abundant of the 64 miRNAs was miR-750, followed by three C. floridanus-specific microRNAs. The role of miR-750 is unknown, but the expression of an orthologous miRNA in the Asian tiger shrimp, Penaeus monodon, is regulated by immune stress (Kaewkascholkul et al., 2016). Notably, 16 of the miRNAs detected in the C. floridanus TF are also present in A. mellifera worker jelly and/or royal jelly (Guo et al., 2013), which are oral secretions fed to larvae to bias them toward a worker or queen fate.
Trophallactic fluid contains long-chain hydrocarbons
Trophallaxis has long been suggested to contribute to the exchange and homogenization of colony odor (Crozier and Dix, 1979). The observation that different ant species have distinct blends of non-volatile, cuticular hydrocarbons (CHCs) – which also display quantitative variation within species – have made these chemicals prime candidates for conveying nestmate recognition cues (Sharma et al., 2015; Bos et al., 2011; Hefetz, 2007; Lahav et al., 1999; Ozaki et al., 2005; van Zweden and d’Ettorre, 2010). CHC profiles of ants within a colony are similar, while individuals isolated from their home colony often vary in their CHC make-up (Boulay and Lenoir, 2001), suggesting that the profiles are constantly unified between nestmates (Boulay et al., 2000a). However, it remains unclear if CHCs are exchanged principally by trophallaxis, or through other mechanisms such as allo-grooming and contact with the nest substrate (Boulay et al., 2000a; Soroker et al., 1995; Bos et al., 2011). We therefore investigated whether CHCs are present in TF, and how they relate to those present on the cuticle. To do this, we analyzed by gas chromatography mass spectrometry (GC-MS) the TF and cuticular extracts collected from the same five groups of 20–38 workers.
We identified 63 molecules in TF (Table 1), including eight fatty acids and fatty acid esters, 13 unbranched and 36 branched hydrocarbons with one to five methyl branches. The majority (40/63) of the TF compounds comprised 27 or more carbons. Cuticular extracts also contained predominantly multiply branched alkanes with 27 or more carbons, corresponding to previous observations of CHCs in this species (Sharma et al., 2015; Endler et al., 2004) (Figure 3). All the highly abundant hydrocarbons in TF (marked in Figure 3) were also present on the cuticle of the workers analyzed, supporting the potential for trophallaxis to mediate inter-individual CHC exchange. Moreover, there was a significantly greater (t-test, p<0.0003) similarity across colonies in the hydrocarbon profiles of TF than hydrocarbon profiles of the cuticle (similarity was measured by cross-correlation between GC-MS profiles, Figure 3D–E). Altogether, these observations indicate that while CHCs are likely to be exchanged by trophallaxis, additional mechanisms are probably involved in generating the colony-specific bouquet of these compounds.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.20375.008
-
Figure 3—source data 2
- https://doi.org/10.7554/eLife.20375.009
Juvenile hormone is exchanged through trophallaxis and can influence larval development
Given the abundance of the expanded family of JHE-like proteins in TF, we asked whether JH itself is also present in this fluid. The primary JH found in Hymenoptera, JH III, is thought to circulate in the hemolymph after being produced by the corpora allata (Wigglesworth, 1936). To detect and quantify JH in both TF and hemolymph, we employed a derivatization and purification process prior to GC-MS analysis (Brent and Vargo, 2003; Shu et al., 1997; Bergot et al., 1981). While both sets of measurements were highly variable across samples, we found high levels of JH in TF with concentrations of the same order of magnitude as those found in hemolymph (Figure 4A).
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.20375.013
-
Figure 4—source data 2
- https://doi.org/10.7554/eLife.20375.014
-
Figure 4—source data 3
- https://doi.org/10.7554/eLife.20375.015
-
Figure 4—source data 4
- https://doi.org/10.7554/eLife.20375.016
Because JH is an important regulator of development, reproduction and behavior (Nijhout and Wheeler, 1982), we investigated whether the dose of JH that a larva receives by trophallaxis with a nursing worker might be physiologically relevant. If an average worker has approximately 0.34 μL of TF in her crop at a given time (as measured in our initial experiments; see above), this corresponds to a dose of approximately 31 pg of JH. The analysis of 35 larvae collected midway through development (i.e., third instar out of four worker instars, mean ± SD: 4.0 ± 0.19 mm long and 1.4 ± 0.12 mm wide) revealed that they contained 100–700 pg of JH (Figure 4B). Thus, the amount of JH received during an average trophallaxis event amounts to 5–31% of the JH content of a third instar larva. While it is difficult to determine how much of the JH fed to larvae remains in the larval digestive tract, these results indicate that there is potentially sufficient JH in a single trophallaxis-mediated feeding to shift the titer of a recipient larva.
We next determined whether adding exogenous JH to the food of nursing workers could change the growth of the reared larvae. We created groups of 25–30 workers and allowed them to each rear 5–10 larvae to pupation, in the presence of food and sucrose solution that was supplemented with either JH or only a solvent. Larvae reared by JH-supplemented workers grew into larger adults than those reared by solvent-supplemented controls (Figure 4C, GLMM, p<9.01e−06). Moreover, larvae reared by JH-supplemented caretakers were twice as likely to successfully undergo metamorphosis relative to controls (Figure 4D, binomial GLMM, p<7.39e−06). These results are consistent with previous studies in other species of ants and bees using methoprene (a non-hydrolyzable JH analog), whose external provision to the colony can lead to larvae developing into larger workers and even queens (Nijhout and Wheeler, 1982; Libbrecht et al., 2013; Wheeler and Nijhout, 1981).
Comparative proteomics reveals species-specific growth-regulatory proteins in trophallactic fluid
To expand our survey of TF, we collected this fluid from other species of social insects: a closely related ant (C. fellah), an ant from another sub-family (the fire ant Solenopsis invicta) and the honey bee A. mellifera. Nano-LC-MS/MS analyses identified 79, 350 and 136 proteins in these three species, respectively (84 were identified in C. floridanus). We assigned TF proteins from all four analyzed species to 138 distinct groups of predicted orthologous proteins (see Materials and methods); of these, 72 proteins were found in the TF of only one species (Figure 5, Figure 5—source data 1).
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.20375.018
Only eight ortholog groups contained representatives present in the TF of all four species (Figure 5B). Most of these appear to have functions related to digestion, except for apolipophorin, which is involved in lipid/nutrient transport (Kutty et al., 1996). For ortholog groups found in the TF of the three ant species, most were also digestion-related with the notable exception of CREG1, a secreted glycoprotein that has been implicated in cell growth control (Di Bacco et al., 2003) and insect JH response (Li Y et al., 2007; Barchuk et al., 2007; Zhang et al., 2014). Genus- and species-specific proteins were frequently associated with growth or developmental roles. For example, A. mellifera TF contained 12 distinct major royal jelly proteins (MRJP), which are thought to be involved in nutrient storage and developmental fate determination (Drapeau et al., 2006; Kamakura, 2011), while S. invicta TF contained a highly abundant JH-binding protein and a vitellogenin. The two Camponotus species shared many orthologous groups, consistent with their close phylogenetic relationship. Five of the seven JHE/Est-6 proteins found by Scaffold in C. floridanus TF also had orthologs present in C. fellah TF. Additionally, these two species shared a MRJP/Yellow homolog, and an NPC2-related protein, which, in D. melanogaster, is involved in sterol binding and ecdysteroid biosynthesis (Huang et al., 2007). Finally, 26 ortholog groups contained representatives from multiple species, but not the most closely related ones (e.g., A. mellifera and S. invicta, or S. invicta and only one of the two Camponotus species). One-third of these were associated with growth and developmental processes (e.g., three hexamerins, two MRJP/yellow proteins, imaginal disc growth factor 4, vitellogenin-like Vhdl and an additional NPC2). Together these analyses indicate that approximately half of all TF protein ortholog groups appear to be digestion-related, consistent with TF being composed of the contents of the foregut. However, many TF proteins have putative roles in growth, nutrient storage, or the metabolism and transport of JH, vitellogenin or ecdysone.
Discussion
We have characterized the fluid that is orally exchanged during trophallaxis, a distinctive behavior of eusocial insects generally considered as a means of food sharing. Our results reveal that the transmitted liquid contains much more than food and digestive enzymes, and includes non-proteinaceous and proteinaceous molecules implicated in chemical discrimination of nestmates, growth and development, and behavioral maturation. These findings suggest that trophallaxis underlies a private communication mechanism that can have multiple phenotypic consequences. More generally, our observations open the possibility that exchange of oral fluids (e.g., saliva) in other animals might also serve functions not previously suspected (Humphrey and Williamson, 2001; Ribeiro et al., 1995).
In ants, trophallaxis has long been thought to be a mode of transfer for the long-chain hydrocarbons that underlie nestmate recognition (Boulay et al., 2000a Dahbi et al., 1999; Lahav et al., 1999; van Zweden and d’Ettorre, 2010; Soroker et al., 1995; Lalzar et al., 2010). However, previous work has analyzed only the passage of radio-labeled hydrocarbons between individuals (Boulay et al., 2000a; Dahbi et al., 1999; Soroker et al., 1995) and between individuals and substrates (Bos et al., 2011). Without explicitly sampling the contents of the crop, it is not possible to differentiate between components passed by trophallaxis or by physical contact. Our study is, to our knowledge, the first to demonstrate that TF contains endogenous long-chain hydrocarbons. Interestingly, we found differences in the hydrocarbon profiles of the TF and the cuticle of the same ants, suggesting the existence of other processes regulating the relative proportion of CHCs and/or additional sources of these compounds. Further work manipulating and comparing hydrocarbon profiles of TF and cuticle of multiple individual ants within the same colony is necessary to understand the impact of the observed variation in TF and cuticular hydrocarbons.
Consistent with the view that the gut is one of the first lines of defense in the body’s interface with the outside world (Lemaitre and Hoffmann, 2007; Söderhäll and Cerenius, 1998), the four TF proteomes analyzed in this study include many potential defense-related proteins. Homologs of the Cathepsin D family of proteins, which have been implicated in growth, defense and digestion (Saikhedkar et al., 2015; Ahn and Zhu-Salzman, 2009), were present in the TF of both Camponotus species. Prophenoloxidase 2, an enzyme responsible for the melanization involved in the insect immune response (Söderhäll and Cerenius, 1998), was found in the TF of A. mellifera and S. invicta. All four species contained several members of the serine protease and serpin families, which have well-documented roles in the prophenoloxidase cascade and more general immune responses in D. melanogaster and other animals (Lemaitre and Hoffmann, 2007; Lucas et al., 2009). The TF of two species had orthologs of the recognition lectin GNBP3 which is also involved in the prophenoloxidase cascade. Finally, a handful chromatin-related proteins (e.g., histones, CAF1) were found in S. invicta, and to a lesser extent in the other two ants. Their presence may simply reflect a few cells being sloughed from the lumen of the foregut, or could be indicative of a defense process termed ETosis, whereby chromatin is released from the nuclei of inflammatory cells to form extracellular traps that kill pathogenic microbes (Robb et al., 2014; Brinkmann et al., 2004).
The presence of development- and growth-related components in the TF of diverse social insects suggests that this fluid may play a role in directing larval development. Previous work has shown that the developmental fate of larvae and the process of caste determination in various ant and bee species can be influenced by treating colony members with JH or JH analogs (Wheeler, 1986; Wheeler and Nijhout, 1981; Cnaani et al., 1997; Cnaani et al., 2000; Rajakumar et al., 2012). Moreover, in some species, workers play an important role in regulating caste determination (Kamakura, 2011; Linksvayer et al., 2011; Villalta et al., 2016; Mutti et al., 2011) but it is unknown how workers might influence the JH titers of larvae. Our finding that TF contains JH raises the possibility of trophallaxis as a direct means by which larval hormone levels and developmental trajectories can be manipulated. This type of mechanism has some precedent: honey bee workers bias larval development toward queens by feeding them royal jelly, an effect that might be mediated by MRJPs (Kamakura, 2011; Buttstedt et al., 2016; Kucharski et al., 2015; Kamakura, 2016). Across the TF of four species of social insect, we have found diverse molecules intimately involved in insect growth regulation: JH, JHE, JH-binding protein, vitellogenin, hexamerin, apolipophorin, and MRJPs. Several of the proteins and microRNAs identified in ant and bee TF are also found in royal and worker jelly (Guo et al., 2013; Zhang et al., 2014). Notably, there are also some similarities between proteins in social insect TF and mammalian milk (Zhang et al., 2015), such as the cell growth regulator CREG1. The between-genera variation in the most abundant growth-related proteins (e.g., MRJPs in honey bee, JH-binding protein and hexamerins in fire ant, and JHEs in Camponotus) indicates that there might be multiple evolutionary origins and/or rapid divergence in trophallaxis-based signals potentially influencing larval growth.
The finding that TF contains many proteins, miRNAs, CHCs and JH raises the question of how they come to be present in this fluid. Many non-food-derived components in TF are likely to be either directly deposited into the alimentary canal: over two-thirds of the proteins present in TF had a predicted signal peptide suggesting that they are probably secreted by the cells lining the foregut or glands connected to the alimentary canal (e.g., postpharyngeal, labial, mandibular, salivary glands). Additionally, TF molecules are likely to be acquired by transfer among nestmates through licking, grooming, and trophallaxis. For example, CHCs are produced by oenocytes, transmitted to the cuticle, ingested, and sequestered in the postpharyngeal gland (Bos et al., 2011; van Zweden and d’Ettorre, 2010; Soroker and Hefetz, 2000). Interestingly, JH is synthesized in the corpora allata, just caudal to the brain (Goodman and Granger, 2009). Although the mechanisms of JH uptake by tissues are still incompletely understood (Parra-Peralbo and Culi, 2011; Rodríguez-Vázquez et al., 2015; Engelmann and Mala, 2000; Suzuki et al., 2011), the TF proteins apolipophorin, hexamerin, JH-binding protein and vitellogenin have all been implicated in JH binding and transport between hemolymph and tissues (Goodman and Granger, 2009; Rodríguez-Vázquez et al., 2015; Engelmann and Mala, 2000; Suzuki et al., 2011; Amsalem et al., 2014), and these factors may be responsible for transporting JH into the foregut. Unfortunately, even in Drosophila mechanisms of JH uptake by tissues are still unclear (Rodríguez-Vázquez et al., 2015; Engelmann and Mala, 2000; Suzuki et al., 2011; Parra-Peralbo and Culi, 2011). Extracellular miRNAs are secreted and transported through a variety of pathways, but the functional relevance of such molecules is still controversial (Sarkies and Miska, 2013; Turchinovich et al., 2016; Masood et al., 2016; Søvik et al., 2015; Rayner and Hennessy, 2013).
While the simultaneous presence of JH and a set of putative JH degrading enzymes (the JHE/Est-6 family) in C. floridanus TF may appear surprising, there are at least two possible explanations. First, this might reflect a regulatory mechanism for JH levels at the individual and colony levels. Many biological systems use such negative feedback mechanisms to buffer signals and enable rapid responses to environmental change (Alon, 2007). Alternatively, these enzymes may have evolved a different function along with their novel localization in TF. In other insects JHE-related enzymes are typically expressed in the fat body and circulate in the hemolymph (Hammock et al., 1975; Ward et al., 1992; Campbell et al., 1998)., The expansion of JHE-like proteins in C. floridanus was accompanied by a high specificity in their localization, with two (E2AM67, E2AM68) being present exclusively in the hemolymph, and four (E2ANU0, E2AI90, E2AJM0, E2AJL9) exclusively in TF. This coexistence of JH and JHEs reveals a striking parallel with the constituents of a different socially exchanged fluid in D. melanogaster. Drosophila Est-6 is highly abundant in seminal fluid, together with its presumed substrate, the male-specific pheromone 11-cis-vaccenyl acetate (cVA). cVA has multiple roles in influencing sexual and aggressive interactions and its transfer to females during copulation potently diminishes her attractiveness to future potential mates (Poiani, 2006; Chertemps et al., 2012; Costa, 1989).
Given that larvae are fed JH and food, and both are necessary for development, a future goal will be to dissect the relative contribution and potential synergy of these components. Furthermore, JH has many other functions in social insects, including behavioral modulation, longevity, fecundity and immunity (Flatt et al., 2005; Jindra et al., 2013; Dolezal et al., 2012; Wang et al., 2012); the relevance of trophallaxis-mediated JH exchange between adults remains to be explored. Considering that most, if not all, individuals in a colony must share food through trophallaxis, it will also be of interest to understand if and how individual-, or caste-specific TF-based information signals emerge. Recent work on trophallaxis networks (Greenwald et al., 2015) and interaction networks (Mersch et al., 2013), indicate that ants preferentially interact with others of similar behavioral type (e.g., nurses with nurses, and foragers with foragers). This raises the possibility that different pools of TF with different qualities exist. A key challenge will be to develop specific genetic and/or pharmacological tools to test the biological relevance of molecules transmitted by trophallaxis.
Materials and methods
Insect source and rearing
Request a detailed protocolCamponotus floridanus workers came from 16 colonies established in the laboratory from approximately 1-year-old founding queens and associated workers collected from the Florida Keys in 2006, 2011 and 2012. Ants were provided once a week with fresh sugar water, an artificial diet of honey or maple syrup, eggs, agar, canned tuna and a few D. melanogaster. Maple syrup was substituted for honey and no Drosophila were provided in development and proteomic experiments to avoid contamination with other insect proteins. Colonies were maintained at 25°C with 60% relative humidity and a 12 hr light:12 hr dark cycle. Camponotus fellah colonies were established from queens collected after a mating flight in March 2007 in Tel Aviv, Israel (Colonies #5, 28, 33). The ant colonies were maintained at 32°C with 60% relative humidity and a 12 hr light:12 hr dark cycle. Fire ant workers (Solenopsis invicta) were collected from three different colonies, two polygyne, one monogyne, maintained at 32°C with 60% relative humidity and a 12 hr light:12 hr dark cycle. Honeybee workers (Apis mellifera, Carnica and Buckfast) were collected from six different hives maintained with standard beekeeping practices.
Fluid collection
Request a detailed protocol‘Voluntary’ TF was collected from individuals who had been starved and socially isolated for 1–3 weeks, then fed 25% sucrose solution, and promptly re-introduced to other separately isolated nestmates. Ants were monitored closely for trophallaxis events; when one had begun, a pulled glass pipette was brought between the mouthparts of the 2–4 individuals and fluid was collected. Typically this stopped the trophallaxis event, but on rare occasions it was possible to acquire up to sub-microliter volumes of TF.
For ‘forced’ collection, ants were anesthetized by CO2 (on a CO2 pad; Flypad, FlyStuff, San Diego, CA) and promptly flipped ventral-side up. The abdomen of each ant was lightly squeezed with wide insect forceps to prompt the regurgitation of fluids from the social stomach. Ants that underwent anesthesia and light squeezing yielded on average 0.34 µL of fluid and recovered in approximately 5 min. TF was collected with graduated borosilicate glass pipettes, and transferred immediately to either buffer (100 mM NaH2PO4, 1 mM EDTA, 1 mM DTT, pH 7.4) for proteomic experiments, to RNase-free water for RNA analyses, or to pure EtOH for GC-MS measurements. TF was collected from C. fellah and S. invicta in the same manner as from C. floridanus; for S. invicta 100–300 ants were used due to their smaller body size and crop volume.
Honey bee TF was collected from bees that were first cold-anesthetized and then transferred to a CO2 pad to ensure continual anesthesia during collection of TF as described above. Honey bees yielded higher volumes of TF than did ants (0.94 µL TF per honey bee, SD = 0.54 relative to 0.34 µL TF per C. floridanus ant, SD = 0.27).
Hemolymph was collected from CO2-anesthetized ants by puncturing the junction between the foreleg and the distal edge of the thoracic plate with a pulled glass pipette. This position was chosen over the abdomen in order to ensure that hemolymph was collected and not the contents of the crop. Approximately 0.1 µL was taken from each ant.
Midgut samples were collected by first anesthetizing an ant, immobilizing it in warm wax ventral-side up, covering the preparation in 1x PBS, and opening the abdomen with dissection forceps and iris scissors. The midgut was punctured by a sharp glass pipette and its contents collected. Because the pipette also contacts the surrounding fluid, some hemolymph contamination was unavoidable.
Proteomic analyses
Gel separation and protein digestion
Request a detailed protocolProtein samples were loaded on a 12% mini polyacrylamide gel and migrated about 2.5 cm. After Coomassie staining, regions of gel lanes containing visible bands (generally > 18 kDa) were excised into 2–4 pieces, depending on the gel pattern of considered experiment. Gel pieces were digested with sequencing-grade trypsin (Promega, Switzerland) as described (Shevchenko et al., 2006). Extracted tryptic peptides were dried and resuspended in 0.1% formic acid, 2% (v/v) acetonitrile for mass spectrometry analyses.
Proteomic mass spectrometry analyses
Request a detailed protocolTryptic peptide mixtures were injected on a Dionex RSLC 3000 nanoHPLC system (Dionex, Sunnyvale, CA) interfaced via a nanospray source to a high resolution mass spectrometer based on Orbitrap technology (Thermo Fisher, Bremen, Germany): LTQ-Orbitrap XL (‘voluntary’ and hemolymph samples), LTQ-Orbitrap Velos (A. mellifera TF samples) or QExactive Plus (‘voluntary’ and all other samples). Peptides were loaded onto a trapping microcolumn Acclaim PepMap100 C18 (20 mm x 100 μm ID, 5 μm, Dionex) before separation on a C18 reversed-phase analytical nanocolumn at a flowrate of 0.3 μL/min. Q-Exactive Plus instrument was interfaced with an Easy Spray C18 PepMap nanocolumn (25 or 50 cm x 75 µm ID, 2 µm, 100 Å, Dionex) using a 35-min gradient from 4% to 76% acetonitrile in 0.1% formic acid for peptide separation (total time: 65 min). Full MS surveys were performed at a resolution of 70,000 scans. In data-dependent acquisition controlled by Xcalibur software (Thermo Fisher), the 10 most intense multiply charged precursor ions detected in the full MS survey scan were selected for higher energy collision-induced dissociation (HCD, normalized collision energy NCE = 27%) and analyzed in the orbitrap at 17,500 resolution. The window for precursor isolation was of 1.5 m/z units around the precursor and selected fragments were excluded for 60 s from further analysis.
The LTQ-Orbitrap Velos mass spectrometer was interfaced with a reversed-phase C18 Acclaim Pepmap nanocolumn (75 μm ID x 25 cm, 2.0 μm, 100 Å, Dionex) using a 65 min gradient from 5% to 72% acetonitrile in 0.1% formic acid for peptide separation (total time: 95 min). Full MS surveys were performed at a resolution of 60,000 scans. In data-dependent acquisition controlled by Xcalibur software, the 20 most intense multiply charged precursor ions detected in the full MS survey scan were selected for CID fragmentation (NCE = 35%) in the LTQ linear trap with an isolation window of 3.0 m/z and then dynamically excluded from further selection for 120 s.
LTQ-Orbitrap XL instrument was interfaced with a reversed-phase C18 Acclaim Pepmap (75 μm ID x 25 cm, 2.0 μm, 100 Å, Dionex) or Nikkyo (75 μm ID x 15 cm, 3.0 μm, 120 Å, Nikkyo Technos, Tokyo, Japan) nanocolumn using a 90 min gradient from 4% to 76% acetonitrile in 0.1% formic acid for peptide separation (total time: 125 min). Full MS surveys were performed at a resolution of 60,000 scans. In data-dependent acquisition controlled by Xcalibur software, the 10 most intense multiply charged precursor ions detected in the full MS survey scan were selected for CID fragmentation (NCE = 35%) in the LTQ linear trap with an isolation window of 4.0 m/z and then dynamically excluded from further selection for 60 s.
Proteomic data analysis
Request a detailed protocolMS data were analyzed primarily using Mascot 2.5 (RRID:SCR_014322, Matrix Science, London, UK) set up to search either the UniProt (RRID:SCR_002380, www.uniprot.org) or NCBI (RRID:SCR_003496, www.ncbi.nlm.nih.gov) database restricted to C. floridanus (UniProt, August 2014 version: 14,801 sequences; NCBI, July 2015 version: 34,390 sequences), S. invicta (UniProt, August 2014 version: 14,374 sequences; NCBI, January 2015 version: 21,118 sequences) or A. mellifera (UniProt, September 2015 version: 15,323 sequences; NCBI, February 2016 version: 21,772 sequences) taxonomy. For Mascot search of C. fellah samples, we used the database provided by the C. fellah transcriptome (27,062 sequences) and the UniProt C. floridanus reference proteome (January 2016 version, 14,287 sequences). Trypsin (cleavage at K, R) was used as the enzyme definition, allowing two missed cleavages. Mascot was searched with a parent ion tolerance of 10 ppm and a fragment ion mass tolerance of 0.50 Da (LTQ-Orbitrap Velos / LTQ-Orbitrap XL) or 0.02 Da (QExactive Plus). The iodoacetamide derivative of cysteine was specified in Mascot as a fixed modification. N-terminal acetylation of protein, deamidation of asparagine and glutamine, and oxidation of methionine were specified as variable modifications.
Scaffold software (version 4.4, RRID:SCR_014345, Proteome Software Inc., Portland, OR) was used to validate MS/MS based peptide and protein identifications, and to perform dataset alignment. Peptide identifications were accepted if they could be established at greater than 90.0% probability as specified by the Peptide Prophet algorithm (Keller et al., 2002) with Scaffold delta-mass correction. Protein identifications were accepted if they could be established at greater than 95.0% probability and contained at least two identified peptides. Protein probabilities were assigned by the Protein Prophet algorithm (Nesvizhskii et al., 2003). Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony and in most cases these proteins appear identical in annotations.
Quantitative spectral counting was performed using the normalized spectral abundance factor (NSAF), a measure that takes into account the number of spectra matched to each protein, the length of each protein and the total number of proteins in the input database (Zybailov et al., 2006). To compare relative abundance across samples with notably different protein abundance, we divided each protein’s NSAF value by the total sum of NSAF values present in the sample.
Across forced TF samples, ~10% of detected spectra were identified as peptides by Mascot and Scaffold, while in hemolymph samples ~16% of spectra were identified. Of those identified protein spectra in TF, about 2/3 mapped to the C. floridanus proteome, 1/3 to typical proteomic experiment contaminants (e.g., human keratin, trypsin added for digestion of proteins), < 5% to elements from the ants’ diet (egg proteins, D.melanogaster), and < 0.5% to microbial components including endosymbiont Blochmannia floridanus. Of the identified hemolymph protein spectra, 94% mapped to the C. floridanus proteome and the remaining 6% to keratins, trypsins, etc. We used PEAKS proteomic software (version 7.5, Bioinformatics Solutions Inc., Waterloo, Canada) for a more in-depth analysis of a typical TF sample, making less stringent de novo protein predictions allowing substitutions and post-translational modifications, and we could identify approximately 50% of spectra as peptides.
All proteomics data are available through ProteomeXchange at PXD004825.
Secretion signals were predicted using SignalP 4.1 (Petersen et al., 2011).
Social isolation
View detailed protocolGroups of 25–30 C. floridanus worker ants were taken from colonies C3, C6, C7, C8, C9 and C10 (queens all collected in the Florida Keys December 2012). Upon collection, TF was collected under CO2 anesthesia. Immediately after TF collection ants were collectively isolated from their home colony and queen in fluon-coated plastic boxes with shelter, insect-free food (maple syrup, chicken eggs, tuna and agar), and water. Ants were kept in these group-isolated queenless conditions for 14 days, after which, their TF was collected again.
Small RNA analysis
Total RNA isolation and quantification
Request a detailed protocolTotal RNA was isolated (Sapetschnig et al., 2015) from approximately 100 µL of TF from C. floridanus ants by the Trizol-LS method according to the manufacturer’s protocols (Life Technologies, Inc., Grand Island, NY). A Qubit RNA HS Assay Kit (high sensitivity, 5 to 100 ng) was used with a Qubit 3.0 fluorometer according to the manufacturer’s protocols (Life Technologies, Carlsbad); a sample volume of 1 μL was added to 199 μL of a Qubit working solution. The RNA concentration was 17.7 ng/µL.
Library preparation
Request a detailed protocolWe used 5 µL of RNA (approximately 88.5 ng) to generate a small RNA sequencing library using reagents and methods provided with TruSeq Small RNA Sample Prep Kit (Illumina, San Diego). Briefly, T4 RNA ligase was used to ligate RA5 and RA3 RNA oligonucleotides to 5' and 3' ends of RNA, respectively. Adapter-ligated RNA was reverse-transcribed using a RTP primer and the resulting cDNA was amplified in an 11-cycle PCR that used RP1 and indexed RP1 primers. We size-selected cDNA libraries using 6% TBE PAGE gels (Life Technologies, Carlsbad, CA) and ethidium bromide staining. Quality of the generated small RNA sequencing library was confirmed using electropherograms from a 2200 TapeStation System (Agilent Technologies, Santa Clara, CA). Desired sizes of cDNA bands were cut from the gel (between 147 and 157 nt), the gel matrix broken by centrifugation through gel breaker tubes (IST Engineering Inc., Milpitas, CA), and cDNA eluted with 400 μL of 0.3M Na-Chloride. Further purification of cDNA was by centrifugation through Spin-X 0.22 μm cellulose acetate filter columns (Corning Costar, Corning, NY) followed by ethanol precipitation. Libraries were sequenced on a HiSeq 2000 Sequencer (Illumina). Small RNA sequencing data are available through SRA database at SRP082161.
Pre-processing of RNA libraries
Request a detailed protocolRaw sequenced reads from small RNA libraries were submitted to quality filtering and adaptor trimming using cutadapt (version 1.8.1, RRID:SCR_011841, https://pypi.python.org/pypi/cutadapt/1.8.1). Small RNA reads with Phred quality below 20, and fewer than 18 nucleotides after trimming of adaptors, were discarded.
MicroRNA prediction using miRDeep2
Request a detailed protocolRemaining sequences from small RNA libraries were used with A. melifera microRNA database (miRbase version 21) and C. floridanus genome (NCBI version 1.0) to predict microRNAs precursors using miRDeep2 (RRID:SCR_012960) with default parameters and GFF (General Feature Format) file extracted by Perl scripts.
Small RNA reads mapped to C. floridanus and D. melanogaster genomes
Request a detailed protocolApproximately half the reads that mapped to the C. floridanus genome also mapped to D. melanogaster, a component of the ants’ laboratory diet. While some of these RNAs are likely to be endogenous C. floridanus RNAs, we eliminated all reads identical between C. floridanus and D. melanogaster. Remaining sequences from small RNA libraries were mapped to reference sequences from C. floridanus and D. melanogaster genomes using Bowtie (version 1.1.1, RRID:SCR_005476, one mismatch allowed). The C. floridanus genome (version 1.0) was downloaded from NCBI. The D. melanogaster genome (version v5.44) was downloaded from flybase.org. Remaining sequenced reads that did not map to the C. floridanus or D. melanogaster genome were assembled into contigs with VelvetOptimiser (version 2.2.5; http://bioinformatics.net.au/software.velvetoptimiser.shtml), and BLASTed against non-redundant NCBI bacterial and viral databases to assess their source organism. Hits with an E-value smaller than 1e−5 for nucleotide comparison were considered significant.
Automatic annotation, penalization, size distribution and gene expression
Request a detailed protocolTo perform automatic annotation, we used BedTools (version v2.17.0, RRID:SCR_006646) to compare genomic coordinates from mapped reads against predicted microRNAs, mRNA, tRNA and ncRNA (represented by lncRNA). Reads mapping to both the C. floridanus and D. melanogaster genomes were deemed ambiguous and were eliminated. The remaining reads were normalized by reads per million (RPM). The gene expression was measured and normalized by RPM and plotted in heatmap and barplot graphs.
Gas chromatography mass spectrometry and related sample preparation
Request a detailed protocolHydrocarbon analysis was performed on trophallactic fluid from five groups of 20–38 ants, each collected from one of five different colonies (C11-C15). TF samples were placed directly into 3:1 hexane:methanol. Immediately after TF collection, body surface CHCs were collected by placing anesthetized ants in hexane for 1 min before removal with cleaned forceps. Methanol was added to the cuticular-extraction hexane (maintaining the 3:1 proportion of the TF samples). The TF and body samples were vortexed for 30 s and centrifuged for 7 min. Hexane fractions were collected using a thrice-washed Hamilton syringe. Samples were kept at −20°C until further analysis.
A Trace 1300 GC chromatograph interfaced with a TSQ 8000 Evo Triple Quadrupole Mass Spectrometer (Thermo Scientific, Bremen, Germany) was used for the study. Hydrocarbons were separated on a 30 m x 0.25 mm I.D. (0.25 mm film thickness) Zebron ZB-5 ms capillary column (Phenomenex, Torrance, CA) using the following program: initial temperature 70°C held for 1 min, ramped to 210°C at 8 °C/min, ramped to 250°C at 2 °C/min, ramped to 300°C at 8 °C/min and held for 5 min. Helium was used as carrier gas at a constant flow of 1 mL/min. Injections of 1 µL of ants’ TF or body extracts were made using splitless mode. The injection port and transfer line temperature were kept at 250°C, and the ion source temperature set at 200°C. Ionization was done by electron-impact (EI, 70 eV) and acquisition performed in full scan mode in the mass range 50–550 m/z (scan time 0.2 s). Identification of hydrocarbons was done using XCalibur and NIST 14 library. The TIC MS was integrated and the area reported as a function of Retention time (Rt, min) for each peak.
Identification of trophallactic fluid hydrocarbons
Request a detailed protocolCharacterization of branched alkanes by GC-MS remains a challenge due to the similarity of their electron impact (EI) mass spectra and the paucity of corresponding spectra listed in EI mass spectra databases. A typical GC-MS chromatogram (Figure 3—figure supplement 1A) reveals the complexity of the TF sample. The workflow described here was systematically used to characterize the linear and branched hydrocarbons present in TF samples summarized in Table 1. The parent ion was first determined for each peak after background subtraction. Ambiguities remained in some cases due to the low intensity or absence of the molecular ion. In order to determine the number of ramifications for alkanes and alkenes, we examined the distribution of intensities of all fragment ions in the spectrum. From the experimental mass spectrum (Figure 3—figure supplement 1B–C), the intensity of all fragment ions was extracted and fitted with an exponential decay function (Figure 3—figure supplement 1D–E). Figure 3—figure supplement 1 clearly shows two different resulting EI mass spectra profiles: on the left a linear alkane corresponding to n-hexacosane (Rt = 34.69 min and MF C26H54); on the right, a ramified structure most likely corresponding to 9,20-Dimethyl nonacosane (Rt = 43.29 min and MF C31H64) with fragment ions 141 (loss of 295) and 309 (loss of 127) from the molecular ion at m/z 436. This demonstrates that extracting and fitting the fragment profile provides a reliable method, not only to discriminate between linear and branched hydrocarbons, but also to determine the number of ramifications. Once the number of ramifications and the parent ion mass were known, the position of the different ramifications could be deduced.
To determine the positions of the ramifications, we next estimated the Kovats retention index (RI) for each proposed structure. Since most branched alkanes are not present in the main libraries of EI mass spectra, a predicted RI was determined depending on: (i) the number of carbons and (ii) the number of ramifications. Based on RI values (for similar GC stationary phase) given in the NIST Chemistry WebBook (Linstrom et al., 2000) for known alkanes from C15 to C38, six different curves of RI vs. number of carbons were constructed for each class of alkane, from linear to pentamethylated alkanes (Figure 3—figure supplement 1F). For the same number of carbons and ramifications, an average value of RI was calculated. For example, the RI value given for C25H52 and two ramifications (RI = 2409) is the average of 3,(7/9/11/13)-dimethyltricosane (RI = 2405), 3,(3/5)-dimethyltricosane (RI = 2444) and 5,(9/11)-dimethyltricosane (RI = 2380) values listed in the NIST Chemistry WebBook. The curve corresponding to linear alkanes was strictly superimposed with the one obtained in our conditions with a standard mixture of C8-C40 alkanes. Those curves were used to determine RI values for all proposed hydrocarbons in TF samples.
As a final validation for the class of hydrocarbon proposed, the experimental retention time of each compound was plotted as a function of the calculated RI index (Figure 3—figure supplement 1G). The experimental plot was overlaid with values we recorded for linear alkane standards C8-C40. For each proposed hydrocarbon, the Kovats index should fit this curve.
The table in Table 1 summarizes all the proposed structures for hydrocarbons and other compounds detected in TF samples.
JH quantification by GC-MS
Request a detailed protocolFor each sample, a known quantity of TF or hemolymph was collected into a graduated glass capillary tube and blown into an individual glass vial containing 5 µL of 100% ethanol. Samples were kept at −20°C until further processing. This biological sample was added to a 1:1 mixture of isooctane and methanol, vortexed for 30 s, and centrifuged for 7 min at maximum speed. Avoiding the boundary between phases, the majority of the isooctane layer and the methanol layer were removed separately, combined and stored at –80°C until analysis. Before analysis, 50% acetonitrile (HPLC grade) was added. Prior to purification, farnesol (Sigma-Aldrich, St Louis, MO) was added to each sample to serve as an internal standard. Samples were extracted three times with hexane (HPLC grade). The hexane fractions were recombined in a clean borosilicate glass vial and dried by vacuum centrifugation. JH III was quantified using the gas chromatography mass spectrometry (GC–MS) method of (Bergot et al., 1981) as modified by (Shu et al., 1997) and (Brent and Vargo, 2003). The residue was washed out of the vials with three rinses of hexane and added to borosilicate glass columns filled with aluminum oxide. In order to filter out contaminants, samples were eluted through the columns successively with hexane, 10% ethyl ether–hexane and 30% ethyl ether–hexane. After drying, samples were derivatized by heating at 60°C for 20 min in a solution of methyl-d alcohol (Sigma-Aldrich) and trifluoroacetic acid (Sigma-Aldrich, St Louis, MO). Samples were dried down, resuspended in hexane, and again eluted through aluminum oxide columns. Non-derivatized components were removed with 30% ethyl ether. The JH derivative was collected into new vials by addition of 50% ethyl-acetate–hexane. After drying, the sample was resuspended in hexane. Samples were then analyzed using an HP 7890A Series GC (Agilent Technologies, Santa Clara, CA) equipped with a 30 m x 0.25 mm Zebron ZB-WAX column (Phenomenex, Torrence, CA) coupled to an HP 5975C inert mass selective detector. Helium was used as a carrier gas. JH form was confirmed by first running test samples in SCAN mode for known signatures of JH 0, JH I, JH II, JH III and JH III ethyl; JH III was confirmed as the primary endogenous form in this species. Subsequent samples were analyzed using the MS SIM mode, monitoring at m/z 76 and 225 to ensure specificity for the d3-methoxyhydrin derivative of JH III. Total abundance was quantified against a standard curve of derivatized JH III, and adjusted for the starting volume of TF. The detection limit of the assay is approximately one pg.
Long-term development
Request a detailed protocolTo determine the effect of exogenous JH on larval development, ants were taken from laboratory C. floridanus colonies (Expt 1: C2, C3, C5, C9, C16, C17; Expt 2: C4, C5, C6, C18; Expt 3: C1, C5, C6, C7, C11, C16, C19). Approximately 90% of the ants were taken from inside the nest on the brood, while the remaining 10% were taken from outside the nest. Each colony explant had 20–30 workers (each treatment had the same number of replicates of any given colony) and was provided with five to ten second or third instar larvae from their own colony of origin (staged larvae were equally distributed across replicates). Each explant was provided with water, and either solvent- or JH III-supplemented 30% sugar water and maple-syrup-based ant diet (1500 ng of JH III in 5 µL of ethanol was applied to each 3×3×3 mm food cube and sucrose solution had 83 ng JH III/µL). No insect-based food was provided. Food sources were refreshed twice per week. JH was found to transition to JH acid gradually at room temperature, where after 1 week ~ 50% was JH acid (as measured by radio-assay as in (Kamita et al., 2011), data not shown).
Twice weekly before feeding, each explant was checked for pupae, and developing larvae were measured and counts using a micrometer in the reticle of a stereomicroscope. Upon pupation, or cocoon spinning, larvae/pupae were removed from the care of workers and kept in a clean humid chamber until metamorphosis. Cocoons were removed using dissection forceps. The head width of the pupae was measured using a micrometer in the reticle of a stereomicroscope 1–4 days after metamorphosis (immediately after removal of the larval sheath, head width is not stable). Long-term development experiments were stopped when fewer than three larvae remained across all explants and these larvae had not changed in size for 2 weeks. Of larvae that did not successfully undergo metamorphosis, approximately 75% were eaten by nursing workers at varying developmental stages over the course of the experiment. The remaining non-surviving larvae were split between larvae that finished the cocoon spinning phase (Wallis, 1960) but did not complete metamorphosis and larvae that had ceased to grow by the end of the approximately 10-week experiment.
C. fellah transcriptomics
Transcriptome sequencing
Request a detailed protocolWe sequenced the C. fellah transcriptome (RNAseq) of workers from a single colony initiated from a queen collected during a mating flight in 2007 in Tel Aviv, Israel. Total RNA was extracted from the whole body of four minor workers with RNeasy Plus micro kit and RNase-Free DNase Set (QIAGEN, Hilden, Germany) and then 350 ng cDNA from each ant were pooled and sequenced. Illumina cDNA library was constructed and sequenced using strand-specific, paired-end sequencing of 100 bp reads. The library was sequenced on an Illumina HiSeq machine, which generated a total of 115 million pairs of reads.
Transcriptome assembly
Request a detailed protocolWe ran Trinity (Grabherr et al., 2011) (version r2013-02–25, RRID:SCR_013048) on these sequence reads to assemble the C. fellah transcriptome. Reads were filtered according to the Illumina Chastity filter and then trimmed and filtered using Trimmomatic (Bolger et al., 2014) (RRID:SCR_011848, parameters: ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 SLIDINGWINDOW:4:15 MINLEN:50). Trinity was run with default parameters for strand-specific, paired-read data. The assembled transcriptome consists of 66,156 genes (‘components’) with 99,402 transcripts, with putative open-reading frames found in 9987 and 27,062 of them, respectively. The total sequence length is 109,526,661 bases (including alternative splice variants) and the N50 size is 2243 bases. The C. fellah transcriptome dataset is available under accession number PRJNA339034.
Protein orthology
Request a detailed protocolIn order to identify orthologous proteins across the species TF, we first needed to determine orthology across the four species. Compiling RefSeq, UniProt, and transcriptome protein models from the four species yielded 131,122 predicted protein sequences. For C. floridanus, A. mellifera and S. invicta, this was done for both NCBI RefSeq and UniProt databases because there are discrepancies in annotation and thus in protein identification between databases. We determined 21,836 groups of one-to-one orthologs using OMA stand-alone (Altenhoff et al., 2015) v.1.0.3, RRID:SCR_011978, although only 4538 orthologous groups had members in all four species. Default parameters were used with the exception of minimum sequence length, which was lowered to 30 aa.
The 40 proteins with the highest average NSAF value across TF samples of that species were selected. For each of the top proteins, we checked across the other species and databases for orthologous proteins. If an orthologous protein was identified in the proteome, we then checked if that protein was also present in that species’ TF, even at low abundance.
Dendrogram
Request a detailed protocolThe protein sequences of our proteomically identified TF C. floridanus JHE/Est-6 proteins, functionally validated JHEs (Tribolium castaneum JHE (UniProt D7US45), A. mellifera JHE (Q76LA5), Manduca sexta JHE (Q9GPG0), D. melanogaster JHE (A1ZA98) and Culex quinquefasciatus (R4HZP1)) and Est-6 proteins from A. mellifera (B2D0J5) and D. melanogaster (P08171) were aligned using PROMALS3D with the crystal structure of the M. sexta JHE (2fj0) used as a guide. Phylogeny was inferred using Randomized Axelerated Maximum Likelihood (RAxML, RRID:SCR_006086 [Stamatakis, 2006]), with 100 bootstrapped trees. The dendrogram was visualized in FigTree (v1.4.2, RRID:SCR_008515).
Sample sizes, data visualization and statistics
Request a detailed protocolThe proteomic sample sizes were determined by the variation observed in protein IDs and abundances in unmanipulated samples. Visualization together with hierarchical clustering was done in R version 3.0.2 (www.r-project.org, RRID:SCR_001905) using the ‘heatmap.3’ package in combination with the pvclust package. Heatmap visualization without clustering was done in MATLAB (2012b, RRID:SCR_001622). For experiments in Figures 3 and 4, the number of colonies and number of replicates or ants per colony were determined by both preliminary trials to assess sample variation and the health/abundance of the C. floridanus ant colonies available in the lab. Hydrocarbon GC-MS traces were analyzed in MATLAB normalizing the abundance (area under the curve) of each point by maximum and minimum value within the CHC retention time window. Peaks were found using ‘findpeaks’ with a lower threshold abundance of 7% of the total abundance for that sample. To compare across samples, peaks were filtering into 0.03 min bins by retention time. Cross-correlation was computed using the ‘xcorr’ function. For long-term development experiments, the number of same-staged larvae per colony was the limiting factor for the number of replicates per experiment. Because of this limitation, the long-term development experiment testing the effect of JH was fully repeated three times. Colonies were hibernated approximately 1 month prior to each of these experiments to maximize number of same-staged brood. General linear mixed models (GLMM) were used so that colony, replicate and experiment could be considered as random factors. Models were done in R using ‘lmer’ and ‘glmer’ functions of the lme4 package, and p values were calculated with the ‘lmerTest’ package in R. Overall, no data points were excluded as outliers and all replicates discussed are biological not technical replicates.
Data availability
-
Social exchange of chemical cues, growth proteins and hormones through trophallaxisPublicly available at ProteomeXchange (accession no. PXD004825).
-
Camponotus fellah Transcriptome or Gene expressionPublicly available at the NCBI BioProject database (accession no: PRJNA339034).
-
Camponotus floridanus Transcriptome or Gene expressionPublicly available at the NCBI BioProject database (accession no: PRJNA338939).
References
-
CmCatD, a cathepsin D-like protease has a potential role in insect defense against a phytocystatinJournal of Insect Physiology 55:678–685.https://doi.org/10.1016/j.jinsphys.2009.04.016
-
BookAn Introduction to Systems Biology: Design Principles of Biological CircuitsBoca Raton, Florida: Chapman and Hall/CRC .
-
Human milk as a carrier of biochemical messagesActa Paediatrica 88:27–41.https://doi.org/10.1111/j.1651-2227.1999.tb01298.x
-
Trimmomatic: a flexible trimmer for Illumina sequence dataBioinformatics 30:2114–2120.https://doi.org/10.1093/bioinformatics/btu170
-
Self-organization in social insectsTrends in Ecology & Evolution 12:188–193.https://doi.org/10.1016/S0169-5347(97)01048-3
-
Octopamine reverses the isolation-induced increase in trophallaxis in the carpenter ant Camponotus fellahThe Journal of Experimental Biology 203:513–520.
-
Neutrophil extracellular traps kill bacteriaScience 303:1532–1535.https://doi.org/10.1126/science.1092385
-
Purification and kinetic characterisation of juvenile hormone esterase from Drosophila melanogasterInsect Biochemistry and Molecular Biology 28:501–515.https://doi.org/10.1016/S0965-1748(98)00037-X
-
A duration constant for worker-to-larva trophallaxis in fire antsInsectes Sociaux 43:149–166.https://doi.org/10.1007/BF01242567
-
The critical period for caste determination in Bombus terrestris and its juvenile hormone correlatesJournal of Comparative Physiology A: Sensory, Neural, and Behavioral Physiology 186:1089–1094.https://doi.org/10.1007/s003590000163
-
Esterase-6 and the pheromonal effects of cis-vaccenyl acetate in Drosophila melanogasterJournal of Evolutionary Biology 2:395–407.https://doi.org/10.1046/j.1420-9101.1989.2060395.x
-
Analysis of two genetic models for the innate components of colony odor in social HymenopteraBehavioral Ecology and Sociobiology 4:217–224.https://doi.org/10.1007/BF00297645
-
Trophallaxis mediates uniformity of colony odor in Cataglyphis iberica ants (Hymenoptera, Formicidae)Journal of Insect Behavior 12:559–567.https://doi.org/10.1023/A:1020975009450
-
Worker division of labor and endocrine physiology are associated in the harvester ant, Pogonomyrmex californicusJournal of Experimental Biology 215:454–460.https://doi.org/10.1242/jeb.060822
-
ConferenceThe evolution and social significance of the ant proventriculusProceedings Tenth International Congress of Entomology. pp. 503–508.
-
The interactions between juvenile hormone (JH), lipophorin, vitellogenin, and JH esterases in two cockroach speciesInsect Biochemistry and Molecular Biology 30:793–803.https://doi.org/10.1016/s0965-1748(00)00051-5
-
Honeybees learn floral odors while receiving nectar from foragers within the hiveDie Naturwissenschaften 94:55–60.https://doi.org/10.1007/s00114-006-0157-3
-
Hormonal pleiotropy and the juvenile hormone regulation of Drosophila development and life historyBioEssays: News and Reviews in Molecular, Cellular and Developmental Biology 27:999–1010.https://doi.org/10.1002/bies.20290
-
BookInsect Development: Morphogenesis, Molting and MetamorphosisGilbert L. I, editors. Elsevier Science.
-
Full-length transcriptome assembly from RNA-Seq data without a reference genomeNature Biotechnology 29:644–652.https://doi.org/10.1038/nbt.1883
-
Behavioral plasticity in honey bees is associated with differences in brain microRNA transcriptomeGenes, Brain, and Behavior 11:660–670.https://doi.org/10.1111/j.1601-183X.2012.00782.x
-
Propagation of olfactory information within the honeybee hiveBehavioral Ecology and Sociobiology 60:707–715.https://doi.org/10.1007/s00265-006-0214-0
-
The evolution of hydrocarbon pheromone parsimony in ants (Hymenoptera : Formicidae) – interplay of colony odor uniformity and odor idiosyncrasy. A reviewMyrmecological News 10:59–68.
-
A review of saliva: Normal composition, flow, and functionThe Journal of Prosthetic Dentistry 85:162–169.https://doi.org/10.1067/mpr.2001.113778
-
BookNourishment and the evolution of the social VespidaeIn: Ross K. G, Matthews R. W, editors. The Social Biology of Wasps. Ithaca, NewYork: Cornell University Press. pp. 426–450.
-
The juvenile hormone signaling pathway in insect developmentAnnual Review of Entomology 58:181–204.https://doi.org/10.1146/annurev-ento-120811-153700
-
Shrimp miRNAs regulate innate immune response against white spot syndrome virus infectionDevelopmental & Comparative Immunology 60:191–201.https://doi.org/10.1016/j.dci.2016.03.002
-
Juvenile hormone esterase: biochemistry and structureJournal of Pesticide Science 35:265–274.https://doi.org/10.1584/jpestics.R10-09
-
Molecular characterization and developmental expression of a retinoid- and fatty acid-binding glycoprotein from Drosophila. A putative lipophorinThe Journal of Biological Chemistry 271:20641–20649.https://doi.org/10.1074/jbc.271.34.20641
-
Direct behavioral evidence for hydrocarbons as Ant recognition discriminatorsNaturwissenschaften 86:246–249.https://doi.org/10.1007/s001140050609
-
A bacterial filter protects and structures the gut microbiome of an insectThe ISME Journal 10:1866–1876.https://doi.org/10.1038/ismej.2015.264
-
Primer effects of a brood pheromone on honeybee behavioural developmentProceedings. Biological Sciences 268:163–168.https://doi.org/10.1098/rspb.2000.1345
-
Me and we: the interplay between individual and group behavioral variation in social collectivesCurrent Opinion in Insect Science 5:16–24.https://doi.org/10.1016/j.cois.2014.09.010
-
The host defense of Drosophila melanogasterAnnual Review of Immunology 25:697–743.https://doi.org/10.1146/annurev.immunol.25.022106.141615
-
Identification and characterization of a juvenile hormone response element and its binding proteinsJournal of Biological Chemistry 282:37605–37617.https://doi.org/10.1074/jbc.M704595200
-
Larval and nurse worker control of developmental plasticity and the evolution of honey bee queen-worker dimorphismJournal of Evolutionary Biology 24:1939–1948.https://doi.org/10.1111/j.1420-9101.2011.02331.x
-
WebsiteNIST Chemistry WebBook, NIST Standard Reference1258 Database Number 69NIST Chemistry WebBook.
-
The serpin saga; development of a new class of virus derived anti-inflammatory protein immunotherapeuticsAdvances in Experimental Medicine and Biology 666:132–156.https://doi.org/10.1007/978-1-4419-1601-3_11
-
Natural selection and the communal rearing of pups in African wild dogs (Lycaon pictus)Behavioral Ecology and Sociobiology 10:1–13.https://doi.org/10.1007/BF00296390
-
Effect of a brood pheromone on honeybee hypopharyngeal glandsComptes Rendus De l'Academie Des Sciences. Serie III, Sciences De La Vie 319:769–772.
-
IRS and TOR nutrient-signaling pathways act via juvenile hormone to influence honey bee caste fateThe Journal of Experimental Biology 214:3977–3984.https://doi.org/10.1242/jeb.061499
-
A statistical model for identifying proteins by tandem mass spectrometryAnalytical Chemistry 75:4646–4658.https://doi.org/10.1021/ac0341261
-
The developmental control of size in insectsWiley Interdisciplinary Reviews. Developmental Biology 3:113–134.https://doi.org/10.1002/wdev.124
-
Juvenile hormone and the physiological basis of insect polymorphismsThe Quarterly Review of Biology 57:109–133.https://doi.org/10.1086/412671
-
The seminal symphony: how to compose an ejaculateTrends in Ecology & Evolution 28:414–422.https://doi.org/10.1016/j.tree.2013.03.005
-
Complexity of seminal fluid: a reviewBehavioral Ecology and Sociobiology 60:289–310.https://doi.org/10.1007/s00265-006-0178-0
-
How did Arthropod Sesquiterpenoids and Ecdysteroids Arise? Comparison of hormonal pathway genes in noninsect arthropod genomesGenome Biology and Evolution 7:1951–1959.https://doi.org/10.1093/gbe/evv120
-
Extracellular communication via microRNA: lipid particles have a new messageJournal of Lipid Research 54:1174–1181.https://doi.org/10.1194/jlr.R034991
-
BookInsect Saliva: Function, Biochemistry, and PhysiologyIn: Chapman R. F, Boer G, editors. Regulatory Mechanisms in Insect Feeding. Chapman & Hall. pp. 74–97.
-
Cathepsins of lepidopteran insects: Aspects and prospectsInsect Biochemistry and Molecular Biology 64:51–59.https://doi.org/10.1016/j.ibmb.2015.07.005
-
Tertiary siRNAs mediate paramutation in C. elegansPLoS Genetics 11:e1005078.https://doi.org/10.1371/journal.pgen.1005078
-
Hydrocarbon site of synthesis and circulation in the desert ant Cataglyphis nigerJournal of Insect Physiology 46:1097–1102.https://doi.org/10.1016/s0022-1910(99)00219-x
-
Hydrocarbon dynamics within and between nestmates inCataglyphis niger (Hymenoptera: Formicidae)Journal of Chemical Ecology 21:365–378.https://doi.org/10.1007/BF02036724
-
Role of the prophenoloxidase-activating system in invertebrate immunityCurrent Opinion in Immunology 10:23–28.https://doi.org/10.1016/S0952-7915(98)80026-5
-
Function and evolution of microRNAs in eusocial HymenopteraFrontiers in Genetics 6:193.https://doi.org/10.3389/fgene.2015.00193
-
Evolution of digestive systems of insectsAnnual Review of Entomology 35:181–200.https://doi.org/10.1146/annurev.en.35.010190.001145
-
Extracellular miRNA: A collision of two paradigmsTrends in Biochemical Sciences 41:883–892.https://doi.org/10.1016/j.tibs.2016.08.004
-
Extracellular miRNAs: the mystery of their origin and functionTrends in Biochemical Sciences 37:460–465.https://doi.org/10.1016/j.tibs.2012.08.003
-
Insect Hydrocarbons: Biology, Biochemistry and Chemical Ecology222–243, Nestmate recognition in social insects and the role of hydrocarbons, Insect Hydrocarbons: Biology, Biochemistry and Chemical Ecology, Cambridge University Press, 10.1017/CBO9780511711909.012.
-
Early developmental processes limit socially mediated phenotypic plasticity in an antBehavioral Ecology and Sociobiology 70:285–291.https://doi.org/10.1007/s00265-015-2052-4
-
Spinning movements in the larvae of the ant, Formica fuscaInsectes Sociaux 7:187–199.https://doi.org/10.1007/BF02224080
-
Hierarchical chemosensory regulation of male-male social interactions in DrosophilaNature Neuroscience 14:757–762.https://doi.org/10.1038/nn.2800
-
Analysis of the catalytic mechanism of juvenile hormone esterase by site-directed mutagenesisInternational Journal of Biochemistry 24:1933–1941.https://doi.org/10.1016/0020-711X(92)90289-D
-
The microRNA spectrum in 12 body fluidsClinical Chemistry 56:1733–1741.https://doi.org/10.1373/clinchem.2010.147405
-
Developmental and physiological determinants of caste in social Hymenoptera: Evolutionary implicationsThe American Naturalist 128:13–34.https://doi.org/10.1086/284536
-
A study of some ant larvae, with a consideration of the origin and meaning of the social habit among insectsProceedings of the American Philosophical Society 57:293–339.
-
The function of the corpus allatum in the growth and reproduction of Rhodnius Prolixus (Hemiptera)Journal of Cell Science 79:91–121.
-
Quantitative studies of liquid food transmission in antsInsectes Sociaux 4:157–166.https://doi.org/10.1007/BF02224149
-
Identification of candidate odorant degrading gene/enzyme systems in the antennal transcriptome of Drosophila melanogasterInsect Biochemistry and Molecular Biology 53:30–43.https://doi.org/10.1016/j.ibmb.2014.07.003
-
Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiaeJournal of Proteome Research 5:2339–2347.https://doi.org/10.1021/pr060161n
Article and author information
Author details
Funding
European Research Council (Advanced Grant 249375)
- Laurent Keller
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung
- Richard Benton
- Laurent Keller
European Research Council (Starting Independent Researcher 205202)
- Richard Benton
European Research Council (Consolidator Grant 615094)
- Richard Benton
Wellcome (Wellcome Trust grant 104640/Z/14/Z)
- Eric A Miska
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior
- Zamira G Soares
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Robert I’Anson Price and Thomas Richardson for access to A. mellifera colonies, and Raphaël Braunschweig for access to C. fellah colonies. We thank George Shizuo Kamita for advice on JHEs. We would like to thank Sylviane Moss and Kay Harnish for help with high-throughput sequencing and Alexandra Bezler, Thomas Auer, Roman Arguello, Juan Alcaniz-Sanchez, Pavan Ramdya, Lucia Prieto-Godino, Katja Hoedjes, Raphaël Braunschweig and Sean McGregor for their useful comments on the manuscript. This work was funded by an ERC Advanced Grant (249375) and a Swiss NSF grant to LK, an ERC Starting Independent Researcher and Consolidator Grants (205202 and 615094) and a Swiss NSF Grant to R.B. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. USDA is an equal opportunity provider and employer. ZSG was supported by CAPES Brazil. EAM was supported by Wellcome Trust (104640/Z/14/Z).
Copyright
© 2016, LeBoeuf et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 7,361
- views
-
- 1,385
- downloads
-
- 103
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Chromosomes and Gene Expression
- Evolutionary Biology
Gene regulation is essential for life and controlled by regulatory DNA. Mutations can modify the activity of regulatory DNA, and also create new regulatory DNA, a process called regulatory emergence. Non-regulatory and regulatory DNA contain motifs to which transcription factors may bind. In prokaryotes, gene expression requires a stretch of DNA called a promoter, which contains two motifs called –10 and –35 boxes. However, these motifs may occur in both promoters and non-promoter DNA in multiple copies. They have been implicated in some studies to improve promoter activity, and in others to repress it. Here, we ask whether the presence of such motifs in different genetic sequences influences promoter evolution and emergence. To understand whether and how promoter motifs influence promoter emergence and evolution, we start from 50 ‘promoter islands’, DNA sequences enriched with –10 and –35 boxes. We mutagenize these starting ‘parent’ sequences, and measure gene expression driven by 240,000 of the resulting mutants. We find that the probability that mutations create an active promoter varies more than 200-fold, and is not correlated with the number of promoter motifs. For parent sequences without promoter activity, mutations created over 1500 new –10 and –35 boxes at unique positions in the library, but only ~0.3% of these resulted in de-novo promoter activity. Only ~13% of all –10 and –35 boxes contribute to de-novo promoter activity. For parent sequences with promoter activity, mutations created new –10 and –35 boxes in 11 specific positions that partially overlap with preexisting ones to modulate expression. We also find that –10 and –35 boxes do not repress promoter activity. Overall, our work demonstrates how promoter motifs influence promoter emergence and evolution. It has implications for predicting and understanding regulatory evolution, de novo genes, and phenotypic evolution.