Vitellogenin Underwent Subfunctionalization to Acquire Caste and Behavioral Specific Expression in the Harvester Ant Pogonomyrmex barbatus.

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.

Extended family structure in the ant Formica paralugubris: the role of the breeding system

In populations of various ant species, many queens reproduce in the same nest (polygyny), and colony boundaries appear to be absent with individuals able to move fi eely between nests (unicoloniality). Such societies depart strongly from a simple family structure and pose a potential challenge to kin selection theory, because high queen number coupled with unrestricted gene flow among nests should result in levels of relatedness among nestmates close to zero. This study investigated the breeding system and genetic structure of a highly polygynous and largely unicolonial population of the wood ant Formica paralugubris. A microsatellite analysis revealed that nestmate workers, reproductive queens and reproductive males (the queens' mates) are all equally related to each other, with relatedness estimates centring around 0.14. This suggests that most of the queens and males reproducing in the study population had mated within or close to their natal nest, and that the queens did not disperse far after mating. We developed a theoretical model to investigate how the breeding system affects the relatedness structure of polygynous colonies. By combining the model and our empirical data, it was estimated that about 99.8% of the reproducing queens and males originated from within the nest, or from a nearby nest. This high rate of local mating and the rarity of long-distance dispersal maintain significant relatedness among nestmates, and contrast with the common view that unicoloniality is coupled with unrestricted gene flow among nests.

Single loci with major effects in fire ants : the effects of the complementary sex-determining locus and the green beard supergene on gene expression

Etant données la complexité et la redondance des réseaux de gènes influençant de nombreux phénotypes, l'étude des rares cas d'un locus unique ayant des effets importants sur de nombreux phénotypes peut fournir des informations cruciales sur l'évolution des traits complexes. Nous avons séquencé le génome de la fourmi de feu Solenopsis invicta pour étudier comment l'expression des gènes détermine les effets majeurs et étendus de deux loci uniques sur le phénotype. Le premier locus concerne la détermination du sexe par le modèle des allèles complémentaires. Ce locus est connu pour déterminer le sexe chez tous les hyménoptères mais n'a été caractérisé que chez les abeilles. Les hétérozygotes pour ce locus se développent en reines diploïdes (ou ouvrières stériles) alors que les homozygotes se développent en mâles diploïdes incapables de produire du sperme et les hémizygotes en mâles haploïdes fertiles. Nous avons comparé l'expression des gènes entre les reines et les deux types de mâles au stade pupe, ainsi que 1 et 11 jours après l'émergence. Nous avons trouvé un changement prononcé de l'expression des gènes chez les mâles diploïdes, passant de très proche de celle des reines au stade pupe à identique aux mâles haploïdes 11 jours après l'émergence. Cela signifie que les mâles diploïdes sont condamnés à être stériles parce que les effets après émergence du locus de détermination du sexe ne per¬mettent pas d'effacer les effets de la ploïdie sur l'expression des gènes pendant le stade pupe, quand la spermatogénèse prend place. Le second locus aux effets majeurs que nous avons étudié est le supergène dit "green beard", qui consiste en 616 gènes couvrant 55% d'un chromosome (13 Mb) et est caractérisé par une absence de recombinaison entre les deux variants du supergène : "Social B" et "Social b" (SB et Sb). Au travers de l'effet "green beard", par lequel les ouvrières avec le supergène Sb discriminent favorablement les reines qui partagent ce supergène de façon perceptible, le génotype des reines fondatrices au niveau de ce supergène détermine l'organisation de la colonie : soit elle contient une seule reine SB/SB, soit plusieurs reines SB/Sb. Nous avons montré que le chromosome Sb a évolué comme le chromosome Y, accumulant probablement des allèles favorables dans des colonies avec plusieurs reines mais défavorables dans des colonies avec une seule reine (cf. gènes sexuellement antagonistes), ainsi que des transposons et des séquences répéti¬tives. Nous avons également montré que le polymorphisme du supergène cause de grandes différences d'expression chez les ouvrières et particulièrement les reines mais pas chez les mâles. Pour comprendre comment le polymorphisme du supergène chez les reines peut affecter l'organisation de la colonie, nous avons comparé l'expression entre les génotypes SB/SB et SB/Sb chez des reines vierges (1 et 11 jours) et des reines matures. Nous avons montré que les reines SB/SB sur-régulent des gènes impliqués dans la reproduction, expli-quant pourquoi elle grandissent plus rapidement et peuvent fonder des colonies de façon indépendante, tandis que les reines SB/Sb (qui ne peuvent fonder une nouvelle colonie) sur-régulent des gènes de signalement chimique qui affectent l'organisation des colonies par l'effet "green beard". - Given the complexity and redundancy of the gene networks that underlie many pheno- types, the study of rare cases of a single locus having major effects on many phenotypes can give powerful insights into the evolution of complex traits. We sequenced the genome of Solenopsis invicta fire ants to study how gene expression mediates the widespread major effects of two single loci on phenotype. The first is the complementary sex-determining locus, which is known to exist in most Hymenoptera despite being characterized only for honeybees. Heterozygotes at this locus become diploid queens (or sterile workers), homozy¬gotes become aspermic diploid males, and hemizygotes become fertile haploid males. We compared gene expression between queens and both types of males in pupae and 1 and 11 days after eclosion. We found a pronounced shift in gene expression in diploid males, from being nearly identical to queens as pupae to identical to haploid males 11 days after eclosion. This means that diploid males are condemned to sterility because the overriding effects of the sex locus after eclosion cannot undo the ploidy effects on expression during the pupal stage, when spermatogenesis must be completed. The second locus with major ef¬fects that we studied was the so-called "green beard" supergene, which consists of 616 genes encompassing 55% of one chromosome (13 Mb), without recombination between the two variants "Social B" and "Social b" (SB and Sb) supergene. Through the green beard effect, i.e. workers with the Sb supergene discriminating in favor of queens who perceptibly share this supergene, the founding queen's genotype at the supergene determines colony organi¬zation: either headed by a single SB/SB queen or many SB/Sb queens. We show that the Sb chromosome evolved like a Y-chromosome, probably accumulating alleles beneficial in multi-queen colonies but disadvantageous in single-queen colonies (cf. sexually antagonistic genes), as well as transposons and repetitive sequences. We also show that the polymor¬phism of the supergene causes widespread expression differences in workers and especially queens but not in males. To understand how the polymorphism at the supergene in queen can transform colony organization, we compared the expression between SB/SB and SB/Sb virgin queens (1 and 11 days) and mother queens. We show that SB/SB queens up-regulate genes involved in reproduction, explaining why they mature faster and can found colonies independently, while SB/Sb queens (which cannot found colonies) up-regulate chemical signaling genes that can transform colonies through the green beard effect.

Riqueza de formigas de solo na praia da Pedreira, Parque Estadual de Itapuã, Viamão, RS, Brasil

Neste trabalho são listadas as espécies de formigas de solo encontradas em três ambientes (mata nativa, barreira pedregosa e areias da orla) na praia da Pedreira, localizada na Zona de Uso Intensivo do Parque Estadual de Itapuã, município de Viamão, RS. No total, foram identificadas 60 espécies representantes de 24 gêneros, 18 tribos e oito subfamílias (Dolichoderinae, Ecitoninae, Ectatomminae, Formicinae, Heteroponerinae, Myrmicinae, Ponerinae e Pseudomyrmecinae). Dos três ambientes amostrados, de acordo com o estimador de riqueza jackknife de primeira ordem, a mata nativa apresentou a maior riqueza (Sest= 37,5), seguindo-se a barreira pedregosa (Sest= 8,9) e as areias da orla (Sest= 5,9). Apenas Acromyrmex laticeps, Crematogaster sp. e Solenopsis invicta foram comuns aos três ambientes. Neste trabalho é feito o primeiro registro de ocorrência de Pachycondyla crenata e Pachycondyla laevigata (Ponerinae) para o Rio Grande do Sul.

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

Measuring ecological niche overlap from occurrence and spatial environmental data

1. Aim - Concerns over how global change will influence species distributions, in conjunction with increased emphasis on understanding niche dynamics in evolutionary and community contexts, highlight the growing need for robust methods to quantify niche differences between or within taxa. We propose a statistical framework to describe and compare environmental niches from occurrence and spatial environmental data.¦2. Location - Europe, North America, South America¦3. Methods - The framework applies kernel smoothers to densities of species occurrence in gridded environmental space to calculate metrics of niche overlap and test hypotheses regarding niche conservatism. We use this framework and simulated species with predefined distributions and amounts of niche overlap to evaluate several ordination and species distribution modeling techniques for quantifying niche overlap. We illustrate the approach with data on two well-studied invasive species.¦4. Results - We show that niche overlap can be accurately detected with the framework when variables driving the distributions are known. The method is robust to known and previously undocumented biases related to the dependence of species occurrences on the frequency of environmental conditions that occur across geographic space. The use of a kernel smoother makes the process of moving from geographical space to multivariate environmental space independent of both sampling effort and arbitrary choice of resolution in environmental space. However, the use of ordination and species distribution model techniques for selecting, combining and weighting variables on which niche overlap is calculated provide contrasting results.¦5. Main conclusions - The framework meets the increasing need for robust methods to quantify niche differences. It is appropriate to study niche differences between species, subspecies or intraspecific lineages that differ in their geographical distributions. Alternatively, it can be used to measure the degree to which the environmental niche of a species or intraspecific lineage has changed over time.

Ecology and evolution of social-organization: insights from fire ants and other highly eusocial insects

Social organisms exhibit conspicuous intraspecific variation in all facets of their social organization. A prominent example of such variation in the highly eusocial Hymenoptera is differences in the number of reproductive queens per colony, Differences in queen number in ants are associated with differences in a host of reproductive and social traits, including queen phenotype and breeding strategy, mode of colony reproduction, and pattern of sex allocation. We examine the causes and consequences of changes in colony queen number and associated traits using the fire ant Solenopsis invicta as a principal model. Ecological constraints on mode of colony founding may act as important selective forces causing the evolution of queen number in this and many other ants, with social organization generally perpetuated across generations by means of the social environment molding appropriate queen phenotypes and reproductive strategies. Shifts in colony queen number have profound effects on genetic structure within nests and may also influence genetic structure at higher levels (aggregations of nests or local demes) because of the association of queen number with particular mating and dispersal habits. Divergence of breeding habits between populations with different social organizations has the potential to promote genetic differentiation between these social variants. Thus, evolution of social organization can be important in generating intrinsic selective regimes that channel subsequent social evolution and in initiating the development of significant population genetic structure, including barriers to gene flow important in cladogenesis.

Sociogenetics of fire ants

RÉSUMÉ GRAND PUBLIC La complexité des sociétés d'insectes (telles que les abeilles, les termites ou les fourmis) a depuis longtemps fasciné l'Homme. Depuis le débfit du XIXème siècle, de nombreux travaux observationnels, comportementaux et théoriques leur on été consacrés afin de mieux les décrire et comprendre. L'avènement de la biologie moléculaire à la fin du XXèrne siècle a offert de nouveaux outils scientifiques pour identifier et étudier les gènes et molécules impliqués dans le développement et le comportement des êtres vivants. Alors que la majorité de ces études s'est focalisée sur des organismes de laboratoire tel que la mouche ou les nématodes, l'utilisation de ces outils est restée marginale jusqu'à présent dans l'étude des sociétés d'insectes. Lors de ma thèse, j'ai développé des outils moléculaires permettant de déterminer le niveau d'activité de zo,ooo gènes chez la fourmi de feu, Solenopsis invicta, ainsi qu'une base de données et un portail en ligne regroupant les informations relatives à l'étude génétique des fourmis: Fourmidable. J'ai ensuite utilisé ces outils dans le cadre d'une étude comportementale chez la fourmis S. invicta. Dans les sociétés d'insectes, une hiérarchie peut déterminer le statut reproducteur des individus. Suite à la mort d'un dominant, les subordonnés entrent en compétition en vue d'améliorer leur statut. Un tel phénomène se produit au sein des colonies de S. invicta contenant une unique reine mère, des milliers d'ouvrières et des centaines de reines vierges ailées. A la mort de la reine mère, un grand nombre de reines vierges tentent de la remplacer en arrachant leurs ailes et en activant leurs organes reproducteurs plutôt que de partir en vol nuptial. Ces tentatives sont le plus souvent arrêtées par les ouvrières qui exécutent la plupart de ces reines sur la base de signaux olfactifs produits lors de l'activation des organes reproducteurs. Afin de mieux comprendre les mécanismes moléculaires impliqués, j'ai étudié l'activité de gènes au sein des reines au début de ce processus. J'ai ainsi déterminé que des gènes impliqués dans communication olfactive, le développement des organes reproducteurs et la métabolisation de l'hormone juvénile sont activês à ce moment là. La vitesse à laquelle les reines perdent leurs ailes ainsi que les niveaux d'expression de gènes sont ensuite liés à leur probabilité de survie. ABSTRACT : Honeybees, termites and ants occupy the "pinnacle of social evolution" with societies of a complexity that rivals our own. Humans have long been fascinated by social insects, but studying them has been mostly limited to observational and behavioral experiments. The advent of molecular biology first made it possible to investigate the molecular-genetic basis of development in model systems such as the fruit fly Drosophila melarcogaster or the roundworm Caenorhabditis elegans and subsequently their behavior. Molecular and genomic tools are now becoming available for the study of social insects as well. To permit genomic research on the fire ant, Solenopsis invicta, we developed a cDNA microarray that can simultaneously determine the expression levels of approximately 1oooo genes. These genes were assembled and bioinformatically annotated using custom pipelines. The obtained data formed the cornerstones for Fourmidable, a web portal centralizing sequence, gene annotation and gene expression data as well as laboratory protocols for research on ants. In many animals living in groups the reproductive status of individuals is determined by their social status. In species with social hierarchies, the death of dominant individuals typically upheaves the social hierarchy and provides an opportunity for subordinate individuals to improve their social status. Such a phenomenon occurs in the monogyne form of S. invicta, where colonies typically contain a single wingless reproductive queen, thousands of workers and hundreds of winged non-reproductive virgin queens. Upon the death of the mother queen, many virgin queens shed their wings and initiate reproductive development instead of departing on a mating flight. Workers progressively execute almost all of them over the following weeks. The workers base their collective decision on pheromonal cues associated with the onset of reproductive development of the virgin queens which occurs after orphaning. We used the aforementioned tools to determine that genes putatively involved in processes including olfactory signaling, reproductive development and Juvenile Hormone metabolism are differentially expressed at the onset of competition. Additionally, we found that queens that initiate reproductive development faster and, to a certain extent, shed their wings faster after orphaning are more likely to become replacement queens. These results provide candidate genes that are putatively linked to competition outcome. To determine the extent to which specific genes affect different aspects of life in ant colonies, functional tests such as gene activation and silencing will still be required. We conclude by discussing some of the challenges and opportunities for molecular-genetic research on ants. RÉSUMÉ Les sociétés d'abeilles, de termites et de fourmis sont d'une complexité proche de celle de la nôtre et ont depuis longtemps fasciné l'Homme. Cependant, leur étude était jusqu'à présent limitée aux observations et expériences comportementales. L'avènement de la biologie moléculaire a d'abord rendu possible l'étude moléculaire et génétique du développement d'organismes modèles tels que la mouche Drosophila melanogaster ou le nématode Caenorhabditis elegans, puis dans un second temps de leur comportement. De telles études deviennent désormais possibles pour les insectes sociaux. Nous avons développé une puce à ADN permettant de déterminer simultanément les niveaux d'expression de 1oooo gènes de la fourmi de feu, Solenopsís invicta. Ces gènes ont été séquencés puis assemblés et annotés à l'aide de pipelines que nous avons développés. En se basant sur les informations obtenues, nous avons créé un portail web, Fourmidable. Ce portail vise à centraliser toutes les informations de séquence, d'annotation et d'expression de gènes, ainsi que les protocoles de laboratoire utilisés pour la recherche sur les fourmis. Par la suite, nous avons utilisé les outils développés pour étudier un aspect particulier de S. invicta. Chez les animaux grégaires, une hiérarchie sociale peut déterminer le statut reproducteur des individus. Suite à la mort d'un individu dominant, les individus subordonnés peuvent entrer en compétition en vue d'améliorer leur statut. Un tel phénomène se produit au sein des colonies monogynes de S. invicta, qui contiennent habituellement une unique reine mère, des milliers d'ouvrières et des centaines de reines vierges ailées. Suite à la mort de la reine mère, dominante, un grand nombre de reines vierges, subordonnées, perdent leurs ailes et activent leurs organes reproducteurs au lieu de partir en vol nuptial. Au cours des semaines suivantes, les ouvrières exécutent la plupart de ces reines sur la base de signaux olfactifs produits lors de l'activation des organes reproducteurs. Afin de mieux comprendre les mécanismes moléculaires impliqués, nous avons étudié l'expression de gènes au début de cette compétition. Nous avons identifié 297 gènes différemment exprimés, dont l'annotation indique qu'ils seraient impliqués dans des processus biologiques dont la communication olfactive, le développement des organes reproducteurs et la métabolisation de l'hormone juvénile. Par la suite, nous avons déterminé que la vitesse à laquelle les reines perdent leurs ailes en début de compétition ainsi que les niveaux d'expression de gènes sont corrélés à la probabilité de survie des reines. Nous concluons en discutant des opportunités offertes par la recherche génétique sur les fourmis ainsi que les défis qu'elle devra surmonter.

Convergent genetic architecture underlies social organization in ants.

Complex adaptive polymorphisms are common in nature, but what mechanisms maintain the underlying favorable allelic combinations [1-4]? The convergent evolution of polymorphic social organization in two independent ant species provides a great opportunity to investigate how genomes evolved under parallel selection. Here, we demonstrate that a large, nonrecombining "social chromosome" is associated with social organization in the Alpine silver ant, Formica selysi. This social chromosome shares architectural characteristics with that of the fire ant Solenopsis invicta [2], but the two show no detectable similarity in gene content. The discovery of convergence at two levels-the phenotype and the genetic architecture associated with alternative social forms-points at general genetic mechanisms underlying transitions in social organization. More broadly, our findings are consistent with recent theoretical studies suggesting that suppression of recombination plays a key role in facilitating coordinated shifts in coadapted traits [5, 6].

Current status of a model system: the gene Gp-9 and its association with social organization in fire ants.

The Gp-9 gene in fire ants represents an important model system for studying the evolution of social organization in insects as well as a rich source of information relevant to other major evolutionary topics. An important feature of this system is that polymorphism in social organization is completely associated with allelic variation at Gp-9, such that single-queen colonies (monogyne form) include only inhabitants bearing B-like alleles while multiple-queen colonies (polygyne form) additionally include inhabitants bearing b-like alleles. A recent study of this system by Leal and Ishida (2008) made two major claims, the validity and significance of which we examine here. After reviewing existing literature, analyzing the methods and results of Leal and Ishida (2008), and generating new data from one of their study sites, we conclude that their claim that polygyny can occur in Solenopsis invicta in the U.S.A. in the absence of expression of the b-like allele Gp-9(b) is unfounded. Moreover, we argue that available information on insect OBPs (the family of proteins to which GP-9 belongs), on the evolutionary/population genetics of Gp-9, and on pheromonal/behavioral control of fire ant colony queen number fails to support their view that GP-9 plays no role in the chemosensory-mediated communication that underpins regulation of social organization. Our analyses lead us to conclude that there are no new reasons to question the existing consensus view of the Gp-9 system outlined in Gotzek and Ross (2007).

Genetic regulation of colony social organization in fire ants: an integrative overview

Expression of colony social organization in fire ants appears to be under the control of a single Mendelian factor of large effect. Variation in colony queen number in Solenopsis invicta and its relatives is associated with allelic variation at the gene Gp-9, but not with variation at other unlinked genes; workers regulate queen identity and number on the basis of Gp-9 genotypic compatibility. Nongeneticfactors, such as prior social experience, queen reproductive status, and local environment, have negligible effects on queen number which illustrates the nearly complete penetrance of Gp-9. As predicted, queen number can be manipulated experimentally by altering worker Gp-9 genotype frequencies. The Gp-9 allele lineage associated with polygyny in South American fire? ants has been retained across multiple speciation events, which may signal the action of balancing selection to maintain social polymorphism in these species. Moreover positive selection is implicated in driving the molecular evolution of Gp-9 in association with the origin of polygyny. The identity of the product of Gp-9 as an odorant-binding protein suggests plausible scenarios for its direct involvement in the regulation of queen number via a role in chemical communication. While these and other lines of evidence show that Gp-9 represents a legitimate candidate gene of major effect, studies aimed at determining (i) the biochemical pathways in which GP-9 functions; (ii) the phenotypic effects of molecular variation at Gp-9 and other pathway genes; and (iii) the potential involvement of genes in linkage disequilibrium with Gp-9 are needed to elucidate the genetic architecture underlying social organization in fire ants. Information that reveals the links between molecular variation, individual phenotype, and colony-level behaviors, combined with behavioral models that incorporate details of the chemical communication involved in regulating queen number will yield a novel integrated view of the evolutionary changes underlying a key social adaptation.

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.