Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

Suppression subtractive hybridization analysis reveals expression of conserved and novel genes in male accessory glands of the ant Leptothorax gredleri

Angelika Oppelt1*, Fernanda C Humann2, Marion Fuessl1, Sergio V Azevedo3, David S Marco Antonio3, Jürgen Heinze1 and Klaus Hartfelder2

Author Affiliations

1 Biologie I, Universität Regensburg, Universitätsstr. 31, D-93040 Regensburg, Germany

2 Departamento de Biologia Celular e Molecular e Bioagentes Patogênicos, Faculdade de Medicina de Ribeirão Preto, Universidade de São Paulo, Av. Bandeirantes 3900, 14049-900, Ribeirão Preto, Brazil

3 Departamento de Genética, Faculdade de Medicina de Ribeirão Preto, Universidade de São Paulo, Av. Bandeirantes 3900, 14049-900, Ribeirão Preto, Brazil

For all author emails, please log on.

BMC Evolutionary Biology 2010, 10:273  doi:10.1186/1471-2148-10-273

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/10/273


Received:4 May 2010
Accepted:8 September 2010
Published:8 September 2010

© 2010 Oppelt et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

During mating, insect males eject accessory gland proteins (Acps) into the female genital tract. These substances are known to affect female post-mating behavior and physiology. In addition, they may harm the female, e.g., in reducing its lifespan. This is interpreted as a consequence of sexual antagonistic co-evolution. Whereas sexual conflict abounds in non-social species, the peculiar life history of social insects (ants, bees, wasps) with lifelong pair-bonding and no re-mating aligns the reproductive interests of the sexes. Harming the female during mating would negatively affect male fitness and sexual antagonism is therefore not expected. Indeed, mating appears to increase female longevity in at least one ant species. Acps are presumed to play a role in this phenomenon, but the underlying mechanisms are unknown. In this study, we investigated genes, which are preferentially expressed in male accessory glands of the ant Leptothorax gredleri, to determine which proteins might be transferred in the seminal fluid.

Results

By a suppression subtractive hybridization protocol we obtained 20 unique sequences (USs). Twelve had mutual best matches with genes predicted for Apis mellifera and Nasonia vitripennis. Functional information (Gene Ontology) was available only for seven of these, including intracellular signaling, energy-dependent transport and metabolic enzyme activities. The remaining eight USs did not match sequences from other species. Six genes were further analyzed by quantitative RT-PCR in three life cycle stages of male ants. A gene with carboxy-lyase activity and one of unpredicted function were significantly overexpressed in accessory glands of sexually mature males.

Conclusions

Our study is the first one to investigate differential gene expression in ants in a context related to mating. Our findings indicate that male accessory glands of L. gredleri express a series of genes that are unique to this species, possibly representing novel genes, in addition to conserved ones for which functions can be predicted. Identifying differentially expressed genes might help to better understand molecular mechanisms involved in reproductive processes in eusocial Hymenoptera. While the novel genes could account for rapidly evolving ones driven by intra-sexual conflict between males, conserved genes imply that rather beneficial traits might get fixed by a process described as inter-sexual cooperation between males and females.

Background

Substances produced by the male accessory glands (MAGs) of insects and transferred into the female genital tract during mating are known to reduce pathogen transmission, to form mating plugs or spermatophores, and to be important in sperm competition. In addition, they trigger fundamental changes in female physiology, behavior, and reproduction [1-4]. The major biologically active components of MAG secretions are carbohydrates, lipids, and in particular accessory gland proteins (Acps) [1].

In Drosophila melanogaster, more than 100 Acps have been identified [5]. They play roles in the regulation of female receptivity [6-9], defense against bacteria [10], the activation of the female immune system [11,12], the efficiency of sperm storage [13,14], and the stimulation of ovulation and egg maturation [6,8,9]. They are also decisive for the fertilization of eggs, as shown in experiments in which the activity of sperm of infertile mutant males lacking MAGs could be restored by subsequently re-mating females with males that provided Acps but no sperm [15].

In addition to stimulating immediate female fertility, MAG products appear to negatively affect female lifespan [16,17]. This reflects inter-sexual conflict [18]: the reproductive interests of males and females concerning the consequences of a mating act may differ, in particular when females later re-mate with other partners. Males prefer the female to maximize short-term investment into the offspring resulting from the present copulation, while females may benefit more from saving resources for later reproduction with other males. Sexual conflict may result in sexually antagonistic co-evolution and Acps indeed appear to evolve very rapidly, as predicted from the ensuing arms race [5,19-23].

Compared to the wealth of information on accessory gland products of drosophilids, studies on other insect species are seriously lagging behind [1]. This is particularly the case for eusocial Hymenoptera, despite the special implications that their life style and evolutionary history have on sexual conflict. In contrast to Drosophila and other non-social insects, eusocial Hymenoptera are characterized by life-long pair-bonding, even though the males frequently die after mating [24-26]. Investigations have shown that monogamy is the ancestral state of all truly eusocial organisms [27], and lifetime monogamy seems to be a universal precondition for the evolution of obligate eusociality [24,25]. The queens of ants, eusocial bees, and wasps indeed mate only once or a few times early in their adult lives without ever re-mating again. Sperm received during these matings is stored for months, years, or even decades, and utilized to fertilize eggs throughout the queens' lives. Furthermore, queens first rear large numbers of sterile workers from fertilized eggs before switching to the production of female sexuals. Hence, in contrast to males of non-social insects, eusocial insect males are expected to avoid harming the female, as this would reduce the total amount of female sexuals produced and, consequently, male reproductive success. Instead, males should have evolved mechanisms that promote female performance. Even where multiple mating has evolved after the monogamy window has been passed, competition among ejaculates [28] is not expected to negatively affect female life span and future fecundity. Inter-sexual conflict might thus be replaced by "inter-sexual cooperation", though intra-sexual conflict might still evolve. Indeed, ant males appear to positively affect queen life span: regardless of being mated with a fertile or sterilized male, mated queens of Cardiocondyla obscurior lived considerably longer than virgin queens [29].

Like in other insects, MAG products of eusocial Hymenoptera are involved in the formation of mating plugs [30-32] and ejaculate (or eventually sperm) competition [28,33]. In addition, they presumably elicit rapid changes in the physiology and behavior of freshly mated queens, e.g., the shedding of wings, reluctance of mating again, and changes in cuticular hydrocarbons [34]. It is likely that MAG products, and in particular Acps, are also involved in the life-prolonging effects of mating [29]. Nevertheless, very little is known about Acps in eusocial Hymenoptera other than the honey bee, Apis mellifera [35-38]. We therefore examined, which genes are differentially expressed in the male accessory glands of the ant Leptothorax gredleri.

Virgin queens of L. gredleri attract males through "female calling", i.e., they climb up grass stems or similar, elevated positions near their maternal nests and release droplets of a sex pheromone from their poison glands [39-41]. During mating, sperm is transferred within a gelatinous substance, functionally similar to a spermatophore. The latter stays inside the female's genital tract for several hours, while the sperm is slowly transferred into the spermatheca [42]. Females have been observed to copulate with multiple males but it is unclear whether ejaculate is transferred during all these matings. Genetic analyses document single paternity of their offspring [40].

In order to obtain information on gene products specifically produced by the accessory glands of L. gredleri males we generated a suppression subtractive library for transcripts using a Representational Difference Analysis (RDA) protocol. Such an approach enriches tissue-specific transcripts (here: MAGs) by removing from the library sequences that are shared with control tissue (here: the body of the males after excision of MAGs).

Results

RDA library characteristics and bioinformatic analyses

To detect differentially expressed genes in L. gredleri MAGs, we sequenced 192 clones from the subtracted gland library, resulting in 288 sequencing reads, of which 158 reads representing 130 different clones were of appropriate quality (Phred quality ≥ 20). After quality analysis and following CAP3 assembly we obtained 12 contigs and 8 singlets as Unique Sequences (USs). The sequences were submitted to GenBank (accession numbers GT897811 until GT897966).

BLASTx analyses of the USs (Table 1) revealed significant matches (E-scores < e-5, a cut-off level frequently used in EST projects; [43]) for 12 USs. The remaining 8 USs showed no significant matches to genes of other species.

Table 1. EST sequences from RDA libraries of accessory glands cDNA from Leptothorax gredleri

Most of the ESTs (40) clustered to a sequence corresponding to an Apis mellifera predicted gene similar to CG1486-PA. Based on domain predictions, the fly ortholog of this gene is annotated in Flybase as being involved in carboxylic acid metabolism. Another US represented by 12 ESTs showed similarity to a heat shock factor binding protein 1-like from Apis mellifera, with the corresponding CG5446-PA in Drosophila melanogaster. Further USs were also formed by clustering considerable numbers of ESTs (Table 1), such as US4 consisting of 28 ESTs, and US3 composed of 10 ESTs. None of these transcripts, however, presented matches when submitted to BLASTx analysis.

In functional terms, the searches for orthologs revealed that, seven of the L. gredleri USs corresponded to functionally defined genes (Table 1), with best hits to predicted A. mellifera proteins. By judging from Gene Ontology attributes, two are expected to be involved in energy-dependent transport processes, such as the ATP synthase gamma chain and Nervana 2, a sodium/potassium-transporting ATPase subunit. The Leukocyte-antigen-related-like protein may function as a protein tyrosine phosphatase, while the kokopelli gene has a cyclin-dependent protein kinase regulator activity and is involved in male germ-line stem cell division. The final two sequences with predicted molecular functions are CG6910-PA, with inositol oxygenase activity, and an Apis mellifera ribosomal protein LP1 (RpLP1) given as being a structural constituent of ribosomes.

Still unknown are biological processes or molecular functions for the five remaining sequences with BLASTx matches. SJCHGC05576 and LOC100123166 are both predicted proteins of the parasitic wasp Nasonia vitripennis, and Ciao-1 and CG32432-PA and heat shock factor binding protein 1-like are predicted proteins for A. mellifera. Neither of these had matches with Drosophila genes.

Quantitative RT-PCR analysis of RDA library genes

The RDA method is applicable to systems where genomic information is either completely lacking or still in an incipient stage. RDAs are able to detect both known and so far unknown transcripts, including regulatory RNAs. However, being a high throughput approach, transcripts indicated by RDA as differentially expressed should be further analyzed by specific gene-expression analyses, such as quantitative RT-PCR (qRT-PCR). Thus, we designed gene-specific primers for 19 of the 20 RDA products. One product contained repeats, and so it was not considered suitable for primer design. For six of the 19 products we obtained primers, shown in Table 2, that reliably amplified specific products and, thus, were suitable for qRT-PCR. These were primers for the genes leukocyte-antigen-related-like (leukocyte) and CG1486, representing two genes with BLASTx matches and predicted functions, for ciao-1 and CG32432, being genes with matches but of unknown function, and the final two being for US3 and US7, representing no-match sequences.

Table 2. qRT-PCR primers sequences with their target genes (when known)

For five of the six genes we could show that their mean transcript levels are at least two-times higher in the L. gredleri MAGs relative to body carcasses, this holding for all three of the analyzed male life cycle stages, except for ciao-1 in the males running outside the nest that were sexually mature (males outside nests) samples (Figure 1). For the sixth gene, US7, overexpression was not validated. This may, thus, represent false positive RDA identification, represented by a single sequencing read (Table 1).

thumbnailFigure 1. Gene expression in male accessory glands. Relative expression of leukocyte-antigen-related-like, CG1486, ciao-1, CG32432, US3 and US7 transcripts in male accessory glands of the ant Leptothorax gredleri relative to body carcasses. The y-axis, drawn through 1, corresponds to the expression level in the body carcass (control tissue) of the respective life cycle stage, these being late male pupae, males inside nests and males outside nests. Shown are fold change expression ± SEM for five biological replicates. Statistically significant differences (REST randomization significance test, P < 0.05) are indicated by asterisks.

When testing for statistical significance of differential expression by randomization tests built into REST software we could show significant overexpression in two genes. CG1486, a gene with predicted carboxy-lyase activity, was not only represented by the highest number of ESTs (Table 1), but was also confirmed as significantly overexpressed in MAGs of dark-pigmented male pupae (male pupae; P = 0.010) and sexually mature males (males outside nests; P = 0.015). The other gene, CG32432, was shown to be significantly overexpressed in sexually mature males (P = 0.037) (Figure 1). For all of these five genes biological variation in expression levels was considerable, a problem frequently encountered in qRT-PCR analyses, even with highly controlled cell cultures [44].

Discussion and conclusions

Using the RDA method we successfully identified a set of 20 genes as being differentially expressed in MAGs of the ant L. gredleri. Unfortunately, a fully sequenced ant genome is not yet available, so an annotation against a genome was not feasible. Such an annotation would be needed to identify proteins that are secreted into the lumen of L. gredleri MAGs, as their RNAs should have hydrophobic amino-terminal signal sequences [45]. For example, Swanson et al. [22] found that 24% of the 212 identified genes in their study of Drosophila male reproductive proteins have such a putative signal sequence.

All best hit matches were against A. mellifera or N. vitripennis predicted genes, emphasizing the importance of these hymenopteran reference genomes in future gene expression studies in eusocial insects. Nevertheless, most of our ESTs represented functionally undefined genes or ones without matches to known genes. At first sight, this might indicate rapid evolution and species specificity in MAG secretion products, similar to what has been shown for MAGs of Drosophila [22], where about 50% of the Drosophila Acps have no homologs outside the drosophilid species complex. When comparing the genes predicted to be differentially expressed in MAGs of L. gredleri with those identified in a proteome analysis of honey bee seminal fluid [38] we did not find any matching gene products for the set with GeneOntology functional attributes. This finding would also be in favor of a rapid evolution hypothesis for MAG products. Queens of Leptothorax are predominantly singly-mated [40,46], but occasional multiple mating might introduce intra-sexual conflict and thus increase the evolutionary tempo of seminal fluid proteins. In eusocial species with multiple mating, the beneficial effect of seminal fluid on sperm survival is reduced when seminal fluid is mixed with alien sperm [28]. However, this may actually serve the interests of the queen rather than inflicting damage. Furthermore, the queen appears to be in control of ejaculate competition, since the negative impact of seminal fluid on foreign sperm can be neutralized by secretions from the queen sperm-storage organ [28]. In this case, intra-sexual conflict seems to be controlled by inter-sexual cooperation.

Concerning our studies, nothing is as yet known on the molecular nature of MAG secretions of other ant species, and with Apis mellifera and Nasonia vitripennis being the closest species for genomic comparisons, the percentage of functionally undefined genes and of no-match ESTs might simply reflect the large phylogenetic distance separating ants, bees, and parasitic wasps. Speculating about the speed of the evolution of L. gredleri Acps would, at this point, be premature, especially since, 2D-electrophoreses of Acps suggest that they are conserved in L. gredleri and related taxa relative to, e.g., muscle proteins (M. Fuessl, A. Oppelt, J. Heinze, unpublished).

Furthermore, not all seminal fluid genes can be subject to rapid evolution, because successful reproduction must be assured. In a cross-species comparison of Drosophila MAG protein encoding genes, Mueller et al. [23] found not only rapidly evolving but also conserved genes, and functionally similar proteins may be represented and relevant in the seminal fluid of other animals [4]. Members of protein classes conserved across organisms in seminal fluids may reflect their general importance in reproductive physiology [4] and might therefore be beneficial for both males and females. In this context, immune defense proteins are of special interest. MAGs secrete proteins of antimicrobial activity [10] and proteins, which regulate the expression of genes with antimicrobial activity in females after sperm transfer [11,12,47]. Such Acps might successively protect sperm from microbial attack in the reproductive tracts of both sexes, as well as coat eggs when they are laid. A L. gredleri EST, whose ortholog CG6910 is known to be upregulated in D. melanogaster females after mating [47], might play such a role. The Drosophila gene is annotated as possibly showing oxidoreductase activity and being downregulated after bacterial infection. For the elucidation of its function in ants, it should be of interest to see whether the gene similar to CG6910 is also differentially expressed in virgin and mated L. gredleri queens. Females are not 'silent partners' when it comes to seminal fluids and it is always an interplay of male and female factors that coordinates the gametes for fertilization [1,28,48].

The CG32432 protein predicted in the honey bee genome as a Drosophila ortholog with undefined function was found to be significantly overexpressed in sexually mature L. gredleri males, and one gene similar to CG1486-PA was found to be overexpressed in two male life cycle stages. This latter protein is predicted to have a carboxy-lyase activity and might be involved in carboxylic acid metabolism, catalyzing the non-hydrolytic addition or removal of a carboxyl group to or from a compound. Carboxylic acids are widespread in nature, giving rise to fatty acids. Fatty acids, such as palmitic-, linoleic-, oleic-, and stearic-acids, have been described as components of the mating plug in the bumblebee Bombus terrestris, which form a physical barrier [30] and prevent re-mating [31]. Mating plugs were also found in fire ants [32], queenless Dinoponera ants [49], attine ants [50], stingless bees [51], and in the honey bee. In the latter, however, they do not play a function as a plug, but rather as a mating sign that may even promote multiple matings [52]. In L. gredleri, the existence of a spermatophore that could be comparable to a mating plug was already shown [42]. The overexpression of a gene with carboxy-lyase activity in the late pupal stage, when the accessory glands are already formed, and in sexually mature males found outside the nest supports the idea that accessory gland products are crucial for spermatophore formation in this species.

Taken together, our current study represents a first step addressing at a molecular level the reproductive interplay between L. gredleri males and their mates. MAG products identified in this study may play different roles, either as conserved or as species-specific proteins.

Understanding the molecular underpinnings of differences in reproductive strategies among non-social and eusocial species, and in the latter among monogamous (the majority of ants bees and wasps) and polygamous ones (e.g. the honey bees) will still require sequencing of MAG products in a range of species. Nevertheless, we believe that at this interface of "molecular communication" between the sexes, much can be gained towards comprehending the evolution of eusociality at the passage through the monogamy window [25,28], with inter-sexual conflict being replaced by inter-sexual cooperation.

Methods

Colony collection and male dissection

Colonies of Leptothorax gredleri Mayr 1855 (Hymenoptera: Formicidae) were collected from nest-sites at the edge of an abandoned army drill ground in Erlangen, Germany (49° 35' 09" N, 11° 02' 02" E). Each single colony was transferred into a three-chamber plastic box with a cavity between two microscope slides serving as a nest. Colonies were kept under standard rearing conditions [53,54] until they produced sexuals.

Males were dissected in chilled Beadle solution (128.3 mM NaCl, 4.7 mM KCl, and 2.3 mM CaCl2; [55]). The accessory glands were separated and the seminal ducts were removed. The intestine was removed to avoid contamination of libraries with microorganismal RNA and then the rest of the body served as reference tissue. Glands and the remaining body carcass were stored separately in RNAlater® (Ambion) for subsequent RNA extraction. For library preparations, 40 pairs of glands were pooled for a single sample, to be subtracted against the carcasses of 20 bodies.

For qRT-PCR analyses we dissected male pupae ready to eclose, when they are already dark pigmented (male pupae), newly eclosed males collected from within the nest (males inside nests), and males running outside the nest that were sexually mature (males outside nests). Each sample consisted of at least 20 pairs of glands, to be compared against corresponding body carcasses. As the real time RT-PCR analyses were carried out as quintuplicates, this represented biological material of more than 100 males per life cycle stage.

Representational difference analysis

RNA of both the glands and the remaining body from males (collected from within nests) was extracted with TRIzol (Invitrogen) according to the manufacturer's standard protocol. We used 2 μg of each RNA sample to generate the cDNA library. Reverse transcription and long distance PCR were carried out by employing the SMART PCR cDNA synthesis kit (Clontech), resulting in double-stranded cDNA.

For suppression subtractive hybridization we employed a cDNA RDA [56] adapted for applications in insects [57]. This strategy removes transcripts shared by the driver and tester populations and enriches the tester library for differentially expressed genes. An accessory gland library was generated using cDNA of accessory glands as the tester population. The latter was hybridized against cDNA of the remaining body carcass as the driver population. Double-stranded cDNAs (1 μg) were restriction digested with MboI (New England BioLabs), ligated to R-adapters [57], and PCR amplified following the protocol of Hubank and Schatz [58] to generate the respective cDNA representations. Enrichment of differentially expressed transcripts was achieved in two successive rounds of PCR amplifications that employed different adapters (J and N) and sequential subtractive hybridizations of the tester to an excess of driver cDNA, first in a 1:100 and then in a 1:800 ratio [58]. After each of these successive steps, cDNAs were purified (Illustra GFX PCR purification kit, GE).

cDNA from the accessory gland library after the second round of selection was used for ligation into pGEM®-T Easy Vector System (Promega) and transformation of E. coli DH5α chemocompetent cells. The cells were reared on solid LB containing ampicillin, X-Gal and IPTG. Selected clones were grown in liquid LB medium containing ampicillin. For the sequencing, the cells were lysed and the plasmids extracted and sequenced using Big Dye Terminator Cycle Sequencing Ready Reaction (Applied Biosystems) with M13 primers on an ABI-PRISM 3100 (Applied Biosystems) automated gene analyzer.

Bioinformatic analysis

By means of the E-Gene annotation pipeline [59], sequencing reads were first filtered to detect and remove ribosomal RNAs and mitochondrial DNA. Subsequently, vector sequences were trimmed using Crossmatch. Read quality was checked and reads were assembled by the Phred-Phrap program module. Reads that had passed the quality check were next submitted to CAP3 assembly [60] to obtain USs. All contigs, and singlets were dynamically translated and compared by BLASTx to a non-redundant (nr) database (GenBank). ESTs with similarities to genes of known function in other organisms were clustered using Gene Ontology (GO) terms [61] attributed to their respective Drosophila orthologs. GO attributes were retrieved from the Flybase server http://flybase.org webcite.

Differential gene expression analysis by quantitative RT-PCR

For real time qRT-PCR analysis, gene-specific primers were designed for a set of selected USs using Gene Runner version 3.05 and Primer3 softwares.

As reference gene we chose elongation factor 1-alpha (elf1-α), a recommended control gene for honey bee qRT-PCR studies [62]. To obtain an L. gredleri sequence for elf1-α we aligned the corresponding sequences of Leptothorax muscorum (a direct sister species to L. gredleri) and Temnothorax rugatulus (a more distantly related ant species) (GenBank accession numbers ABK54792.1 and ABK54711.1, respectively) and designed primers to a conserved region. An L. gredleri elf1-α fragment was cloned and sequenced, and qRT-PCR suitable primers were designed (Table 2). For the qRT-PCR analyses, we used male pupae, males inside nests and males outside nests. Each sample consisted of a least 20 male bodies or glands. RNA of the glands and remaining bodies carcasses was extracted using TRIzol. The RNA samples were treated with 0.1 U DNaseI (Invitrogen) for 40 min to eliminate possible DNA contamination. First strand cDNA was produced using a SuperScript II (Invitrogen) protocol at 42°C for 50 min and 70°C for 15 min.

Quality and annealing temperatures of the gene-specific primers were tested in a temperature-gradient PCR protocol run in a PTC200 thermal cycler (MJ Research). Product length varied from 305 bp for elf1-α to 98 bp for the shortest of the L. gredleri singlets.

Subsequently, a dilution series was analyzed for each primer pair using a composite cDNA sample from bodies and glands of all stages in order to obtain standard curves for determining primer efficiencies. These and subsequent qRT-PCR amplifications were performed by using a SYBR Green (Applied Biosystems) protocol in an ABI Prism 7500 system (Applied Biosystems). The amplification protocol was 50°C for 2 min, 95°C for 10 min, 40 cycles of 95°C for 15 s and 60°C for 1 min. Subsequently, dissociation curves were acquired to check melting peak quality. For each life-cycle time point we ran five biological replicates, each assayed in technical triplicates, both for accessory glands and body carcasses.

Threshold cycle (Ct) values for each technical replicate were used to calculate 2-ΔΔCt values [63] using the means of the five body carcass replicates of each life cycle stage for normalization of the respective gland samples. For statistical analysis, Ct values were used as input to randomization tests implemented in Relative Expression Software Tool (REST - [64]). P values ≤ 0.05 were considered as statistically significant.

Authors' contributions

JH and KH conceived the study. KH, JH and AO designed the study. AO and MF collected and reared the males, dissected them and extracted the RNA. FCH and AO conducted the RDA, cloned and sequenced the fragments and conducted the bioinformatics analysis. SVA and AO designed the primers, MF, AO, FCH, SVA and DSMA tested the primers and conducted qRT-PCR. FCH analyzed the qRT-PCR results. AO, FCH, KH, and JH wrote the manuscript. AO and FCH contributed equally to the study. All authors read and approved the final manuscript.

Acknowledgements

We thank Alexandra Schrempf and Andreas Trindl for helpful advice, as well as the multiuser sequencing and real-time PCR facility of the Departamento de Biologia Celular e Molecular e Bioagentes Patogênicos da FMRP-USP for technical support. This study was supported by DFG (He 1623/19) and by a bilateral CAPES/DAAD grant (PROBRAL project 261/07).

FCH and SVA are sponsored by FAPESP and DSMA by CNPq doctoral fellowships.

References

  1. Gillott C: Male accessory gland secretions: Modulators of female reproductive physiology and behavior.

    Annual Review of Entomology 2003, 48:163-184. PubMed Abstract | Publisher Full Text OpenURL

  2. Wolfner MF: The gifts that keep on giving: Physiological functions and evolutionary dynamics of male seminal proteins in Drosophila.

    Heredity 2002, 88:85-93. PubMed Abstract | Publisher Full Text OpenURL

  3. Poiani A: Complexity of seminal fluid: a review.

    Behavioral Ecology and Sociobiology 2006, 60:289-310. Publisher Full Text OpenURL

  4. Ram KR, Wolfner MF: Sustained post-mating response in Drosophila melanogaster requires multiple seminal fluid proteins.

    PLoS Genetics 2007, 3:2428-2438. Publisher Full Text OpenURL

  5. Ram KR, Wolfner MF: Seminal influences: Drosophila Acps and the molecular interplay between males and females during reproduction.

    Integrative and Comparative Biology 2007, 47:427-445. Publisher Full Text OpenURL

  6. Chen PS, Stumm-Zollinger E, Aigaki T, Balmer J, Bienz M, Bohlen P: A male accessory gland peptide that regulates reproductive behavior of female Drosophila melanogaster.

    Cell 1988, 54:291-298. PubMed Abstract | Publisher Full Text OpenURL

  7. Kalb JM, Dibenedetto AJ, Wolfner MF: Probing the function of Drosophila melanogaster accessory glands by directed cell ablation.

    Proceedings of the National Academy of Sciences USA 1993, 90:8093-8097. Publisher Full Text OpenURL

  8. Chapman T, Bangham J, Vinti G, Seifried B, Lung O, Wolfner MF, Smith HK, Partridge L: The sex peptide of Drosophila melanogaster: Female post-mating responses analyzed by using RNA interference.

    Proceedings of the National Academy of Sciences USA 2003, 100:9923-9928. Publisher Full Text OpenURL

  9. Liu H, Kubli E: Sex-peptide is the molecular basis of the sperm effect in Drosophila melanogaster.

    Proceedings of the National Academy of Sciences USA 2003, 100:9929-9933. Publisher Full Text OpenURL

  10. Lung O, Kuo L, Wolfner MF: Drosophila males transfer antibacterial proteins from their accessory gland and ejaculatory duct to their mates.

    Journal of Insect Physiology 2001, 47:617-622. PubMed Abstract | Publisher Full Text OpenURL

  11. McGraw LA, Gibson G, Clark AG, Wolfner MF: Genes regulated by mating, sperm, or seminal proteins in mated female Drosophila melanogaster.

    Current Biology 2004, 14:1509-1514. PubMed Abstract | Publisher Full Text OpenURL

  12. Peng J, Zipperlen P, Kubli E: Drosophila sex-peptide stimulates female innate immune system after mating via the Toll and Imd pathways.

    Current Biology 2005, 15:1690-1694. PubMed Abstract | Publisher Full Text OpenURL

  13. Neubaum DM, Wolfner MF: Mated Drosophila melanogaster females require a seminal fluid protein, Acp36DE, to store sperm efficiently.

    Genetics 1999, 153:845-857. PubMed Abstract | PubMed Central Full Text OpenURL

  14. Tram U, Wolfner MF: Male seminal fluid proteins are essential for sperm storage in Drosophila melanogaster.

    Genetics 1999, 153:837-844. PubMed Abstract | PubMed Central Full Text OpenURL

  15. Xue L, Noll M: Drosophila female sexual behavior induced by sterile males showing copulation complementation.

    Proceedings of the National Academy of Sciences USA 2000, 97:3272-3275. Publisher Full Text OpenURL

  16. Chapman T, Liddle LF, Kalb JM, Wolfner MF, Partridge L: Cost of mating in Drosophila melanogaster females is mediated by male accessory gland products.

    Nature 1995, 373:241-244. PubMed Abstract | Publisher Full Text OpenURL

  17. Lung O, Tram U, Finnerty CM, Eipper-Mains MA, Kalb JM, Wolfner MF: The Drosophila melanogaster seminal fluid protein Acp62F is a protease inhibitor that is toxic upon ectopic expression.

    Genetics 2002, 160:211-224. PubMed Abstract | PubMed Central Full Text OpenURL

  18. Chapman T, Arnqvist G, Bangham J, Rowe L: Sexual conflict.

    Trends in Ecology & Evolution 2003, 18:41-47. OpenURL

  19. Rice WR, Holland B: The enemies within: Intergenomic conflict, interlocus contest evolution (ICE), and the intraspecific Red Queen.

    Behavioral Ecology and Sociobiology 1997, 41:1-10. Publisher Full Text OpenURL

  20. Clark AG, Begun DJ, Prout T: Female x male interactions in Drosophila sperm competition.

    Science 1999, 283:217-220. PubMed Abstract | Publisher Full Text OpenURL

  21. Rice WR: Dangerous liaisons.

    Proceedings of the National Academy of Sciences USA 2000, 97:12953-12955. Publisher Full Text OpenURL

  22. Swanson WJ, Clark AG, Waldrip-Dail HM, Wolfner MF, Aquadro CF: Evolutionary EST analysis identifies rapidly evolving male reproductive proteins in Drosophila.

    Proceedings of the National Academy of Sciences USA 2001, 98:7375-7379. Publisher Full Text OpenURL

  23. Mueller JL, Ram KR, McGraw LA, Qazi MCB, Siggia ED, Clark AG, Aquadro CF, Wolfner MF: Cross-species comparison of Drosophila male accessory gland protein genes.

    Genetics 2005, 171:131-143. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Boomsma JJ: Kin selection versus sexual selection: why the ends do not meet.

    Current Biology 2007, 17:R673-R678. PubMed Abstract | Publisher Full Text OpenURL

  25. Boomsma JJ: Lifetime monogamy and the evolution of eusociality.

    Philosophical Transactions of the Royal Society of London B 2009, 364:3191-3207. Publisher Full Text OpenURL

  26. Boomsma JJ, Baer B, Heinze J: The evolution of male traits in social insects.

    Annual Review of Entomology 2005, 50:395-420. PubMed Abstract | Publisher Full Text OpenURL

  27. Hughes WHO, Oldroyd BP, Beekman M, Ratnieks FLW: Ancestral monogamy shows kin selection is key to the evolution of eusociality.

    Science 2008, 320:1213-1216. PubMed Abstract | Publisher Full Text OpenURL

  28. den Boer SPA, Baer B, Boomsma JJ: Seminal fluid mediates ejaculate competition in social insects.

    Science 2010, 327:1506-1509. PubMed Abstract | Publisher Full Text OpenURL

  29. Schrempf A, Heinze J, Cremer S: Sexual cooperation: Mating increases longevity in ant queens.

    Current Biology 2005, 15:267-270. PubMed Abstract | Publisher Full Text OpenURL

  30. Baer B, Maile R, Schmid-Hempel P, Morgan ED, Jones GR: Chemistry of a mating plug in Bumblebees.

    Journal of Chemical Ecology 2000, 26:1869-1875. Publisher Full Text OpenURL

  31. Baer B, Morgan ED, Schmid-Hempel P: A nonspecific fatty acid within the Bumblebee mating plug prevents females from remating.

    Proceedings of the National Academy of Sciences USA 2001, 98:3926-3928. Publisher Full Text OpenURL

  32. Mikheyev AS: Evidence for mating plugs in the fire ant Solenopsis invicta.

    Insectes Sociaux 2003, 50:401-402. Publisher Full Text OpenURL

  33. den Boer SPA, Boomsma JJ, Baer B: Seminal fluid enhances sperm viability in the leafcutter ant Atta colombica.

    Behavioral Ecology and Sociobiology 2008, 62:1843-1849. Publisher Full Text OpenURL

  34. Oppelt A, Heinze J: Mating is associated with immediate changes of the hydrocarbon profile of Leptothorax gredleri ant queens.

    Journal of Insect Physiology 2009, 55:624-628. PubMed Abstract | Publisher Full Text OpenURL

  35. Collins AM, Caperna TJ, Williams V, Garrett WM, Evans JD: Proteomic analyses of male contributions to honey bee sperm storage and mating.

    Insect Molecular Biology 2006, 15:541-549. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Frattini-Colonello NA, Hartfelder K: Differential gene expression profiling in mucus glands of honey bee (Apis mellifera) drones during sexual maturation.

    Apidologie 2009, 40:481-495. Publisher Full Text OpenURL

  37. Frattini-Colonello NA, Hartfelder K: She's my girl - male accessory gland products and their function in the reproductive biology of social bees.

    Apidologie 2005, 36:231-244. Publisher Full Text OpenURL

  38. Baer B, Heazlewood J, Taylor NL, Eubel H, Millar AH: The seminal fluid proteome of the honeybee Apis mellifera.

    Proteomics 2009, 9:2085-2097. PubMed Abstract | Publisher Full Text OpenURL

  39. Heinze J, Lipski N, Hölldobler B: Reproductive competition in colonies of the ant Leptothorax gredleri.

    Ethology 1992, 90:265-278. Publisher Full Text OpenURL

  40. Oberstadt B, Heinze J: Mating biology and population structure of the ant, Leptothorax gredleri.

    Insectes Sociaux 2003, 50:340-345. Publisher Full Text OpenURL

  41. Oppelt A, Spitzenpfeil N, Kroiss J, Heinze J: The significance of intercolonial variation of cuticular hydrocarbons for inbreeding avoidance in ant sexuals.

    Animal Behaviour 2008, 76:1029-1034. Publisher Full Text OpenURL

  42. Oppelt A, Heinze J: Dynamics of sperm transfer in the ant Leptothorax gredleri.

    Naturwissenschaften 2007, 94:781-786. PubMed Abstract | Publisher Full Text OpenURL

  43. Nunes FMF, Valente V, Sousa JF, Cunha MAV, Pinheiro DG, Maia RM, Araujo DD, Costa MCR, Martins WK, Carvalho AF, Monesi N, Nascimento AM, Peixoto PMV, Silva MFR, Ramos RGP, Reis LFL, Dias-Neto E, Souza SJ, Simpson AJG, Zago MA, Soares AEE, Bitondi MMG, Espreafico EM, Espindola FS, Paco-Larson ML, Simoes ZLP, Hartfelder K, Silva WA Jr: The use of Open Reading frame ESTs (ORESTES) for analysis of the honey bee transcriptome.

    BMC Genomics 2004, 5:e84. BioMed Central Full Text OpenURL

  44. Willems E, Leyns L, Vandesompele J: Standardization of real-time PCR gene expression data from independent biological replicates.

    Analytical Biochemistry 2008, 379:127-129. PubMed Abstract | Publisher Full Text OpenURL

  45. Harwood R: Protein transfer across membranes: the role of signal sequences and signal peptidase activity. In The enzymology of post-translational modification of proteins. Edited by Freedman RB, Hawkins HC. London: Academic Press; 1980:3-52. OpenURL

  46. Hammond RL, Bourke AFG, Bruford MW: Mating frequency and mating system of the polygynous ant, Leptothorax acervorum.

    Molecular Ecology 2001, 10:2719-2728. PubMed Abstract | Publisher Full Text OpenURL

  47. Lawniczak MKN, Begun DJ: A genome-wide analysis of courting and mating responses in Drosophila melanogaster females.

    Genome 2004, 47:900-910. PubMed Abstract | Publisher Full Text OpenURL

  48. den Boer SPA, Boomsma JJ: Honey bee males and queens use glandular secretions to enhance sperm viability before and after storage.

    Journal of Insect Physiology 2009, 55:538-543. PubMed Abstract | Publisher Full Text OpenURL

  49. Monnin T, Peeters C: Monogyny and regulation of worker mating in the queenless ant Dinoponera quadriceps.

    Animal Behavior 1998, 55:299-306. Publisher Full Text OpenURL

  50. Baer B, Boomsma JJ: Male reproductive investment and queen mating-frequency in fungus-growing ants.

    Behavioral Ecology 2004, 15:426-432. Publisher Full Text OpenURL

  51. Kerr WE, Zucchi R, Nakadaira JT, Butolo JE: Reproduction in the social bees (Hymenoptera: Apidae).

    Journal of New York Entomological Society 1962, 70:265-276. OpenURL

  52. Woyciechowski M, Kabat L, Król E: The function of the mating signs in honeybees, Apis mellifera L.: new evidence.

    Animal Behaviour 1994, 47:733-735. Publisher Full Text OpenURL

  53. Buschinger A: Experimente und Beobachtungen zur Gründung und Entwicklung neuer Sozietäten der sklavenhaltenden Ameise Harpagoxenus sublaevis (Nyl.).

    Insectes Sociaux 1974, 21:381-406. Publisher Full Text OpenURL

  54. Heinze J, Ortius D: Social organization of Leptothorax acervorum from Alaska (Hymenoptera: Formicidae).

    Psyche (Cambridge) 1991, 98:227-240. OpenURL

  55. Darrouzet E, Huignard J, Chevrier C: Effect of differential host exposure on reproduction in the solitary ectoparasitoid Eupelmus vuilleti.

    Entomologia Experimentalis et Applicata 2002, 103:73-81. Publisher Full Text OpenURL

  56. Pastorian K, Hawel III L, Byus CV: Optimization of cDNA representational difference analysis for the identification of differentially expressed mRNAs.

    Analytical Biochemistry 2000, 283:89-98. PubMed Abstract | Publisher Full Text OpenURL

  57. Judice CC, Carazzole MF, Festa F, Sogayar MC, Hartfelder K, Pereira GAG: Gene expression profiles underlying alternative caste phenotypes in a highly eusocial bee, Melipona quadrifasciata.

    Insect Molecular Biology 2006, 15:33-44. PubMed Abstract | Publisher Full Text OpenURL

  58. Hubank M, Schatz DG: Representational difference analysis of cDNA. In Functional Genomics. Edited by Hunt SP, Livesey FJ. Oxford, New York: University Oxford Press; 2000:45-80. OpenURL

  59. Durham AM, Kashiwabara AY, Matsunaga FTG, Ahagon PH, Rainone F, Gruber A: EGene: a configurable pipeline generation system for automated sequence analysis.

    Bioinformatics 2005, 21:2812-1813. PubMed Abstract | Publisher Full Text OpenURL

  60. Huang X, Madan A: CAP3: A DNA sequencing assembly program.

    Genome Research 1999, 9:868-877. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. The Gene Ontology Consortium: Gene Ontology: tool for the unification of biology.

    Nature Genetic 2000, 25:25-29. Publisher Full Text OpenURL

  62. Lourenço AP, Mackert A, Cristino AS, Simoes ZLP: Validation of reference genes for gene expression studies in the honey bee, Apis mellifera, by quantitative real-time RT-PCR.

    Apidologie 2008, 39:378-385. Publisher Full Text OpenURL

  63. Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR.

    Nucleic Acids Research 2001, 29:e45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST©) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.

    Nucleic Acids Research 2002, 30:1-10. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL