Draft genome of the red harvester ant Pogonomyrmex barbatus

Genes and genomic processes underpinning the social lives of ants

The >15000 ant species are all highly social and show great variation in colony organization, complexity and behavior. The mechanisms by which such sociality evolved, as well as those underpinning the elaboration of ant societies since their ∼140 million year old common ancestor, have long been pondered. Here, we review recent insights generated using various genomic approaches. This includes understanding the molecular mechanisms underlying caste differentiation and the diversity of social structures, studying the impact of eusociality on genomic evolutionary rates, and investigating gene expression changes associated with differences in lifespan between castes. Furthermore, functional studies involving RNAi and CRISPR have recently been successfully applied to ants, opening the door to exciting research that promises to revolutionize the understanding of the evolution and diversification of social living.

The interactions of ants with their biotic environment

This special feature results from the symposium 'Ants 2016: ant interactions with their biotic environments' held in Munich in May 2016 and deals with the interactions between ants and other insects, plants, microbes and fungi, studied at micro- and macroevolutionary levels with a wide range of approaches, from field ecology to next-generation sequencing, chemical ecology and molecular genetics. In this paper, we review key aspects of these biotic interactions to provide background information for the papers of this special feature After listing the major types of biotic interactions that ants engage in, we present a brief overview of ant/ant communication, ant/plant interactions, ant/fungus symbioses, and recent insights about ants and their endosymbionts. Using a large molecular clock-dated Formicidae phylogeny, we map the evolutionary origins of different ant clades' interactions with plants, fungi and hemiptera. Ants' biotic interactions provide ideal systems to address fundamental ecological and evolutionary questions about mutualism, coevolution, adaptation and animal communication.

Building a new research framework for social evolution: intralocus caste antagonism

The breeding and non‐breeding ‘castes’ of eusocial insects provide a striking example of role‐specific selection, where each caste maximises fitness through different morphological, behavioural and physiological trait values. Typically, queens are long‐lived egg‐layers, while workers are short‐lived, largely sterile foragers. Remarkably, the two castes are nevertheless produced by the same genome. The existence of inter‐caste genetic correlations is a neglected consequence of this shared genome, potentially hindering the evolution of caste dimorphism: alleles that increase the productivity of queens may decrease the productivity of workers and vice versa, such that each caste is prevented from reaching optimal trait values. A likely consequence of this ‘intralocus caste antagonism’ should be the maintenance of genetic variation for fitness and maladaptation within castes (termed ‘caste load’), analogous to the result of intralocus sexual antagonism. The aim of this review is to create a research framework for understanding caste antagonism, drawing in part upon conceptual similarities with sexual antagonism. By reviewing both the social insect and sexual antagonism literature, we highlight the current empirical evidence for caste antagonism, discuss social systems of interest, how antagonism might be resolved, and challenges for future research. We also introduce the idea that sexual and caste antagonism could interact, creating a three‐way antagonism over gene expression. This includes unpacking the implications of haplodiploidy for the outcome of this complex interaction.

Hymenoptera Genome Database: Using HymenopteraMine to Enhance Genomic Studies of Hymenopteran Insects

The Hymenoptera Genome Database (HGD; http://hymenopteragenome.org ) is a genome informatics resource for insects of the order Hymenoptera, which includes bees, ants and wasps. HGD provides genome browsers with manual annotation tools (JBrowse/Apollo), BLAST, bulk data download, and a data mining warehouse (HymenopteraMine). This chapter focuses on the use of HymenopteraMine to create annotation data sets that can be exported for use in downstream analyses. HymenopteraMine leverages the InterMine platform to combine genome assemblies and official gene sets with data from OrthoDB, RefSeq, FlyBase, Gene Ontology, UniProt, InterPro, KEGG, Reactome, dbSNP, PubMed, and BioGrid, as well as precomputed gene expression information based on publicly available RNAseq. Built-in template queries provide starting points for data exploration, while the QueryBuilder tool supports construction of complex custom queries. The List Analysis and Genomic Regions search tools execute queries based on uploaded lists of identifiers and genome coordinates, respectively. HymenopteraMine facilitates cross-species data mining based on orthology and supports meta-analyses by tracking identifiers across gene sets and genome assemblies.

Context-dependent expression of the foraging gene in field colonies of ants: the interacting roles of age, environment and task

Task allocation among social insect workers is an ideal framework for studying the molecular mechanisms underlying behavioural plasticity because workers of similar genotype adopt different behavioural phenotypes. Elegant laboratory studies have pioneered this effort, but field studies involving the genetic regulation of task allocation are rare. Here, we investigate the expression of the foraging gene in harvester ant workers from five age- and task-related groups in a natural population, and we experimentally test how exposure to light affects foraging expression in brood workers and foragers. Results from our field study show that the regulation of the foraging gene in harvester ants occurs at two time scales: levels of foraging mRNA are associated with ontogenetic changes over weeks in worker age, location and task, and there are significant daily oscillations in foraging expression in foragers. The temporal dissection of foraging expression reveals that gene expression changes in foragers occur across a scale of hours and the level of expression is predicted by activity rhythms: foragers have high levels of foraging mRNA during daylight hours when they are most active outside the nests. In the experimental study, we find complex interactions in foraging expression between task behaviour and light exposure. Oscillations occur in foragers following experimental exposure to 13 L : 11 D (LD) conditions, but not in brood workers under similar conditions. No significant differences were seen in foraging expression over time in either task in 24 h dark (DD) conditions. Interestingly, the expression of foraging in both undisturbed field and experimentally treated foragers is also significantly correlated with the expression of the circadian clock gene, cycle. Our results provide evidence that the regulation of this gene is context-dependent and associated with both ontogenetic and daily behavioural plasticity in field colonies of harvester ants. Our results underscore the importance of assaying temporal patterns in behavioural gene expression and suggest that gene regulation is an integral mechanism associated with behavioural plasticity in harvester ants.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

Computational genome-wide survey of odorant receptors from two solitary bees Dufourea novaeangliae (Hymenoptera: Halictidae) and Habropoda laboriosa (Hymenoptera: Apidae)

Olfactory/odorant receptors (ORs) probably govern eusocial behaviour in honey bees through detection of cuticular hydrocarbons (CHCs) and queen mandibular gland pheromones (QMP). CHCs are involved in nest-mate recognition whereas QMP acts as sex pheromone for drones and as retinue pheromone for female workers. Further studies on the effect of eusociality on the evolution of ORs are hindered by the non-availability of comprehensive OR sets of solitary species. We report complete OR repertoires from two solitary bees Dufourea novaeangliae (112 ORs) and Habropoda laboriosa (151 ORs). We classify these ORs into 34 phylogenetic clades/subfamilies. Differences in the OR sets of solitary and eusocial bees are observed in individual subfamilies like subfamily 9-exon (putative CHC receptors) and L (contains putative QMP receptor group). A subfamily (H) including putative floral scent receptors is expanded in the generalist honey bees only, but not in the specialists. On the contrary, subfamily J is expanded in all bees irrespective of their degree of social complexity or food preferences. Finally, we show species-lineage specific and OR-subfamily specific differences in the putative cis-regulatory DNA motifs of the ORs from six hymenopteran species. Out of these, [A/G]CGCAAGCG[C/T] is a candidate master transcription factor binding site for multiple olfactory genes.

Specialized odorant receptors in social insects that detect cuticular hydrocarbon cues and candidate pheromones

Eusocial insects use cuticular hydrocarbons as components of pheromones that mediate social behaviours, such as caste and nestmate recognition, and regulation of reproduction. In ants such as Harpegnathos saltator, the queen produces a pheromone which suppresses the development of workers' ovaries and if she is removed, workers can transition to a reproductive state known as gamergate. Here we functionally characterize a subfamily of odorant receptors (Ors) with a nine-exon gene structure that have undergone a massive expansion in ants and other eusocial insects. We deorphanize 22 representative members and find they can detect cuticular hydrocarbons from different ant castes, with one (HsOr263) that responds strongly to gamergate extract and a candidate queen pheromone component. After systematic testing with a diverse panel of hydrocarbons, we find that most Harpegnathos saltator Ors are narrowly tuned, suggesting that several receptors must contribute to detection and discrimination of different cuticular hydrocarbons important in mediating eusocial behaviour.Cuticular hydrocarbons (CHC) mediate the interactions between individuals in eusocial insects, but the sensory receptors for CHCs are unclear. Here the authors show that in ants such as H. saltator, the 9-exon subfamily of odorant receptors (HsOrs) responds to CHCs, and ectopic expression of HsOrs in Drosophila neurons imparts responsiveness to CHCs.

orco Mutagenesis Causes Loss of Antennal Lobe Glomeruli and Impaired Social Behavior in Ants

Life inside ant colonies is orchestrated with diverse pheromones, but it is not clear how ants perceive these social signals. It has been proposed that pheromone perception in ants evolved via expansions in the numbers of odorant receptors (ORs) and antennal lobe glomeruli. Here, we generate the first mutant lines in the clonal raider ant, Ooceraea biroi, by disrupting orco, a gene required for the function of all ORs. We find that orco mutants exhibit severe deficiencies in social behavior and fitness, suggesting they are unable to perceive pheromones. Surprisingly, unlike in Drosophila melanogaster, orco mutant ants also lack most of the ∼500 antennal lobe glomeruli found in wild-type ants. These results illustrate that ORs are essential for ant social organization and raise the possibility that, similar to mammals, receptor function is required for the development and/or maintenance of the highly complex olfactory processing areas in the ant brain. VIDEO ABSTRACT.

Genomic exaptation enables Lasius niger adaptation to urban environments

Background

The world is rapidly urbanizing, and only a subset of species are able to succeed in stressful city environments. Efficient genome-enabled stress response appears to be a likely prerequisite for urban adaptation. Despite the important role ants play in the ecosytem, only the genomes of ~13 have been sequenced so far. Here, we present the draft genome assembly of the black garden ant Lasius niger – the most successful urban inhabitant of all ants – and we compare it with the genomes of other ant species, including the closely related Camponotus floridanus.

Results

Sequences from 272 M Illumina reads were assembled into 41,406 contigs with total length of 245 MB, and N50 of 16,382 bp, similar to other ant genome assemblies enabling comparative genomic analysis. Remarkably, the predicted proteome of L. niger is significantly enriched relative to other ant genomes in terms of abundance of domains involved in nucleic acid binding, DNA repair, and nucleotidyl transferase activity, reflecting transposable element proliferation and a likely genomic response. With respect to environmental stress, we note a proliferation of various detoxification genes, including glutatione-S-transferases and those in the cytochrome P450 families. Notably, the CYP9 family is highly expanded with 19 complete and 21 nearly complete members - over twice as many compared to other ants. This family exhibits the signatures of strong directional selection, with eleven positively selected positions in ligand-binding pockets of enzymes. Gene family contraction was detected for several components of the olfactory system, accompanied by instances of both directional selection and relaxation.

Conclusions

Our results suggest that the success of L. niger in urbanized areas may be the result of fortuitous coincidence of several factors, including the expansion of the CYP9 cytochrome family due to coevolution with parasitic fungi, the diversification of DNA repair systems as an answer to proliferation of retroelements, and the reduction of olfactory system and behavioral preadaptations from non-territorial subdominant life strategies found in natural environments. Diversification of cytochromes and DNA repair systems along with reduced odorant communication are the basis of L. niger pollutant resistance and polyphagy, while non-territorial and mobilization strategies allows more efficient exploitation of large but patchy food sources.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0867-x) contains supplementary material, which is available to authorized users.

Identification and characterization of the distinct expression profiles of candidate chemosensory membrane proteins in the antennal transcriptome of Adelphocoris lineolatus (Goeze)

Chemosensory membrane proteins, including odorant receptors (ORs), ionotropic receptors (IRs), gustatory receptors (GRs) and sensory neurone membrane proteins (SNMPs), are supposed to be crucial macromolecules in the insect olfactory signal transduction pathway. The alfalfa plant bug Adelphocoris lineolatus (Goeze) (Hemiptera: Miridae) is highly attracted to high-nitrogen or flowering plants and destroys many important agricultural crops. We assembled the antennal transcriptome of A. lineolatus using Illumina sequencing technology and identified a total of 108 transcripts encoding chemosensory membrane proteins (88 ORs, 12 IRs, four GRs and four SNMPs), amongst which 90 candidates appeared to be full length. Subsequently, both semiquantitative reverse transcription PCR and quantitative real-time PCR experiments were performed to investigate their tissue- and sex-biased expression profiles. The results showed that nearly all of the 108 candidate chemosensory membrane protein genes were largely expressed in adult antennae, and some genes additionally displayed significant differences in the expression levels between sexes. The results of our phylogenetic analysis and the detailed tissue- and sex-biased expression characteristics given here provide an important foundation for further understanding of the complex chemoreception system of the alfalfa plant bug and other Hemiptera species, which also could help us use chemosensory membrane proteins as targets to manipulate insect olfactory behaviour and broaden the applications of available tools for insect pest control.

Transcriptomics and neuroanatomy of the clonal raider ant implicate an expanded clade of odorant receptors in chemical communication

A major aim of sociogenomic research is to uncover common principles in the molecular evolution of sociality. This endeavor has been hampered by the small number of specific genes currently known to function in social behavior. Here we provide several lines of evidence suggesting that ants have evolved a large and novel clade of odorant receptor (OR) genes to perceive hydrocarbon-based pheromones, arguably the most important signals in ant communication. This genomic expansion is also mirrored in the ant brain via a corresponding expansion of a specific cluster of glomeruli in the antennal lobe. We show that in the clonal raider ant, hydrocarbon-sensitive basiconic sensilla are found only on the ventral surface of the female antennal club. Correspondingly, nearly all genes in a clade of 180 ORs within the 9-exon subfamily of ORs are expressed exclusively in females and are highly enriched in expression in the ventral half of the antennal club. Furthermore, we found that across species and sexes, the number of 9-exon ORs expressed in antennae is tightly correlated with the number of glomeruli in the antennal lobe region innervated by odorant receptor neurons from basiconic sensilla. Evolutionary analyses show that this clade underwent a striking gene expansion in the ancestors of all ants and slower but continued expansion in extant ant lineages. This evidence suggests that ants have evolved a large clade of genes to support pheromone perception and that gene duplications have played an important role in the molecular evolution of ant communication.

Identification of Complete Repertoire of Apis florea Odorant Receptors Reveals Complex Orthologous Relationships with Apis mellifera

We developed a computational pipeline for homology based identification of the complete repertoire of olfactory receptor (OR) genes in the Asian honey bee species, Apis florea. Apis florea is phylogenetically the most basal honey bee species and also the most distant sister species to the Western honey bee Apis mellifera, for which all OR genes had been identified before. Using our pipeline, we identified 180 OR genes in A. florea, which is very similar to the number of ORs identified in A. mellifera (177 ORs). Many characteristics of the ORs including gene structure, synteny of tandemly repeated ORs and basic phylogenetic clustering are highly conserved. The composite phylogenetic tree of A. florea and A. mellifera ORs could be divided into 21 clades which are in harmony with the existing Hymenopteran tree. However, we found a few nonorthologous OR relationships between both species as well as independent pseudogenization of ORs suggesting separate evolutionary changes. Particularly, a subgroup of the OR gene clade XI, which had been hypothesized to code cuticular hydrocarbon receptors showed a high number of species-specific ORs. RNAseq analysis detected a total number of 145 OR transcripts in male and 162 in female antennae. Most of the OR genes were highly expressed on the female antennae. However, we detected five distinct male-biased OR genes, out of which three genes (AfOr11, AfOr18, AfOr170P) were shown to be male-biased in A. mellifera, too, thus corroborating a behavioral function in sex-pheromone communication.

Comparative genomics reveals convergent rates of evolution in ant-plant mutualisms

Symbiosis-the close and often long-term interaction of species-is predicted to drive genome evolution in a variety of ways. For example, parasitic interactions have been shown to increase rates of molecular evolution, a trend generally attributed to the Red Queen Hypothesis. However, it is much less clear how mutualisms impact the genome, as both increased and reduced rates of change have been predicted. Here we sequence the genomes of seven species of ants, three that have convergently evolved obligate plant-ant mutualism and four closely related species of non-mutualists. Comparing these sequences, we investigate how genome evolution is shaped by mutualistic behaviour. We find that rates of molecular evolution are higher in the mutualists genome wide, a characteristic apparently not the result of demography. Our results suggest that the intimate relationships of obligate mutualists may lead to selective pressures similar to those seen in parasites, thereby increasing rates of evolution.

Phylogenomics and Divergence Dating of Fungus-Farming Ants (Hymenoptera: Formicidae) of the Genera Sericomyrmex and Apterostigma

Fungus-farming ("attine") ants are model systems for studies of symbiosis, coevolution, and advanced eusociality. A New World clade of nearly 300 species in 15 genera, all attine ants cultivate fungal symbionts for food. In order to better understand the evolution of ant agriculture, we sequenced, assembled, and analyzed transcriptomes of four different attine ant species in two genera: three species in the higher-attine genus Sericomyrmex and a single lower-attine ant species, Apterostigma megacephala, representing the first genomic data for either genus. These data were combined with published genomes of nine other ant species and the honey bee Apis mellifera for phylogenomic and divergence-dating analyses. The resulting phylogeny confirms relationships inferred in previous studies of fungus-farming ants. Divergence-dating analyses recovered slightly older dates than most prior analyses, estimating that attine ants originated 53.6-66.7 million of years ago, and recovered a very long branch subtending a very recent, rapid radiation of the genus Sericomyrmex. This result is further confirmed by a separate analysis of the three Sericomyrmex species, which reveals that 92.71% of orthologs have 99% - 100% pairwise-identical nucleotide sequences. We searched the transcriptomes for genes of interest, most importantly argininosuccinate synthase and argininosuccinate lyase, which are functional in other ants but which are known to have been lost in seven previously studied attine ant species. Loss of the ability to produce the amino acid arginine has been hypothesized to contribute to the obligate dependence of attine ants upon their cultivated fungi, but the point in fungus-farming ant evolution at which these losses occurred has remained unknown. We did not find these genes in any of the sequenced transcriptomes. Although expected for Sericomyrmex species, the absence of arginine anabolic genes in the lower-attine ant Apterostigma megacephala strongly suggests that the loss coincided with the origin of attine ants.

Proteomic analysis of the venom and venom sac of the woodwasp, Sirex noctilio - Towards understanding its biological impact

The European horntail woodwasp, Sirex noctilio, is an invasive insect that attacks conifer hosts, particularly Pinus species. Venom injected by female S. noctilio, together with its symbiotic fungus, damages the normal physiology of Pinus, leading to death of the tree. To identify the proteinaceous components in the venom and uncover the interplay between venom proteins and tree proteins, clarification of the overall profile of proteins produced in the venom gland apparatus was carried out in this work. The venom sac proteome utilised in-solution digested in either a natural or deglycosylated state, prior to nanoHPLC LTQ-Orbitrap under CID/ETD mode. Here, we report the identification of 1454 and 1225 proteins in venom and sac, respectively, with 410 mutual proteins. Approximately 90 proteins were predicted to be secretory, of which 8 have features characteristic of toxins. Chemosensory binding proteins were also identified. Gene ontology and KEGG pathway analysis were employed to predict the protein functions and biological pathways in venom and sac. Protein-protein interaction (PPI) analysis suggested that one-step responses represent the majority of the Sirex-Pinus PPIs, and the proteins representing network hub nodes could be of importance for the development of pest management strategies.

SIGNIFICANCE

The woodwasp Sirex noctilio is an invasive species in many parts of the world, including Australia and North America, where it is considered within the top 10 most serious forest insects. Where they have been introduced, the female woodwasps attack living pine trees, causing significant economic losses. Central to this destruction is the woodwasp's life cycle requirement to bore a hole to deposit eggs and a toxic mucus that disables the tree's network for transporting water and nutrients, yet aids in larval survival. Here we specifically examine the mucus gland apparatus and its contents, revealing the protein components that together with 'noctilisin' facilitate this complex association. The identification of chemosensory binding proteins further supports a role for the woodwasp ovipositor as an instrument for early stages of host tree selection. These findings could provide important clues towards the development of novel control tools against this pest.

Gene expression and variation in social aggression by queens of the harvester ant Pogonomyrmex californicus

A key requirement for social cooperation is the mitigation and/or social regulation of aggression towards other group members. Populations of the harvester ant Pogonomyrmex californicus show the alternate social phenotypes of queens founding nests alone (haplometrosis) or in groups of unrelated yet cooperative individuals (pleometrosis). Pleometrotic queens display an associated reduction in aggression. To understand the proximate drivers behind this variation, we placed foundresses of the two populations into social environments with queens from the same or the alternate population, and measured their behaviour and head gene expression profiles. A proportion of queens from both populations behaved aggressively, but haplometrotic queens were significantly more likely to perform aggressive acts, and conflict escalated more frequently in pairs of haplometrotic queens. Whole-head RNA sequencing revealed variation in gene expression patterns, with the two populations showing moderate differentiation in overall transcriptional profile, suggesting that genetic differences underlie the two founding strategies. The largest detected difference, however, was associated with aggression, regardless of queen founding type. Several modules of coregulated genes, involved in metabolism, immune system and neuronal function, were found to be upregulated in highly aggressive queens. Conversely, nonaggressive queens exhibited a striking pattern of upregulation in chemosensory genes. Our results highlight that the social phenotypes of cooperative vs. solitary nest founding tap into a set of gene regulatory networks that seem to govern aggression level. We also present a number of highly connected hub genes associated with aggression, providing opportunity to further study the genetic underpinnings of social conflict and tolerance.

ArthropodaCyc: a CycADS powered collection of BioCyc databases to analyse and compare metabolism of arthropods

Arthropods interact with humans at different levels with highly beneficial roles (e.g. as pollinators), as well as with a negative impact for example as vectors of human or animal diseases, or as agricultural pests. Several arthropod genomes are available at present and many others will be sequenced in the near future in the context of the i5K initiative, offering opportunities for reconstructing, modelling and comparing their metabolic networks. In-depth analysis of these genomic data through metabolism reconstruction is expected to contribute to a better understanding of the biology of arthropods, thereby allowing the development of new strategies to control harmful species. In this context, we present here ArthropodaCyc, a dedicated BioCyc collection of databases using the Cyc annotation database system (CycADS), allowing researchers to perform reliable metabolism comparisons of fully sequenced arthropods genomes. Since the annotation quality is a key factor when performing such global genome comparisons, all proteins from the genomes included in the ArthropodaCyc database were re-annotated using several annotation tools and orthology information. All functional/domain annotation results and their sources were integrated in the databases for user access. Currently, ArthropodaCyc offers a centralized repository of metabolic pathways, protein sequence domains, Gene Ontology annotations as well as evolutionary information for 28 arthropod species. Such database collection allows metabolism analysis both with integrated tools and through extraction of data in formats suitable for systems biology studies.

Database URL:http://arthropodacyc.cycadsys.org/

The Genome and Methylome of a Subsocial Small Carpenter Bee, Ceratina calcarata

Understanding the evolution of animal societies, considered to be a major transition in evolution, is a key topic in evolutionary biology. Recently, new gateways for understanding social evolution have opened up due to advances in genomics, allowing for unprecedented opportunities in studying social behavior on a molecular level. In particular, highly eusocial insect species (caste-containing societies with nonreproductives that care for siblings) have taken center stage in studies of the molecular evolution of sociality. Despite advances in genomic studies of both solitary and eusocial insects, we still lack genomic resources for early insect societies. To study the genetic basis of social traits requires comparison of genomes from a diversity of organisms ranging from solitary to complex social forms. Here we present the genome of a subsocial bee, Ceratina calcarata This study begins to address the types of genomic changes associated with the earliest origins of simple sociality using the small carpenter bee. Genes associated with lipid transport and DNA recombination have undergone positive selection in C. calcarata relative to other bee lineages. Furthermore, we provide the first methylome of a noneusocial bee. Ceratina calcarata contains the complete enzymatic toolkit for DNA methylation. As in the honey bee and many other holometabolous insects, DNA methylation is targeted to exons. The addition of this genome allows for new lines of research into the genetic and epigenetic precursors to complex social behaviors.

Patterns of molecular evolution of RNAi genes in social and socially parasitic bumblebees

The high frequency of interactions amongst closely related individuals in social insect colonies enhances pathogen transmission. Group-mediated behavior supporting immune defenses tends to decrease selection acting on immune genes. Along with low effective population sizes this might result in relaxed constraint and rapid evolution of immune system genes. Here, we show that antiviral siRNA genes show high rates of molecular evolution with argonaute 2, armitage and maelstrom evolving faster in social bumblebees compared to their socially parasitic cuckoo bumblebees that lack a worker caste. RNAi genes show frequent positive selection at the codon level additionally supported by the occurrence of parallel evolution. Their evolutionary rate is linked to their pathway specific position with genes directly interacting with viruses showing the highest rates of molecular evolution. We suggest that higher pathogen load in social insects indeed drives the molecular evolution of immune genes including antiviral siRNA, if not compensated by behavior.

Ant Genetics: Reproductive Physiology, Worker Morphology, and Behavior

Many exciting studies have begun to elucidate the genetics of the morphological and physiological diversity of ants, but as yet few studies have investigated the genetics of ant behavior directly. Ant genomes are marked by extreme rates of gene turnover, especially in gene families related to olfactory communication, such as the synthesis of cuticular hydrocarbons and the perception of environmental semiochemicals. Transcriptomic and epigenetic differences are apparent between reproductive and sterile females, males and females, and workers that differ in body size. Quantitative genetic approaches suggest heritability of task performance, and population genetic studies indicate a genetic association with reproductive status in some species. Gene expression is associated with behavior including foraging, response to queens attempting to join a colony, circadian patterns of task performance, and age-related changes of task. Ant behavioral genetics needs further investigation of the feedback between individual-level physiological changes and socially mediated responses to environmental conditions.

Expansion of a bitter taste receptor family in a polyphagous insect herbivore

The Insect taste system plays a central role in feeding behaviours and co-evolution of insect-host interactions. Gustatory receptors form the interface between the insect taste system and the environment. From genome and transcriptome sequencing we identified 197 novel gustatory receptor (GR) genes from the polyphagous pest Helicoverpa armigera. These GRs include a significantly expanded bitter receptor family (180 GRs) that could be further divided into three categories based on polypeptide lengths, gene structure and amino acid sequence. Type 1 includes 29 bitter Gr genes that possess introns. Type 2 includes 13 long intronless bitter Gr genes, while Type 3 comprises 131 short intronless bitter Gr genes. Calcium imaging analysis demonstrated that three Type 3 GRs (HarmGR35, HarmGR50 and HarmGR195) can be activated by a crude extract of cotton leaves. HarmGR195, a GR specifically and selectively expressed in adult tarsi, showed a specific response to proline, an amino acid widely present in plant tissues. We hypothesise that the expansion in the H. armigera GR family may be functionally tied to its polyphagous behavior. Understanding the molecular basis of polyphagy may provide opportunities for the development of new environmentally friendly pest control strategies.

Genome-wide analysis of ionotropic receptors provides insight into their evolution in Heliconius butterflies

BACKGROUND

In a world of chemical cues, smell and taste are essential senses for survival. Here we focused on Heliconius, a diverse group of butterflies that exhibit variation in pre- and post-zygotic isolation and chemically-mediated behaviors across their phylogeny. Our study examined the ionotropic receptors, a recently discovered class of receptors that are some of the most ancient chemical receptors.

RESULTS

We found more ionotropic receptors in Heliconius (31) than in Bombyx mori (25) or in Danaus plexippus (27). Sixteen genes in Lepidoptera were not present in Diptera. Only IR7d4 was exclusively found in butterflies and two expansions of IR60a were exclusive to Heliconius. A genome-wide comparison between 11 Heliconius species revealed instances of pseudogenization, gene gain, and signatures of positive selection across the phylogeny. IR60a2b and IR60a2d are unique to the H. melpomene, H. cydno, and H. timareta clade, a group where chemosensing is likely involved in pre-zygotic isolation. IR60a2b also displayed copy number variations (CNVs) in distinct populations of H. melpomene and was the only gene significantly higher expressed in legs and mouthparts than in antennae, which suggests a gustatory function. dN/dS analysis suggests more frequent positive selection in some intronless IR genes and in particular in the sara/sapho and melpomene/cydno/timareta clades. IR60a1 was the only gene with an elevated dN/dS along a major phylogenetic branch associated with pupal mating. Only IR93a was differentially expressed between sexes.

CONCLUSIONS

All together these data make Heliconius butterflies one of the very few insects outside Drosophila where IRs have been characterized in detail. Our work outlines a dynamic pattern of IR gene evolution throughout the Heliconius radiation which could be the result of selective pressure to find potential mates or host-plants.

Genome, transcriptome and methylome sequencing of a primitively eusocial wasp reveal a greatly reduced DNA methylation system in a social insect

Comparative genomics of social insects has been intensely pursued in recent years with the goal of providing insights into the evolution of social behaviour and its underlying genomic and epigenomic basis. However, the comparative approach has been hampered by a paucity of data on some of the most informative social forms (e.g. incipiently and primitively social) and taxa (especially members of the wasp family Vespidae) for studying social evolution. Here, we provide a draft genome of the primitively eusocial model insect Polistes dominula, accompanied by analysis of caste-related transcriptome and methylome sequence data for adult queens and workers. Polistes dominula possesses a fairly typical hymenopteran genome, but shows very low genomewide GC content and some evidence of reduced genome size. We found numerous caste-related differences in gene expression, with evidence that both conserved and novel genes are related to caste differences. Most strikingly, these -omics data reveal a major reduction in one of the major epigenetic mechanisms that has been previously suggested to be important for caste differences in social insects: DNA methylation. Along with a conspicuous loss of a key gene associated with environmentally responsive DNA methylation (the de novo DNA methyltransferase Dnmt3), these wasps have greatly reduced genomewide methylation to almost zero. In addition to providing a valuable resource for comparative analysis of social insect evolution, our integrative -omics data for this important behavioural and evolutionary model system call into question the general importance of DNA methylation in caste differences and evolution in social insects.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Genome Survey Sequencing for the Characterization of the Genetic Background of Rosa roxburghii Tratt and Leaf Ascorbate Metabolism Genes

Rosa roxburghii Tratt is an important commercial horticultural crop in China that is recognized for its nutritional and medicinal values. In spite of the economic significance, genomic information on this rose species is currently unavailable. In the present research, a genome survey of R. roxburghii was carried out using next-generation sequencing (NGS) technologies. Total 30.29 Gb sequence data was obtained by HiSeq 2500 sequencing and an estimated genome size of R. roxburghii was 480.97 Mb, in which the guanine plus cytosine (GC) content was calculated to be 38.63%. All of these reads were technically assembled and a total of 627,554 contigs with a N50 length of 1.484 kb and furthermore 335,902 scaffolds with a total length of 409.36 Mb were obtained. Transposable elements (TE) sequence of 90.84 Mb which comprised 29.20% of the genome, and 167,859 simple sequence repeats (SSRs) were identified from the scaffolds. Among these, the mono-(66.30%), di-(25.67%), and tri-(6.64%) nucleotide repeats contributed to nearly 99% of the SSRs, and sequence motifs AG/CT (28.81%) and GAA/TTC (14.76%) were the most abundant among the dinucleotide and trinucleotide repeat motifs, respectively. Genome analysis predicted a total of 22,721 genes which have an average length of 2311.52 bp, an average exon length of 228.15 bp, and average intron length of 401.18 bp. Eleven genes putatively involved in ascorbate metabolism were identified and its expression in R. roxburghii leaves was validated by quantitative real-time PCR (qRT-PCR). This is the first report of genome-wide characterization of this rose species.

Are feeding preferences and insecticide resistance associated with the size of detoxifying enzyme families in insect herbivores?

The size of gene families associated with xenobiotic detoxification in insects may be associated with the complexity of their diets and their propensities to develop insecticide resistance. We test these hypotheses by collating the annotations of cytochrome P450, carboxyl/cholinesterase and glutathione S-transferase genes in 65 insect species with data on their host use and history of insecticide resistance. We find 2-4 fold variation across the species in the numbers of these genes and, in some orders, especially the Hymenoptera, there is a clear relationship between the numbers of genes and feeding preferences. However in other orders, in particular the Lepidoptera, no such relationship is apparent. The size of these three gene families also tend to correlate with insecticide resistance propensity but this may not be an independent effect because species with broader host ranges are more likely to be pests that are heavily sprayed with insecticides.

The evolution of heat shock protein sequences, cis-regulatory elements, and expression profiles in the eusocial Hymenoptera

Background

The eusocial Hymenoptera have radiated across a wide range of thermal environments, exposing them to significant physiological stressors. We reconstructed the evolutionary history of three families of Heat Shock Proteins (Hsp90, Hsp70, Hsp40), the primary molecular chaperones protecting against thermal damage, across 12 Hymenopteran species and four other insect orders. We also predicted and tested for thermal inducibility of eight Hsps from the presence of cis-regulatory heat shock elements (HSEs). We tested whether Hsp induction patterns in ants were associated with different thermal environments.

Results

We found evidence for duplications, losses, and cis-regulatory changes in two of the three gene families. One member of the Hsp90 gene family, hsp83, duplicated basally in the Hymenoptera, with shifts in HSE motifs in the novel copy. Both copies were retained in bees, but ants retained only the novel HSE copy. For Hsp70, Hymenoptera lack the primary heat-inducible orthologue from Drosophila melanogaster and instead induce the cognate form, hsc70-4, which also underwent an early duplication. Episodic diversifying selection was detected along the branch predating the duplication of hsc70-4 and continued along one of the paralogue branches after duplication. Four out of eight Hsp genes were heat-inducible and matched the predictions based on presence of conserved HSEs. For the inducible homologues, the more thermally tolerant species, Pogonomyrmex barbatus, had greater Hsp basal expression and induction in response to heat stress than did the less thermally tolerant species, Aphaenogaster picea. Furthermore, there was no trade-off between basal expression and induction.

Conclusions

Our results highlight the unique evolutionary history of Hsps in eusocial Hymenoptera, which has been shaped by gains, losses, and changes in cis-regulation. Ants, and most likely other Hymenoptera, utilize lineage-specific heat inducible Hsps, whose expression patterns are associated with adaptive variation in thermal tolerance between two ant species. Collectively, our analyses suggest that Hsp sequence and expression patterns may reflect the forces of selection acting on thermal tolerance in ants and other social Hymenoptera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0573-0) contains supplementary material, which is available to authorized users.

An Introduction to Genome Annotation

Genome projects have evolved from large international undertakings to tractable endeavors for a single lab. Accurate genome annotation is critical for successful genomic, genetic, and molecular biology experiments. These annotations can be generated using a number of approaches and available software tools. This unit describes methods for genome annotation and a number of software tools commonly used in gene annotation.

Hymenoptera Genome Database: integrating genome annotations in HymenopteraMine

We report an update of the Hymenoptera Genome Database (HGD) (http://HymenopteraGenome.org), a model organism database for insect species of the order Hymenoptera (ants, bees and wasps). HGD maintains genomic data for 9 bee species, 10 ant species and 1 wasp, including the versions of genome and annotation data sets published by the genome sequencing consortiums and those provided by NCBI. A new data-mining warehouse, HymenopteraMine, based on the InterMine data warehousing system, integrates the genome data with data from external sources and facilitates cross-species analyses based on orthology. New genome browsers and annotation tools based on JBrowse/WebApollo provide easy genome navigation, and viewing of high throughput sequence data sets and can be used for collaborative genome annotation. All of the genomes and annotation data sets are combined into a single BLAST server that allows users to select and combine sequence data sets to search.

The Genome and Methylome of a Beetle with Complex Social Behavior, Nicrophorus vespilloides (Coleoptera: Silphidae)

Testing for conserved and novel mechanisms underlying phenotypic evolution requires a diversity of genomes available for comparison spanning multiple independent lineages. For example, complex social behavior in insects has been investigated primarily with eusocial lineages, nearly all of which are Hymenoptera. If conserved genomic influences on sociality do exist, we need data from a wider range of taxa that also vary in their levels of sociality. Here, we present the assembled and annotated genome of the subsocial beetle Nicrophorus vespilloides, a species long used to investigate evolutionary questions of complex social behavior. We used this genome to address two questions. First, do aspects of life history, such as using a carcass to breed, predict overlap in gene models more strongly than phylogeny? We found that the overlap in gene models was similar between N. vespilloides and all other insect groups regardless of life history. Second, like other insects with highly developed social behavior but unlike other beetles, does N. vespilloides have DNA methylation? We found strong evidence for an active DNA methylation system. The distribution of methylation was similar to other insects with exons having the most methylated CpGs. Methylation status appears highly conserved; 85% of the methylated genes in N. vespilloides are also methylated in the hymentopteran Nasonia vitripennis. The addition of this genome adds a coleopteran resource to answer questions about the evolution and mechanistic basis of sociality and to address questions about the potential role of methylation in social behavior.

Diversification of the ant odorant receptor gene family and positive selection on candidate cuticular hydrocarbon receptors

Background

Chemical communication plays important roles in the social behavior of ants making them one of the most successful groups of animals on earth. However, the molecular evolutionary process responsible for their chemosensory adaptation is still elusive. Recent advances in genomic studies have led to the identification of large odorant receptor (Or) gene repertoires from ant genomes providing fruitful materials for molecular evolution analysis. The aim of this study was to test the hypothesis that diversification of this gene family is involved in olfactory adaptation of each species.

Results

We annotated the Or genes from the genome sequences of two leaf-cutter ants, Acromyrmex echinatior and Atta cephalotes (385 and 376 putative functional genes, respectively). These were used, together with Or genes from Camponotus floridanus, Harpegnathos saltator, Pogonomyrmex barbatus, Linepithema humile, Cerapachys biroi, Solenopsis invicta and Apis mellifera, in molecular evolution analysis. Like the Or family in other insects, ant Or genes evolve by the birth-and-death model of gene family evolution. Large gene family expansions involving tandem gene duplications, and gene gains outnumbering losses, are observed. Codon analysis of genes in lineage-specific expansion clades revealed signatures of positive selection on the candidate cuticular hydrocarbon receptor genes (9-exon subfamily) of Cerapachys biroi, Camponotus floridanus, Acromyrmex echinatior and Atta cephalotes. Positively selected amino acid positions are primarily in transmembrane domains 3 and 6, which are hypothesized to contribute to the odor-binding pocket, presumably mediating changing ligand specificity.

Conclusions

This study provides support for the hypothesis that some ant lineage-specific Or genes have evolved under positive selection. Newly duplicated genes particularly in the candidate cuticular hydrocarbon receptor clade that have evolved under positive selection may contribute to the highly sophisticated lineage-specific chemical communication in each ant species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1371-x) contains supplementary material, which is available to authorized users.

Cuticular Hydrocarbon Pheromones for Social Behavior and Their Coding in the Ant Antenna

The sophisticated organization of eusocial insect societies is largely based on the regulation of complex behaviors by hydrocarbon pheromones present on the cuticle. We used electrophysiology to investigate the detection of cuticular hydrocarbons (CHCs) by female-specific olfactory sensilla basiconica on the antenna of Camponotus floridanus ants through the utilization of one of the largest family of odorant receptors characterized so far in insects. These sensilla, each of which contains multiple olfactory receptor neurons, are differentially sensitive to CHCs and allow them to be classified into three broad groups that collectively detect every hydrocarbon tested, including queen and worker-enriched CHCs. This broad-spectrum sensitivity is conserved in a related species, Camponotus laevigatus, allowing these ants to detect CHCs from both nestmates and non-nestmates. Behavioral assays demonstrate that these ants are excellent at discriminating CHCs detected by the antenna, including enantiomers of a candidate queen pheromone that regulates the reproductive division of labor.

Chemoreceptor Evolution in Hymenoptera and Its Implications for the Evolution of Eusociality

Eusocial insects, mostly Hymenoptera, have evolved unique colonial lifestyles that rely on the perception of social context mainly through pheromones, and chemoreceptors are hypothesized to have played important adaptive roles in the evolution of sociality. However, because chemoreceptor repertoires have been characterized in few social insects and their solitary relatives, a comprehensive examination of this hypothesis has not been possible. Here, we annotate ∼3,000 odorant and gustatory receptors in recently sequenced Hymenoptera genomes and systematically compare >4,000 chemoreceptors from 13 hymenopterans, representing one solitary lineage (wasps) and three independently evolved eusocial lineages (ants and two bees). We observe a strong general tendency for chemoreceptors to expand in Hymenoptera, whereas the specifics of gene gains/losses are highly diverse between lineages. We also find more frequent positive selection on chemoreceptors in a facultative eusocial bee and in the common ancestor of ants compared with solitary wasps. Our results suggest that the frequent expansions of chemoreceptors have facilitated the transition to eusociality. Divergent expression patterns of odorant receptors between honeybee and ants further indicate differential roles of chemoreceptors in parallel trajectories of social evolution.

Analysis of phylogenomic datasets reveals conflict, concordance, and gene duplications with examples from animals and plants

BACKGROUND

The use of transcriptomic and genomic datasets for phylogenetic reconstruction has become increasingly common as researchers attempt to resolve recalcitrant nodes with increasing amounts of data. The large size and complexity of these datasets introduce significant phylogenetic noise and conflict into subsequent analyses. The sources of conflict may include hybridization, incomplete lineage sorting, or horizontal gene transfer, and may vary across the phylogeny. For phylogenetic analysis, this noise and conflict has been accommodated in one of several ways: by binning gene regions into subsets to isolate consistent phylogenetic signal; by using gene-tree methods for reconstruction, where conflict is presumed to be explained by incomplete lineage sorting (ILS); or through concatenation, where noise is presumed to be the dominant source of conflict. The results provided herein emphasize that analysis of individual homologous gene regions can greatly improve our understanding of the underlying conflict within these datasets.

RESULTS

Here we examined two published transcriptomic datasets, the angiosperm group Caryophyllales and the aculeate Hymenoptera, for the presence of conflict, concordance, and gene duplications in individual homologs across the phylogeny. We found significant conflict throughout the phylogeny in both datasets and in particular along the backbone. While some nodes in each phylogeny showed patterns of conflict similar to what might be expected with ILS alone, the backbone nodes also exhibited low levels of phylogenetic signal. In addition, certain nodes, especially in the Caryophyllales, had highly elevated levels of strongly supported conflict that cannot be explained by ILS alone.

CONCLUSION

This study demonstrates that phylogenetic signal is highly variable in phylogenomic data sampled across related species and poses challenges when conducting species tree analyses on large genomic and transcriptomic datasets. Further insight into the conflict and processes underlying these complex datasets is necessary to improve and develop adequate models for sequence analysis and downstream applications. To aid this effort, we developed the open source software phyparts ( https://bitbucket.org/blackrim/phyparts ), which calculates unique, conflicting, and concordant bipartitions, maps gene duplications, and outputs summary statistics such as internode certainy (ICA) scores and node-specific counts of gene duplications.

Xenobiotic detoxification pathways in honey bees

Relative to most other insect genomes, the western honey bee Apis mellifera has a deficit of detoxification genes spanning Phase I (functionalization), II (conjugation) and III (excretion) gene families. Although honeybees do not display across-the-board greater sensitivity to pesticides, this deficit may render them vulnerable to synergistic interactions among xenobiotics. Diet quality, in terms of protein and phytochemical content, has a pronounced influence on tolerance of toxic compounds. Detoxification gene inventory reduction may reflect an evolutionary history of consuming relatively chemically benign nectar and pollen, as other apoid pollinators display comparable levels of cytochrome P450 gene reduction. Enzymatic detoxification in the eusocial A. mellifera may be complemented by behaviors comprising a 'social detoxification system,' including forager discrimination, dilution by pollen mixing, and colony food processing via microbial fermentation, that reduces the number or quantity of ingested chemicals requiring detoxification.

How Do Genomes Create Novel Phenotypes? Insights from the Loss of the Worker Caste in Ant Social Parasites

A central goal of biology is to uncover the genetic basis for the origin of new phenotypes. A particularly effective approach is to examine the genomic architecture of species that have secondarily lost a phenotype with respect to their close relatives. In the eusocial Hymenoptera, queens and workers have divergent phenotypes that may be produced via either expression of alternative sets of caste-specific genes and pathways or differences in expression patterns of a shared set of multifunctional genes. To distinguish between these two hypotheses, we investigated how secondary loss of the worker phenotype in workerless ant social parasites impacted genome evolution across two independent origins of social parasitism in the ant genera Pogonomyrmex and Vollenhovia. We sequenced the genomes of three social parasites and their most-closely related eusocial host species and compared gene losses in social parasites with gene expression differences between host queens and workers. Virtually all annotated genes were expressed to some degree in both castes of the host, with most shifting in queen-worker bias across developmental stages. As a result, despite >1 My of divergence from the last common ancestor that had workers, the social parasites showed strikingly little evidence of gene loss, damaging mutations, or shifts in selection regime resulting from loss of the worker caste. This suggests that regulatory changes within a multifunctional genome, rather than sequence differences, have played a predominant role in the evolution of social parasitism, and perhaps also in the many gains and losses of phenotypes in the social insects.

Caste-biases in gene expression are specific to developmental stage in the ant Formica exsecta

Understanding how a single genome creates and maintains distinct phenotypes is a central goal in evolutionary biology. Social insects are a striking example of co-opted genetic backgrounds giving rise to dramatically different phenotypes, such as queen and worker castes. A conserved set of molecular pathways, previously envisioned as a set of 'toolkit' genes, has been hypothesized to underlie queen and worker phenotypes in independently evolved social insect lineages. Here, we investigated the toolkit from a developmental point of view, using RNA-Seq to compare caste-biased gene expression patterns across three life stages (pupae, emerging adult and old adult) and two female castes (queens and workers) in the ant Formica exsecta. We found that the number of genes with caste-biased expression increases dramatically from pupal to old adult stages. This result suggests that phenotypic differences between queens and workers at the pupal stage may derive from a relatively low number of caste-biased genes, compared to higher number of genes required to maintain caste differences at the adult stage. Gene expression patterns were more similar among castes within developmental stages than within castes despite the extensive phenotypic differences between queens and workers. Caste-biased expression was highly variable among life stages at the level of single genes, but more consistent when gene functions (gene ontology terms) were investigated. Finally, we found that a large part of putative toolkit genes were caste-biased at least in some life stages in F. exsecta, and the caste-biases, but not their direction, were more often shared between F. exsecta and other ant species than between F. exsecta and bees. Our results indicate that gene expression should be examined across several developmental stages to fully reveal the genetic basis of polyphenisms.

Scrutinizing the immune defence inventory of Camponotus floridanus applying total transcriptome sequencing

BACKGROUND

Defence mechanisms of organisms are shaped by their lifestyle, environment and pathogen pressure. Carpenter ants are social insects which live in huge colonies comprising genetically closely related individuals in high densities within nests. This lifestyle potentially facilitates the rapid spread of pathogens between individuals. In concert with their innate immune system, social insects may apply external immune defences to manipulate the microbial community among individuals and within nests. Additionally, carpenter ants carry a mutualistic intracellular and obligate endosymbiotic bacterium, possibly maintained and regulated by the innate immune system. Thus, different selective forces could shape internal immune defences of Camponotus floridanus.

RESULTS

The immune gene repertoire of C. floridanus was investigated by re-evaluating its genome sequence combined with a full transcriptome analysis of immune challenged and control animals using Illumina sequencing. The genome was re-annotated by mapping transcriptome reads and masking repeats. A total of 978 protein sequences were characterised further by annotating functional domains, leading to a change in their original annotation regarding function and domain composition in about 8% of all proteins. Based on homology analysis with key components of major immune pathways of insects, the C. floridanus immune-related genes were compared to those of Drosophila melanogaster, Apis mellifera, and other hymenoptera. This analysis revealed that overall the immune system of carpenter ants comprises many components found in these insects. In addition, several C. floridanus specific genes of yet unknown functions but which are strongly induced after immune challenge were discovered. In contrast to solitary insects like Drosophila or the hymenopteran Nasonia vitripennis, the number of genes encoding pattern recognition receptors specific for bacterial peptidoglycan (PGN) and a variety of known antimicrobial peptide (AMP) genes is lower in C. floridanus. The comparative analysis of gene expression post immune-challenge in different developmental stages of C. floridanus suggests a stronger induction of immune gene expression in larvae in comparison to adults.

CONCLUSIONS

The comparison of the immune system of C. floridanus with that of other insects revealed the presence of a broad immune repertoire. However, the relatively low number of PGN recognition proteins and AMPs, the identification of Camponotus specific putative immune genes, and stage specific differences in immune gene regulation reflects Camponotus specific evolution including adaptations to its lifestyle.

Evolutionary history of the Azteca-like mariner transposons and their host ants

Three different complete mariner elements were found in the genome of the ant Tapinoma nigerrimum. One (Tnigmar-Mr) was interrupted by a 900-bp insertion that corresponded to an incomplete member of a fourth mariner element, called Azteca. In this work, we isolate and characterize full-length Tnigmar-Az elements in T. nigerrimum. The purpose of this study is to clarify the evolutionary history of Azteca elements and their hosts as well as the possible existence of horizontal transfer processes. For this, Azteca-like elements were also retrieved from the available sequences of various ant genomes, representing four different ant subfamilies: Dolichoderinae, Formicinae, Myrmicinae, and Ponerinae. The tree topology resulting for the Azteca-like elements bore very little resemblance to that of their respective hosts. The pervasive presence of Azteca-like elements in all ant genomes, together with the observation that extant copies are usually younger than the genomes that host them, could be explained either by lineage sorting or by recent horizontal transfer of active elements. However, the finding of closer genetic relationships between elements than between the ants that host them is consistent with the latter scenario. This is clearly observed in Sinvmar-Az, Tnigmar-Az, Acepmar-Az, and Cflomar-Az elements, suggesting the existence of horizontal transfer processes. On the contrary, some elements displayed more divergence than did the hosts harboring them. This may reflect either further horizontal transfer events or random lineage sorting.