Original Research ARTICLE
A genomic comparison of two termites with different social complexity
- 1Department of Evolutionary Biology and Ecology, Institute of Biology I, University of Freiburg, Freiburg, Germany
- 2Section for Ecology and Evolution, Department of Biology, Centre for Social Evolution, University of Copenhagen, Copenhagen, Denmark
- 3China National Genebank, BGI-Shenzhen, Shenzhen, China
- 4Centre for GeoGenetics, Natural History Museum of Denmark, University of Copenhagen, Copenhagen, Denmark
- 5School of Life Sciences, Arizona State University, Tempe, AZ, USA
The termites evolved eusociality and complex societies before the ants, but have been studied much less. The recent publication of the first two termite genomes provides a unique comparative opportunity, particularly because the sequenced termites represent opposite ends of the social complexity spectrum. Zootermopsis nevadensis has simple colonies with totipotent workers that can develop into all castes (dispersing reproductives, nest-inheriting replacement reproductives, and soldiers). In contrast, the fungus-growing termite Macrotermes natalensis belongs to the higher termites and has very large and complex societies with morphologically distinct castes that are life-time sterile. Here we compare key characteristics of genomic architecture, focusing on genes involved in communication, immune defenses, mating biology and symbiosis that were likely important in termite social evolution. We discuss these in relation to what is known about these genes in the ants and outline hypothesis for further testing.
The termites are “social cockroaches,” a monophyletic clade (Infraorder “Isoptera”) nested within the Blattodea (Inward et al., 2007a; Engel et al., 2009; Krishna et al., 2013). They superficially resemble the ants in having wingless worker foragers, but are fundamentally different in a series of ancestral traits that affect the organization of their eusocial colonies (Korb, 2008; Howard and Thorne, 2011). The (eu)social Hymenoptera are haplodiploid holometabolous insects whose males develop from haploid eggs and have transient roles in social life, because they survive only as sperm stored in the spermatheca of queens. Hymenopteran colonies thus consist of female adults that develop from fertilized eggs to differentiate into workers, virgin queens and occasionally soldiers of which only the former care for the helpless grub-like larvae. By contrast, termites are diploid hemimetabolous insects whose colonies usually have workers, soldiers, and reproductives of both sexes. Both have life-time monogamy upon colony founding as ancestral state (Hughes et al., 2008; Boomsma, 2013), but in contrast to the eusocial Hymenoptera, royal pairs regularly remate to produce immatures that increasingly come to resemble the workers, soldiers, and reproductives into which they differentiate. Hence, termite caste differentiation is based on phenotypic plasticity among immatures (Korb and Hartfelder, 2008; Miura and Scharf, 2011), while the eusocial Hymenoptera have castes of adults (Wilson, 1971).
Termites and ants also share many traits that convergently evolved in response to similar selective pressures (Thorne and Traniello, 2003; Korb, 2008; Howard and Thorne, 2011). Both are mostly soil-dwelling and thus continuously exposed to high pathogen loads and their long-lived, populous and genetically homogenous colonies appear to be ideal targets for infections (Schmid-Hempel, 1998). However, both the ants and the termites also evolved impressive disease defense strategies, which have implied that very few pathogens have been able to specialize on infecting perennial ant and termite colonies over evolutionary time (Boomsma et al., 2005). In large part this appears to be due to immune defenses operating both at the individual and the collective (social immunity) level (Cremer et al., 2007; Rosengaus et al., 2011). Another common characteristic of the ants and termites is that both evolved complex communication systems that largely rely on chemical cues, such as cuticular hydrocarbons (CHCs), for nestmate recognition and within-colony communication (e.g., Liebig, 2010; Van Zweden and D'Ettorre, 2010). Strikingly, long-chained CHCs of queens often appear to function as fertility signals for workers of both lineages (Liebig et al., 2009; Weil et al., 2009; Liebig, 2010; van Oystaeyen et al., 2014). Here, we offer the first comparative exploration of the extent to which lineage ancestry has determined these convergent phenotypic similarities based on the first two termite genomes that became recently available (Poulsen et al., 2014; Terrapon et al., 2014).
The two termite genomes represent opposite ends of the social complexity spectrum within the Isoptera (Roisin, 2000) (Table 1) as they exemplify the two fundamental termite life types: the wood-dwelling one-piece nesters and the central place foraging lineages that generally differ in social complexity, feeding ecology, gut symbionts, and developmental plasticity (Abe, 1987; Korb, 2007; Korb and Hartfelder, 2008) (Figure 1). Zootermopsis nevadensis belongs to the former type and Macrotermes natalensis to the latter. Wood-dwelling species (Abe, 1987; Shellman-Reeve, 1997) nest within a single piece of dead wood that serves both as food and nesting habitat so the termites never leave their nest to forage. This social syndrome is widely considered to be ancestral (e.g., Noirot and Pasteels, 1987, 1988; Inward et al., 2007b) and associated with high degrees of developmental plasticity for the individual termites (Figure 2A). Workers remain totipotent immatures throughout several instars that commonly develop further into sterile soldiers, winged sexuals (alates) that found new nests as primary reproductives, or neotenic reproductives that reproduce within the natal nest (Figure 2A).
Figure 1. Simplified phylogeny of the main termite study species with their key traits. Shown is a cladogram of termite genera on which some genomic/molecular genetic research has been done. Added to the right are characteristic social and ecological traits. Social: increasing social complexity from + to +++ (e.g., increasing colony size, division of labor, morphological differentiation between castes); Type: life type, foraging vs. wood dwelling; Region: temperate vs. tropical; Pathogens: soil pathogens vs. wood-decaying fungi, +, present; −, absent. Study species (photo credits): Nasutitermes takasagoensis (Kenji Matsuura), Macrotermes natalensis (Judith Korb), Reticulitermes speratus (Kenji Matsuura), Reticulitermes flavipes (not shown), Coptotermes formosanus (not shown), Prorhinotermes simplex (Judith Korb), Cryptotermes secundus (Judith Korb), Zootermopsis nevadensis (Judith Korb), Hodotermes sjostedti (Toru Miura).
Figure 2. Developmental pathways of (A) wood-dwelling termites such as Z. nevadensis and (B) foraging higher termites such as M. natalensis. Wood-dwelling termites have totipotent immature stages that can explore all caste options, whereas higher termites have a bifurcating caste development pathway splitting into a nymphal line leading to winged dispersing alates and an apterous line leading to workers and soldiers. In M. natalensis this bifurcation is already established in the egg stage. (i) progressive development via nymphal instar(s) into winged sexuals (alates) that disperse and found a new nest elsewhere; (ii) stationary molt remaining in the same instar; (iii) regressive development into an “earlier” instar (gray semi-circle); (iv) development into a soldier, and (v) development into a neotenic replacement reproductive that reproduces within the natal nest. Part (a) is adapted from Korb et al. (2012b). (Photo credits: Judith Korb).
The foraging termite species (also called “multiple piece nesters”; Abe, 1987; Shellman-Reeve, 1997) forage for food outside the nest at some point after colony foundation and bring it back to the colony to feed nestmates. They represent more than 85% of the extant termite species (Kambhampati and Eggleton, 2000). They have true workers and an early separation into distinct developmental pathways (Roisin, 2000; Korb and Hartfelder, 2008) (Figure 2B). In the apterous line, individuals are unable to develop wings and can thus never disperse as reproductives. They become workers and soldiers, but can in some species also advance to become neotenic reproductives in their own nest. In the nymphal line, however, individuals develop wings and dispersing phenotypes that found new colonies elsewhere (Figure 2B). The Macrotermitinae to which Macrotermes natalensis belongs are special examples of foraging termites because their colonies are dependent on nutrition provided by a Termitomyces symbiont (Basidiomycota: Agaricales) (Wood and Thomas, 1989; Nobre et al., 2011). This fungal symbiosis is evolutionarily derived and comes in addition to more fundamental protist (lower termites) and bacterial gut symbionts (all termites), which have played major roles throughout termite evolution. Macrotermes species have two (major/minor) worker castes and two (major/minor) soldier castes (Ruelle, 1970) that may be determined as early as the egg stage (suggested for Macrotermes michaelseni by Okot-Kotber, 1985). Macrotermes colonies often build conspicuous mounds that may harbor several millions of individuals (Noirot and Darlington, 2000; Korb, 2011).
We compare the genomes of these divergent species (Table 1) with those of other insects and outline first hypotheses how sociality and ecological factors left their footprints in the genomes.
Materials and Methods
Construction of Gene Families
To gain insight into the evolution of gene families in termites, we clustered genes from 12 insect genomes (pea aphid: Acyrthosiphon pisum: The International Pea Aphid Genomics Consortium, 2010; body louse: Pediculus humanus: Kirkness et al., 2010; flour beetle: Tribolium castaneum: Richards et al., 2008; fruitfly: Drosophila melanogaster: Adams et al., 2000; jewel wasp: Nasonia vitripennis: Werren et al., 2010; honeybee: Apis mellifera: The Honeybee Genome Sequencing Consortium, 2006; ants: Acromyrmex echinatior: Nygaard et al., 2011, Atta cephalotes: Suen et al., 2011, Camponotus floridanus, Harpegnathos saltator: Bonasio et al., 2010; termites: Z. nevadensis, M. natalensis), the water flea Daphnia pulex (Colbourne et al., 2011), and the round worm Caenorhabditis elegans (Coulson and C. elegans Genome Consortium, 1996). The gene sets of the species that we chose were downloaded from the Ensembl database (Flicek et al., 2014), except for ants and termites which were downloaded from their own reference databases. Then we used Treefam (Li et al., 2006) to construct gene families. For more information see also Terrapon et al. (2014) and Poulsen et al. (2014) (Table S1).
Functional Annotation of Termite Genes
InterproScan v4.8 (Zdobnov and Apweiler, 2001) was used to annotate motifs and domains of translated proteins in two termites. Protein sequences were searched against SUPERFAMILY, Pfam, PRINTS, PROSITE, ProDom, Gene3D, PANTHER, and SMART databases in Interpro with default parameter settings. GO (gene ontology) terms for each gene were obtained from the Interpro database according to the relationship of GO and Interpro terms. The KEGG annotation (Kanehisa and Goto, 2000) was done via the KAAS online server (Moriya et al., 2007) using the SBH method against the eukaryotic species set.
Some gene families were termite-specific and absent from the other investigated genomes. For these genes we performed functional enrichment analyses of GO and IPR (Interpro domain) annotation. P-values for significant difference were obtained by χ2-tests adjusted by FDR (false discovery rate). Similarly, we analyzed differences between the gene sets of Z. nevadensis und M. natalensis by comparing IPR annotation, KEGG pathways, and gene families. We constructed gene families for both genomes using Treefam (Li et al., 2006) and tested for differences in gene numbers using χ2-tests (or Fisher's exact test for small sample sizes). For gene families that were specific to M. natalensis and/or Z. nevadensis, we performed IPR enrichment analyses to obtain information on the putative functions of these genes.
We used the M. natalensis and Z. nevadensis genome assemblies to perform repetitive sequence annotation. First, we did homologous repeat family annotation to identify transposable elements (TEs) using the TE database Repbase v17.06 (Jurka and Kapitonov, 2005) and the programs RepeatMasker (parameter –norna) and RepeatProteinMask v4.0.1 (http://www.RepeatMasker.org) (parameter –p 0.0001) (Smit et al., 1996-2010). De-novo repeat family annotation was done with PILER v1.0 (Edgar and Myers, 2005), LTRfinder v1.05 (Zhao and Wang, 2007) and RepeatModeler v1.05, (http://www.RepeatMasker.org) (Smit et al., 1996-2010) using default parameters. TEs identified by PILER were converted into TE families and aligned with Muscle v3.28 (Edgar, 2004) to obtain consensus sequences from the alignments. In order to reduce redundancy in the results of LTRfinder and PILER, an “all against all” BLASTn (e-value 1e-5) was performed. If sequences overlapped for more than 80% we kept the longer TE.
We combined the TE families with the consensus sequences of LTRfinder and PILER together with those identified using RepeatModeler to obtain the final TE sequence library for the two termites. All TE sequences were classified with RepeatClassifier in the RepeatModeler package against Repbase v17.06 (Jurka and Kapitonov, 2005) (Dataset S1). Finally, we used the de novo TE library to annotate all TEs in the two genomes and combined the results of homologous TE annotation and the de novo annotation. If there were overlapping annotations we kept the longer TE. In addition, we predicted tandem repeats using TRF finder (parameters settings: match = 2, mismatch = 7, delta = 7, PM = 80, PI = 10, Minscore = 50, and MaxPeriod = 12) (Benson, 1999). In total, the non-redundant repetitive sequences accounted for 27.8 and 45.9% of the Z. nevadensis and M. natalensis genome, respectively (Table 2, Dataset S1).
We also checked for Talua elements in both termite species, SINE elements that were first identified in termites (Luchetti, 2005; Luchetti and Mantovani, 2009). Talua reference sequences (Dataset S1) were mapped to the TE annotations using BLASTn (e-value 1e-5). If the alignment contained more than 50% of the Talua domain, the TE was considered to be a Talua containing TE. In total, we found 1575 and 4385 Talua containing TEs in the Z. nevadensis and M. natalensis genome, respectively.
Results and Discussion
Genome Architecture and Repetitive Sequences
A striking difference between ants and termites is that termite genomes are about three times larger (Table S1), which appears to be an ancestral cockroach characteristic (always several Gbs; Koshikawa et al., 2008). Termites actually have smaller genomes than cockroaches and it has been hypothesized that sociality was in fact associated with a reduction in genome size (Koshikawa et al., 2008). Yet the socially more complex M. natalensis has a genome size that is more than twice the genome size of Z. nevadensis (1.31 Gb vs. 562 Mb), which has the smallest genome known for any termite so far (Koshikawa et al., 2008). On the other hand, ant genome size appears to vary relatively little around an average of 300 Mb, with the largest ant genome published so far being 352 Mb (the red fire ant Solenopsis invicta) and smallest genome being 219 Mb (the Argentine ant Linepithema humile) (Table S2).
The two termite assemblies covered over 85% of the genomes, so any differences observed are unlikely to be related to the slightly fewer protein coding genes in Z. nevadensis (15,876 vs. 16,310 in M. natalensis). However, the M. natalensis genome contained a much higher proportion of repeat sequences (67.1 vs. 26.0% in Z. nevadensis) (Table 2). Subtracting these repeat sequences leads to comparable respective genome sizes of 367 and 365 Mb. Further genomic data will be needed to find out whether these ca. 365 Mbs represent a kind of “core genome” for termites and whether additional variation in genome size would then only be due to variation in repeat sequences. It will also be interesting to evaluate the first cockroach genomes to see whether their huge genomes (multiple Gbs) are associated with a higher number of coding or repeat sequences. In ants, genome-wide repeat content so far varies between 11.5 and 28.0% (Gadau et al., 2012) and no overall correlation with genome size appears to exist.
The M. natalensis genome had almost twice as many TEs (transposable elements) than the Z. nevadensis genome (45.9 vs. 27.8%; Table 2) and most of these were LINEs (long interspersed nuclear elements), which accounted for 20% of the M. natalensis genome (Table 2). According to the Rebase classification, most LINEs in M. natalensis resemble BovB retrotransposons, accounting for 16% of the genome, while LINEs contribute only ca. 3% in Z. nevadensis (Table 2). BovBs are relatively well known from vertebrates where they have a patchy distribution in squamates, monotremes, marsupials, ruminants, and several African mammals (Afrotheria), possibly as a consequence of horizontal gene transfer via reptile ticks (Walsh et al., 2013). In ruminants, part of one BovB LINE seems to have been recruited into a functional gene after duplication (Iwashita et al., 2006), but whether similar cooption processes may have occurred in termites remains to be explored.
The M. natalensis genome appears to have fewer SINEs (short interspersed nuclear elements) than the Z. nevadensis genome (3.6 vs. 0.2%). A new SINE retrotransposon, Talua, has recently been described for termites (Luchetti, 2005; Luchetti and Mantovani, 2009). It belongs to a new family of tRNA-derived elements that are very G+C-rich (55–60%) but makes up only a small proportion of the termite genomes (0.25 and 0.19% in Z. nevadensis and M. natalensis, respectively; Table S3). There are multi-copy TEs that are present in both termite genomes that do not resemble any known TEs. They may thus be novel termite-specific TEs, but additional termite and non-termite genomes will be needed to test this against a null hypothesis of being more general TEs that also occur in other hemimetabolous insects.
TE sequence divergence (i.e., percentage of different base pairs) relative to TE consensus sequences showed a peak at about 25% for both M. natalensis and Z. nevadensis (Figure 3), but M. natalensis had an additional divergence rate peak at ca. 7~8% (Figure 3). This might indicate that the lineage leading to M. natalensis has undergone a genome expansion that multiplied TE copies and BovB retrotransposons, which could then explain why the M. natalensis genome is so much larger than the Z. nevadensis genome.
Figure 3. The distribution of sequence divergence rates of transposable elements (TEs) as percentages of the genome size of M. natalensis (left) and Z. nevadensis (right).
Consistent with the high prevalence of repeat sequences, IPR annotation results showed a functional enrichment of DNA/RNA cutting genes in termites (Ribonuclease H domain: 22 genes, Ribonuclease H-like domain: 26 genes, endonuclease/exonuclease/phosphatase: 26 genes) compared to other insects (Table S4). Strikingly, M. natalensis had at least twice as many of such transposon-related genes than Z. nevadensis, supporting the idea that selfish replicating elements played a major role in the evolution of termite genome architecture and size (Tables S5, S6).
Cluster analyses of caste-specific transcriptomes in Z. nevadensis revealed that several of these DNA/RNA-cutting genes are overexpressed in the nymphal stages (i.e., instars with wing buds) compared to all other stages and castes (Terrapon et al., 2014). Nymphs are individuals destined to develop into winged dispersing reproductives, suggesting that TE activity might be linked to maturation processes such as gonad development. Such functions remain speculative at this point, but would be consistent with TEs having been coopted to fulfill host functions and to play fundamental roles in epigenetic regulation in organisms as different as Arabidopsis thaliana plants, Caenorhabditis elegans worms, Drosophila melanogaster flies and Mus musculus house mice (e.g., Lippman et al., 2004; Slotkin and Martienssen, 2007; Fedoroff, 2012). Silenced TEs are often activated through stressful environmental conditions (Slotkin and Martienssen, 2007; Fedoroff, 2012). In wood-dwelling termites, such conditions may arise by reduced food availability or possibly parasite pressure inducing higher rates of nymphal (sexual dispersing) development (Lenz, 1994; Korb and Schmidinger, 2004; Korb and Fuchs, 2006). Hence, it may be interesting to test whether a similar link between TE activity and stressful conditions exists during nymphal development.
Whether TEs can also be linked with epigenetic regulation of gene expression through DNA methylation (Lippman et al., 2004; Slotkin and Martienssen, 2007; Fedoroff, 2012) remains to be seen. DNA methylation has been proposed to regulate caste differentiation (Kucharski et al., 2008; Elango et al., 2009; Gadau et al., 2012; Terrapon et al., 2014) and the complete epigenetic toolbox was indeed identified in Z. nevadensis with orthologs of DNMT1 and DNMT3 (Terrapon et al., 2014). However, in M. natalensis only DNMT1 (and possibly DNMT2) could be confirmed, but not DNMT3.
Termite-specific expansions for gene families were also found among chemoperception genes that are important for communication (Table S4). Given the disparate social systems of Z. nevadensis and M. natalensis, differences in expansions of such genes may be related to divergent communication systems. Chemoperception genes mainly comprise four families: Odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), and odorant binding proteins. ORs mostly control for the specificity and sensitivity of insect olfaction. GRs are primarily involved in contact chemoperception and IRs belong to a recently discovered gene family for olfaction and gustation in Drosophila (Benton et al., 2009; Grosjean et al., 2011; Rytz et al., 2013). Odorant binding proteins primarily shuttle such compounds through the hydrophilic environment of the sensory lymph to the receptors.
The IR family is most consistently expanded in Z. nevadensis, representing the highest known value in insects (Terrapon et al., 2014). This IR number was between 4 and 10-fold higher in Z. nevadensis than in eusocial Hymenopterans, but the 80 intact GR genes remained within the overall range of 10–97 known from ants and honeybees (Zhou et al., 2012). The number of OR genes in Z. nevadensis was between one third and one half of the numbers normally found in the ants (Zhou et al., 2012), consistent with the lifestyle of wood-dwelling termites likely requiring lower levels of olfactory communication.
Overall, we found termite-specific enrichment in all four major gene families relating to olfaction (Table S4). Most IPR enrichment occurred in the ionotropic glutamate receptors that include IR genes (21). Significant enrichment was also found in ORs (7), GRs (7 TM chemoreceptor: 7), and various odorant-binding proteins (9, 7, 5). Direct comparison between Z. nevadensis and M. natalensis (Table S6) showed that Z. nevadensis had significantly more genes related to chemical communication than M. natalensis (Table S7). However, chemoperception genes are notoriously difficult to assemble and annotate (Terrapon et al., 2014), so this difference should be considered with caution, also because these genes were manually annotated in Z. nevadensis (with support from antennal RNAseq data), but automatically in M. natalensis. More work will therefore be needed before solid conclusions on the relative role of ORs in different termite species can be drawn.
Both termite species live in potentially pathogen-rich habitats. Z. nevadensis nests in decaying wood with abundant fungal growth that has probably selected for intensive allogrooming behaviors (Korb et al., 2012a). Also M. natalensis is potentially exposed to many pathogens both from its soil-nesting habitat and across its foraging range. Macrotermes species are known to protect their Termitomyces fungal symbiont from being overgrown by other fungi (Nobre et al., 2011) and termite-specific antimicrobial peptides (AMPs) have been described in another genus of fungus-growing termites (Lamberty et al., 2001).
Relative to ants and other insects, we did not find enrichments for immune defense genes in the two termite genomes and neither were there substantial differences between the two termite genomes (Tables S4, S6). All of the immune-related pathways, including pattern recognition, signaling, and gene regulation (as described for Drosophila melanogaster and other insects; Hoffmann, 2003; Hultmark, 2003; Schmid-Hempel, 2005) are present in both termite genomes (Table S8). Only two differences are noteworthy (Table 3). First, Z. nevadensis has 6 gram-negative binding proteins (GNBPs), whereas only four of these were recovered in M. natalensis. These four GNBPs are all termite-specific (Figure 4) and some of them were previously shown to be under positive selection in several Nasutitermes species, especially in species with arboreal nests (Bulmer and Crozier, 2006). The Macrotermes genome seems to lack the insect-typical GNBP duplicate and one GNBP gene that has so far only been found in Z. nevadensis (Figure 4). Second, while AMPs were not enriched in either termite genome (Table S4), their identities were completely different with Z. nevadensis having 2 AMPs and M. natalensis having 3 other AMPs (Table 3, Table S8). M. natalensis has a termite-specific defensin-like gene termicin, a category of genes that seem to have duplicated repeatedly during the radiation of Nasutitermes termites (Bulmer and Crozier, 2004). After duplication, one copy seems to often be under strong selection, while the other evolves toward neutrality (Bulmer and Crozier, 2004; Bulmer et al., 2010). Also in the soil-foraging Reticulitermes species these genes seem to be under positive selection (Bulmer et al., 2010).
Table 3. Gram-negative binding proteins (GNBPs) and anti-microbial-peptide (termicin) genes known from different termites.
Figure 4. Phylogeny of gram-negative binding proteins (GNBPs) constructed with PhylML v3.0 (LG substitution model with 100 bootstrap replicates) after alignment of the peptide sequences in ClustalW2.
In contrast to other insects where AMP production is normally induced, these genes seem to be constitutively expressed in fungus-growing termites, as has been shown for Pseudacanthotermes spiniger (Lamberty et al., 2001), which might be an adaptation to protect the symbiont against competing fungi. Termicin and other defensins (Table S8) were absent in Z. nevadensis but this species has GNBPs that are differentially expressed between castes (Terrapon et al., 2014) and may thus serve a similar function in protecting the nest from fungal infections. For the arboreal nesting termite Nasutitermes corniger it has been shown that GNBP2 has (1,3)-glucanase effector activity and functions as an antifungal agent (Bulmer et al., 2009). It is incorporated in the nest building material, where it cleaves and releases pathogenic components while priming termites for improved antimicrobial defense (Bulmer et al., 2009). Such a defensive strategy is likely to be most effective for termites with closed nests, consistent with positive selection on GNBP being most pronounced in Nasutitermes that live in arboreal nests (Bulmer and Crozier, 2006). Hence, antifungal stategies might differ in termites with different habitats; with GNBPs and termicin possibly playing complementary roles. This is supported by the fact that GNBPs in subterranean, foraging Reticulitermes species evolve neutrally while termicin was shown to have been under strong positive selection in these species (Table 3).
We can reject the possible alternative hypothesis that different defense strategies are linked to the gut symbionts that need different defense strategies to protect the symbiotic partner. As lower termites harbor protists as well as bacteria, while higher termites only have bacteria, we would then have expected higher termites having more AMPs and lower termites more GNBPs, but this is not the case because lower Reticulitermes termites have positively selected termicins. If there is an association between nesting habit and defense strategy, we expect that GNBPs are under positive selection in other wood-dwelling termites, and termicins are selected in soil-foraging termites. Additional genomic data, particularly for wood-dwelling termites, would be needed to validate this hypothesis.
Reduced numbers of immune defense genes were found in ants and the honeybee (Evans et al., 2006; Gadau et al., 2012) but also here there seems to be selection on some of the AMP genes. Similar to termicin, positive selection was detected on defensin in ants (Viljakainen and Pamilo, 2008), but this gene was not overexpressed after experimental fungal infections of leaf-cutting ant colonies, whereas two other AMPs were (Yek et al., 2013). This contrasts with dipterans (Drosophila and Anopheles) for which no evidence was found for positive selection on any AMPs (Sackton et al., 2007; Simard et al., 2007), but instead for immune recognition and signaling proteins (Schlenke and Begun, 2003; Jiggins and Kim, 2005; Sackton et al., 2007). This provides further support for the hypothesis that social insects have responded differently to selection pressure caused by microbial pathogens than solitary insects (Viljakainen and Pamilo, 2008).
Compared to M. natalensis, the Z. nevadensis genome is enriched in genes that are related to male fertility/spermatogenesis (e.g., KLHL10) (Table 4, Table S7). This suggests that the co-expansion (and co-expression) of these genes in Z. nevadensis is not typical for termite sociality but rather taxon-specific. It might be linked to the seasonal reproduction of this temperate zone species where spermatogenesis is cyclically switched on and off, which contrasts with tropical Macrotermes males that produce offspring all year round. However, some members of two spermatogenesis-related gene families, seven-in-absentia (SINA) proteins and α-tubulins, do not show Z. nevadensis-specific expansions.
Table 4. Number of genes related to spermatogenesis in Z. nevadensis and M. natalensis based on Pfam domains.
An alternative evolutionary explanation could be that males of wood-dwelling termites have low but consistent probabilities to face sperm competition when neighboring colonies merge after colony foundation. Such mergers are impossible in foraging termites where unrelated males never compete for inseminating the same queen (Boomsma, 2013). This hypothesis would predict no difference between temperate and tropical wood-dwelling termites, but a series of termite genomes will be needed to test these contentions.
The ancestral termite gut microbiota was derived from a cockroach ancestor, but major subsequent changes occurred, most notably when the higher termites evolved (Dietrich et al., 2014). The guts of the wood-dwelling termites are dominated by protists that appear to be primarily adapted to break down wood (Cleveland, 1923; Brugerolle and Radek, 2006), with complementary roles of bacteria that are often symbiotic with the gut-flagellates (Dietrich et al., 2014). The common ancestor of the evolutionarily derived Termitidae lost these flagellate symbionts so their gut microbiotas became dominated by bacteria, which may have facilitated their dietary diversification (Brune and Ohkuma, 2011; Dietrich et al., 2014). The single origin of fungiculture by the Macrotermitinae led to Termitomyces taking over primary plant decomposition and the gut microbiota shifting phylogenetically and functionally to perform complementary roles (Liu et al., 2013; Dietrich et al., 2014; Otani et al., 2014; Poulsen et al., 2014).
Changes in symbiont associations are tightly associated with termite life styles (for a recent review on termite gut symbionts, see Brune, 2014), but this may hardly induce structural genomic changes in the termite hosts, consistent with the similar gene repertoires for plant biomass decomposition found in the two termite genomes (Poulsen et al., 2014). A comparison of carbohydrate-active enzyme (CAZy) profiles of the two termite species showed a reduction in the absolute number of glycoside hydrolase enzymes (85) in M. natalensis compared to Z. nevadensis (97) (Table 5), but very similar relative abundances of specific enzyme families (Poulsen et al., 2014). Profile similarities suggest that plant-biomass decomposition genes may be ancestrally conserved across the termites, but additional termite genomes are needed to shed light on this. Such additional genomic work will need to be accompanied by enzyme function validations to test whether differences in absolute numbers reflect changes in the relative importance of termite-derived enzymes.
Table 5. Number of glycoside hydrolases of different GH families identified in Z. nevadensis and M. natalensis (from Table S28; Poulsen et al., 2014).
Despite the striking differences in social complexity between Z. nevadensis and M. natalensis we did not find major differences in gene composition. The gene families underlying chemical communication seem not to be expanded in the more complex fungus-growing termite compared to Z. nevadensis. The major differences between the two termite genomes are related to genome architecture and the presence of transposons that can explain the much larger genome size of M. natalensis. Whether these ancestrally selfish elements have been domesticated for functions related to the increased social complexity of M. natalensis needs further work. Our comparison allowed us to generate hypotheses that can be tested with functional genomic studies and with more advanced comparative analyses as more termite genomes become available.
We have highlighted the contours of further testable predictions concerning TE number and genome size, male fertility, and habitat-specific disease pressure. For any next termite genome to be sequenced (Figure 1), authors should ask questions like: (1) Is the habitat of this (e.g., drywood) termite more disease-ridden than the habitat of a comparable dampwood termite such as Z. nevadensis? (2) Would this tropical new wood-dwelling termite have similar gene family expansions for male fertility as Z. nevadensis? (3) Has this arboreal higher (e.g., Nasutitermes) termite lost specific immune defenses that match the disease pressure of its habitat and is it equally burdened by TEs as Macrotermes natalensis?
While two genomes are a major achievement in some sense, these genomes also leave us with insufficient resolution to move much beyond the crude comparisons that we offer in this paper, because Z. nevadensis and M. natalensis differ in too many evolutionary and ecological factors (Table 1). It has also become clear from comparative ant genomics that gene expression mechanisms may be more informative than structural gene differences (Simola et al., 2013). Finally, apart from obtaining more termite genomes and population genomic studies on gene expression and signatures of selection, it will also be crucially important to obtain a Cryptocercus cockroach sister lineage genome and more distant outgroup genomes for non-social hemimetabolous insects. Many surprises will likely be waiting in the wings, as both the pea aphid and the body louse genomes turned out to be unusual because of the specialized feeding habits of these insects with or without symbionts.
Conflict of Interest Statement
The Associate Editor Júrgen Rudolf Gadau declares that, despite being affiliated with the same institute and having collaborated with the author Jürgen Liebig, the review process was handled objectively. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Judith Korb was supported by research grant from the German Science Foundation (DFG; KO1895/6), Michael Poulsen by a STENO grant from The Danish Council for Independent Research Natural Sciences, Jacobus J. Boomsma by a Danish National Research Foundation grant (DNRF57), Guojie Zhang by a Marie Curie International Incoming Fellowship (300837), and Jürgen Liebig by the Agriculture and Food Research Initiative (2007-35302-18172 to Jürgen Liebig and Colin S. Brent). We thank the three referees for helpful comments and the editors Jürgen Gadau and Greg Hunt for inviting us to contribute to this special journal issue.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fgene.2015.00009/abstract
Adams, M. D., Celniker, S. E., Holt, R. A., Evans, C. A., Gocayne, J. D., Amanatides, P. G., et al. (2000). The genome sequence of Drosophila melanogaster. Science 287, 2185–2195. doi: 10.1126/science.287.5461.2185
Benton, R., Vannice, K. S., Gomez-Diaz, C., and Vosshall, L. B. (2009). Variant ionotropic glutamate receptors as chemosensory receptors in Drosophila. Cell 136, 149–162. doi: 10.1016/j.cell.2008.12.001
Bonasio, R., Zhang, G., Ye, C., Mutti, N. S., Fang, X., Qin, N., et al. (2010). Genomic comparison of the ants Camponotus floridanus and Harpegnathos saltator. Science 329, 1068–1071. doi: 10.1126/science.1192428
Boomsma, J. J., Schmid-Hempel, P., and Hughes, W. O. H. (2005). “Life histories and parasite pressure across the major groups of social insects,” in Insect Evolutionary Ecology, eds M. D. E. Fellowes, G. J. Holloway, and J. Rolff (Wallingford, CT: CABI), 139–175.
Brugerolle, G., and Radek, R. (2006). “Symbiotic protozoa of termites,” in Intestinal Microorganisms of Termites and Other Invertebrates, eds H. König and A. Varma (Berlin, Heidelberg: Springer), 243–269. doi: 10.1007/3-540-28185-1_10
Brune, A., and Ohkuma, M. (2011). “Role of the termite gut microbiota in symbiotic digestion,” in Biology of Termites: A Modern Synthesis, eds D. E. Bignell, Y. Roisin, and N. Lo (Dordrecht; Heidelberg; London; New York: Springer), 439–476.
Bulmer, M. S., Bachelet, I., Raman, R., Rosengaus, R. B., and Sasisekharan, R. (2009). Targeting an antimicrobial effector function in insect immunity as a pest control strategy. Proc. Natl. Acad. Sci. U.S.A. 106, 12652–12657. doi: 10.1073/pnas.0904063106
Dietrich, C., Köhler, T., and Brune, A. (2014). The cockroach origin of the termite gut microbiota: patterns in bacterial community structure reflect major evolutionary events. Appl. Environ. Microb. 80, 2261–2269. doi: 10.1128/AEM.04206-13
Elango, N., Hunt, B. G., Goodisman, M. A. D., and Yi, S. V. (2009). DNA methylation is widespread and associated with differential gene expression in castes of the honeybee, Apis mellifera. Proc. Natl. Acad. Sci. U.S.A. 106, 11206–11211. doi: 10.1073/pnas.0900301106
Evans, J. D., Aronstein, K., Chen, Y. P., Hetru, C., Imler, J. L., Jiang, H., et al. (2006). Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Mol. Biol. 15, 645–656. doi: 10.1111/j.1365-2583.2006.00682.x
Gadau, J., Helmkampf, M., Nygaard, S., Roux, J., Simola, D. F., Smith, C. R., et al. (2012). The genomic impact of 100 million years of social evolution in seven ant species. Trends Genet. 28, 14–21. doi: 10.1016/j.tig.2011.08.005
Grosjean, Y., Rytz, R., Farine, J. P., Abuin, L., Cortot, J., Jefferis, G., et al. (2011). An olfactory receptor for food-derived odours promotes male courtship in Drosophila. Nature 478, 236–240. doi: 10.1038/nature10428
Howard, K. J., and Thorne, B. L. (2011). “Eusocial evolution in termites and hymenoptera,” in Biology of Termites: A Modern Synthesis, eds D. E. Bignell, Y. Roisin, and N. Lo (Dordrecht; Heidelberg; London; New York: Springer), 97–132.
Hughes, W. O. H., Oldroyd, B. P., Beekman, M., and Ratnieks, F. L. W. (2008). Ancestral monogamy shows kin selection is key to the evolution of eusociality. Science 320, 1212–1216. doi: 10.1126/science.1156108
Inward, D., Beccaloni, G., and Eggleton, P. (2007a). Death of an order: a comprehensive molecular phylogenetic study confirms that termites are eusocial cockroaches. Biol. Lett. 3, 331–335. doi: 10.1098/rsbl.2007.0102
Inward, D. J. G., Vogler, A. P., and Eggleton, P. (2007b). A comprehensive phylogenetic analysis of termites (Isoptera) illuminates key aspects of their evolutionary biology. Mol. Phylogenet. Evol. 44, 953–967. doi: 10.1016/j.ympev.2007.05.014
Iwashita, S., Ueno, S., Nakashima, K., Song, S. Y., Ohshima, K., Tanaka, K., et al. (2006). A tandem gene duplictaion followed by recrutitment of a retrotransposon created the paralogous bucentaur gene (bcnt) in the ancestral ruminant Mol. Biol. Evol. 23, 798–806. doi: 10.1093/molbev/msj088
Kambhampati, S., and Eggleton, P. (2000). “Taxonomy and phylogeny of termites,” in Termites: Evolution, Sociality, Symbiosis and Ecology, eds T. Abe, D. E. Bignell, and M. Higashi (Dordrech: Kluwer Academic Publishers), 1–23. doi: 10.1007/978-94-017-3223-9_1
Kirkness, E. F., Haas, B. J., Sun, W. L., Braig, H. R., Perotti, M. A., Clark, J. M., et al. (2010). Genome sequences of the human body louse and its primary endosymbiont provide insights into the permanent parasitic lifestyle. Proc. Natl. Acad. Sci. U.S.A. 107, 12168–12173. doi: 10.1073/pnas.1003379107
Korb, J. (2011). “Termite mound architecture, from function to construction,” in Biology of Termites: A Modern Synthesis, eds D. E. Bignell, Y. Roisin, and N. Lo (Dordrecht; Heidelberg; London; New York: Springer), 349–376.
Korb, J., and Hartfelder, K. (2008). Life history and development - a framework for understanding the ample developmental plasticity in lower termites. Biol. Rev. 83, 295–313. doi: 10.1111/j.1469-185X.2008.00044.x
Korb, J., Hoffmann, K., and Hartfelder, K. (2012b). Molting dynamics and juvenile hormone titer profiles in the nymphal stages of a lower termite, Cryptotermes secundus (Kalotermitidae) - signatures of developmental plasticity. J. Insect Physiol. 58, 376–383. doi: 10.1016/j.jinsphys.2011.12.016
Koshikawa, S., Miyazaki, S., Cornette, R., Matsumoto, T., and Miura, T. (2008). Genome size of termites (Insecta, Dictyoptera, Isoptera) and wood roaches (Insecta, Dictyoptera, Cryptocercidae). Naturwissenschaften 95, 859–867. doi: 10.1007/s00114-008-0395-7
Lamberty, M., Zachary, D., Lanot, R., Bordereau, C., Robert, A., Hoffmann, J. A., et al. (2001). Insect immunity - Constitutive expression of a cysteine-rich antifungal and a linear antibacterial peptide in a termite insect. J. Biol. Chem. 276, 4085–4092. doi: 10.1074/jbc.M002998200
Lenz, M. (1994). “Food resources, colony growth and caste development in wood-feeding termites,” in Nourishment and Evolution in Insect Societies, eds J. Hunt and C. A. Nalepa (New Delhi: Oxford and IBH Publishing Co. Prt. Ltd), 159–209.
Li, H., Coghlan, A., Ruan, J., Coin, L. J., Heriche, J. K., Osmotherly, L., et al. (2006). TreeFam: a curated database of phylogenetic trees of animal gene families. Nucleic Acids Res. 34, 572–580. doi: 10.1093/nar/gkj118
Liebig, J. (2010). “Hydrocarbon profiles indicate fertility and dominance status in ant, bee, and wasp colonies,” in Insect Hydrocarbons. Biology, Biochemistry, and Chemical Ecology, eds G. J. Blomquist and A. G. Bagneres (Cambridge: Cambridge UP), 254–281. doi: 10.1017/CBO9780511711909.014
Liebig, J., Eliyahu, D., and Brent, C. S. (2009). Cuticular hydrocarbon profiles indicate reproductive status in the termite Zootermopsis nevadensis. Behav. Ecol. Sociobiol. 63, 1799–1807. doi: 10.1007/s00265-009-0807-5
Lippman, Z., Gendrel, A.-V., Black, M., Vaughn, M. W., Dedhia, N., McCombie, W. R., et al. (2004). Role of transposable elements in heterochromatin and epigenetic control. Nature 430, 471–473. doi: 10.1038/nature02651
Liu, N., Zhang, L., Zhou, H., Zhang, M., Yan, X., Wang, Q., et al. (2013). Metagenomic insights into metabolic capacities of the gut microbiota in a fungus-cultivating termite (Odontotermes yunnanensis). PLoS ONE 8:e69184. doi: 10.1371/journal.pone.0069184
Luchetti, A., and Mantovani, B. (2009). Talua SINE biology in the genome of the Reticulitermes subterranean termites (Isoptera, Rhinotermitidae). J. Mol. Evol. 69, 589–600. doi: 10.1007/s00239-009-9285-7
Miura, T., and Scharf, M. E. (2011). “Molecular basis underlying caste differentiation in termites,” in Biology of Termites: A Modern Synthesis, eds D. E. Bignell, Y. Roisin, and N. Lo (Dordrecht; Heidelberg; London; New York: Springer), 211–254.
Moriya, Y., Itoh, M., Okuda, S., Yoshizawa, A. C., and Kanehisa, M. (2007). KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 35, W182–W185. doi: 10.1093/nar/gkm321
Nobre, T., Rouland-Lefevre, C., and Aanen, D. K. (2011). “Comparative biology of fungus cultivation in termites and ants,” in Biology of Termites: A Modern Synthesis, eds D. E. Bignell, Y. Roisin, and N. Lo (Dordrecht; Heidelberg; London; New York: Springer), 193–210.
Noirot, C., and Darlington, J. P. E. C. (2000). “Termite nests: architecture, regulation and defence,” in Termites: Evolution, Sociality, Symbiosis and Ecology, eds T. Abe, D. E. Bignell, and M. Higashi (Dordrecht: Kluwer Academic Publishers), 121–139. doi: 10.1007/978-94-017-3223-9_6
Nygaard, S., Zhang, G., Schiott, M., Li, C., Wurm, Y., Hu, H., et al. (2011). The genome of the leaf-cutting ant Acromyrmex echinatior suggests key adaptations to advanced social life and fungus farming. Genome Res. 21, 1339–1348. doi: 10.1101/gr.121392.111
Okot-Kotber, B. M. (1985). “Mechanisms of caste determination in a higher termite, Macrotermes michaelseni (Isoptera, Macrotermitidae),” in Caste Differentiation in Social Insects, eds J. A. L. Watson, B. M. Okot-Kotber, and C. Noirot (Oxford: Pergamon Press), 267–306. doi: 10.1016/B978-0-08-030783-1.50026-4
Otani, S., Mikaelyan, A., Nobre, T., Hansen, L. H., Koné, N. A., Sørensen, S. J., et al. (2014). Identifying the core microbial community in the gut of fungus-growing termites. Mol. Ecol. 23, 4631–4644. doi: 10.1111/mec.12874
Poulsen, M., Hu, H., Li, C., Chen, Z., Xu, L., Otani, S., et al. (2014). Complementary symbiont contributions to plant decomposition in a fungus-farming termite. Proc. Natl. Acad. Sci. U.S.A. 111, 14500–14505. doi: 10.1073/pnas.1319718111
Richards, S., Gibbs, R. A., Weinstock, G. M., Brown, S. J., Denell, R., Beeman, R. W., et al. (2008). The genome of the model beetle and pest Tribolium castaneum. Nature 452, 949–955. doi: 10.1038/nature06784
Roisin, Y. (2000). “Diversity and evolution of caste patterns,” in Diversity and Evolution of Caste Patterns, eds T. Abe, D. E. Bignell, and M. Higashi (Dordrecht: Kluwer Academic Publishers), 95–119.
Rosengaus, R. B., Traniello, J. F. A., and Bulmer, M. S. (2011). “Ecology, behavior and evolution of disease resistance in termites,” in Biology of Termites: A Modern Synthesis, eds D. E. Bignell, Y. Roisin, and L. Lo (Dordrecht; Heidelberg; London; New York: Springer), 165–192.
Rytz, R., Croset, V., and Benton, R. (2013). Ionotropic Receptors (IRs): chemosensory ionotropic glutamate receptors in Drosophila and beyond. Insect Biochem. Mol. Biol. 43, 888–897. doi: 10.1016/j.ibmb.2013.02.007
Sackton, T. B., Lazzaro, B. P., Schlenke, T. A., Evans, J. D., Hultmark, D., and Clark, A. G. (2007). Dynamic evolution in the innate immunesystem in Drosophila. Nat. Genet. 39, 1461–1468. doi: 10.1038/ng.2007.60
Shellman-Reeve, J. S. (1997). “The spectrum of eusociality in termites,” in The Evolution of Social Behavior in Insects and Arachnids, eds J. C. Choe and B. J. Crespi (Cambridge: Cambridge University Press), 52–93.
Simard, F., Licht, M., Besansky, N. J., and Lehmann, T. (2007). Polymorphism at the defensin gene in the Anopheles gambiae complex: testing different selection hypotheses. Infect. Genet. Evol. 7, 285–292. doi: 10.1016/j.meegid.2006.11.004
Simola, D. F., Wissler, L., Donahue, G., Waterhouse, R. M., Helmkampf, M., Roux, J., et al. (2013). Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality. Genome Res. 23, 1235. doi: 10.1101/gr.155408.113
Smit, A. F. A., Hubley, R., and Green, P. (1996-2010). RepeatMasker Open-3.0. Available online at: http://www.repeatmasker.org
Suen, G., Teiling, C., Li, L., Holt, C., Abouheif, E., Bornberg-Bauer, E., et al. (2011). The genome sequence of the leaf-cutter ant Atta cephalotes reveals insights into its obligate symbiotic lifestyle. PLoS Genet. 7:e1002007. doi: 10.1371/journal.pgen.1002007
Terrapon, N., Li, C., Robertson, H. M., Ji, L., Meng, X., Booth, W., et al. (2014). Molecular traces of alternative social organization in a termite genome. Nat. Commun. 5, e3636. doi: 10.1038/ncomms4636
van Oystaeyen, A., Oliveira, R. C., Holman, L., van Zweden, J. S., Romero, C., Oi, C. A., et al. (2014). Conserved class of queen pheromones stops social insect workers from reproducing. Science 343, 287–290. doi: 10.1126/science.1244899
Van Zweden, J. S., and D'Ettorre, P. (2010). “Nestmate recognition in social insects and the role of hydrocarbons,” in Insect Hydrocarbons: Biology, Biochemistry, and Chemical Ecology, eds G. J. Blomquist and A.G. Bagneres (Cambrdige UK, Cambridge University Press), 222–243.
Walsh, A. M., Kortschak, R. D., Gardner, M. G., Bertozzi, T., and Adelson, D. L. (2013). Widespread horizontal transfer of retrotransposons. Proc. Natl. Acad. Sci. U.S.A. 110, 1012–1016. doi: 10.1073/pnas.1205856110
Weil, T., Hoffmann, K., Kroiss, J., Strohm, E., and Korb, J. (2009). Scent of a queen-cuticular hydrocarbons specific for female reproductives in lower termites. Naturwissenschaften 96, 315–319. doi: 10.1007/s00114-008-0475-8
Werren, J. H., Richards, S., Desjardins, C. A., Niehuis, O., Gadau, J., Colbourne, J. K., et al. (2010). Functional and evolutionary insights from the genomes of three parasitoid Nasonia species. Science 327, 343–348. doi: 10.1126/science.1178028
Wood, T. G., and Thomas, R. J. (1989). “The mutualistic association between Macrotermitinae and Termitomyces,” in Insect-Fungus Interactions, eds N. Wilding, N. M. Collins, P. M. Hammond, and J. F. Webber (New York, NY: Academic Press), 69–92. doi: 10.1016/B978-0-12-751800-8.50009-4
Yek, S. H., Boomsma, J. J., and Schiøtt, M. (2013). Differential gene expression in Acromyrmex leaf-cutting ants after challenges with two fungal pathogens. Mol. Ecol. 22, 2173–2187. doi: 10.1111/mec.12255
Zhou, X., Slone, J. D., Rokas, A., Berger, S. L., Liebig, J., Ray, A., et al. (2012). Phylogenetic and transcriptomic analysis of chemosensory receptors in a pair of divergent ant species reveals caste-specific signatures of odor coding. PLoS Genet. 8:e1002930. doi: 10.1371/journal.pgen.1002930
Keywords: chemical communication, genomes, immunity, social organization, social insects, symbiosis, termites, transposable elements
Citation: Korb J, Poulsen M, Hu H, Li C, Boomsma JJ, Zhang G and Liebig J (2015) A genomic comparison of two termites with different social complexity. Front. Genet. 6:9. doi: 10.3389/fgene.2015.00009
Received: 29 September 2014; Paper pending published: 25 November 2014;
Accepted: 09 January 2015; Published online: 04 March 2015.
Edited by:Juergen Rudolf Gadau, Arizona State University, USA
Reviewed by:Seirian Sumner, University of Bristol, UK
Bart Pannebakker, Wageningen University, Netherlands
Michael E. Scharf, Purdue University, USA
Copyright © 2015 Korb, Poulsen, Hu, Li, Boomsma, Zhang and Liebig. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Judith Korb, Department of Evolutionary Biology and Ecology, Institute of Biology I, University of Freiburg, Hauptstrasse 1, D-79104 Freiburg, Germany e-mail: email@example.com