Search

Confidence measures of branch reliability play an important role in phylogenetics as these measures allow to identify trees or parts of a tree that are well supported by the data and thus adequate to serve as basis for evolutionary inference of biological systems. Unreliable branch relationships in phylogenetic analyses are of concern because of their potential to represent incorrect relationships of interest among more reliable branch relationships. The site-concordance factor implemented in the IQ-TREE package is a recently introduced heuristic solution to the problem of identifying unreliable branch relationships on the basis of quartets. We test the performance of the site-concordance measure with simple examples based on simulated data and designed to study its behaviour in branch support estimates related to different degrees of branch length heterogeneities among a ten sequence tree. Our results show that in particular in cases of relationships with heterogeneous branch lengths site-concordance measures may be misleading. We therefore argue that the maximum parsimony optimality criterion currently used by the site-concordance measure may sometimes be poorly suited to evaluate branch support and that the scores reported by the site-concordance factor should not be considered as reliable.

Background Phylogenetic relationships among the myriapod subgroups Chilopoda, Diplopoda, Symphyla and Pauropoda are still not robustly resolved. The first phylogenomic study covering all subgroups resolved phylogenetic relationships congruently to morphological evidence but is in conflict with most previously published phylogenetic trees based on diverse molecular data. Outgroup choice and long-branch attraction effects were stated as possible explanations for these incongruencies. In this study, we addressed these issues by extending the myriapod and outgroup taxon sampling using transcriptome data. Results We generated new transcriptome data of 42 panarthropod species, including all four myriapod subgroups and additional outgroup taxa. Our taxon sampling was complemented by published transcriptome and genome data resulting in a supermatrix covering 59 species. We compiled two data sets, the first with a full coverage of genes per species (292 single-copy protein-coding genes), the second with a less stringent coverage (988 genes). We inferred phylogenetic relationships among myriapods using different data types, tree inference, and quartet computation approaches. Our results unambiguously support monophyletic Mandibulata and Myriapoda. Our analyses clearly showed that there is strong signal for a single unrooted topology, but a sensitivity of the position of the internal root on the choice of outgroups. However, we observe strong evidence for a clade Pauropoda+Symphyla, as well as for a clade Chilopoda+Diplopoda. Conclusions Our best quartet topology is incongruent with current morphological phylogenies which were supported in another phylogenomic study. AU tests and quartet mapping reject the quartet topology congruent to trees inferred with morphological characters. Moreover, quartet mapping shows that confounding signal present in the data set is sufficient to explain the weak signal for the quartet topology derived from morphological characters. Although outgroup choice affects results, our study could narrow possible trees to derivatives of a single quartet topology. For highly disputed relationships, we propose to apply a series of tests (AU and quartet mapping), since results of such tests allow to narrow down possible relationships and to rule out confounding signal.

The life-prolonging effects of antioxidants have long entered popular culture, but the scientific community still debates whether free radicals and the resulting oxidative stress negatively affect longevity. Social insects are intriguing models for analysing the relationship between oxidative stress and senescence because life histories differ vastly between long-lived reproductives and the genetically similar but short-lived workers. Here, we present the results of an experiment on the accumulation of oxidative damage to proteins, and a comparative analysis of the expression of 20 selected genes commonly involved in managing oxidative damage, across four species of social insects: a termite, two bees and an ant. Although the source of analysed tissue varied across the four species, our results suggest that oxidative stress is a significant factor in senescence and that its manifestation and antioxidant defenses differ among species, making it difficult to find general patterns. More detailed and controlled investigations on why responses to oxidative stress may differ across social species may lead to a better understanding of the relations between oxidative stress, antioxidants, social life history and senescence.This article is part of the theme issue ‘Ageing and sociality: why, when and how does sociality change ageing patterns?'

Eusocial insect queens are remarkable in their ability to maximise both fecundity and longevity, thus escaping the typical trade-off between these two traits. In species exhibiting complex eusocial behaviour, several mechanisms have been proposed to underlie the remoulding of the trade-off, such as reshaping of the juvenile hormone pathway, or caste-specific susceptibility to oxidative stress. However, it remains a challenge to disentangle the molecular mechanisms underlying the remoulding of the trade-off in eusocial insects from caste-specific physiological attributes that have subsequently arisen due to their different life histories. Socially plastic species such as the orchid beeEuglossa viridissimarepresent excellent models to address the role of socialityper sein longevity as they allow direct comparisons of solitary and social individuals within a common genetic background. We present data on gene expression and juvenile hormone levels from young and old bees, from both solitary and social nests. We found 940 genes to be differentially expressed with age in solitary females,versusonly 14 genes in social dominant females, and seven genes in subordinate females. We performed a weighted gene co-expression network analysis to further highlight candidate genes related to ageing in this species. Primary “ageing gene” candidates were related to protein synthesis, gene expression, immunity and venom production. Remarkably, juvenile hormone titres did not vary with age or social status. These results represent an important step in understanding the proximate mechanisms underlying the remodeling of the fecundity/longevity trade-off that accompanies the evolutionary transition from solitary life to eusociality.Significance statementThe remarkably long lifespan of the queens of eusocial insects despite their high reproductive output suggests that they are not subject to the widespread trade-off between fecundity and longevity that governs solitary animal life histories, yet surprisingly little is known of the molecular mechanisms underpinning their longevity. Using a socially plastic bee in which some individuals of a population are social whilst others are solitary, we identified hundreds of candidate genes and related gene networks that are involved in the remoulding of the fecundity/longevity tradeoff. As well as identifying candidate ageing genes, our data suggest that even in incipient stages of sociality there is a marked reprogramming of ageing; long live the queen.

The exceptional longevity of social insect queens despite their lifelong high fecundity remains poorly understood in ageing biology. To gain insights into the mechanisms that might underlie ageing in social insects, we compared gene expression patterns between young and old castes (both queens and workers) across different lineages of social insects (two termite, two bee and two ant species). After global analyses, we paid particular attention to genes of the insulin/insulin-like growth factor 1 signalling (IIS)/target of rapamycin (TOR)/juvenile hormone (JH) network, which is well known to regulate lifespan and the trade-off between reproduction and somatic maintenance in solitary insects. Our results reveal a major role of the downstream components and target genes of this network (e.g. JH signalling, vitellogenins, major royal jelly proteins and immune genes) in affecting ageing and the caste-specific physiology of social insects, but an apparently lesser role of the upstream IIS/TOR signalling components. Together with a growing appreciation of the importance of such downstream targets, this leads us to propose the TI–J–LiFe ( T OR/ I IS– J H– Li fespan and Fe cundity) network as a conceptual framework for understanding the mechanisms of ageing and fecundity in social insects and beyond. This article is part of the theme issue ‘Ageing and sociality: why, when and how does sociality change ageing patterns?’

Background Blowflies are ubiquitous insects, often shiny and metallic, and the larvae of many species provide important ecosystem services (e.g., recycling carrion) and are used in forensics and debridement therapy. Yet, the taxon has repeatedly been recovered to be para- or polyphyletic, and the lack of a well-corroborated phylogeny has prevented a robust classification. Results We here resolve the relationships between the different blowfly subclades by including all recognized subfamilies in a phylogenomic analysis using 2221 single-copy nuclear protein-coding genes of Diptera. Maximum likelihood (ML), maximum parsimony (MP), and coalescent-based phylogeny reconstructions all support the same relationships for the full data set. Based on this backbone phylogeny, blowflies are redefined as the most inclusive monophylum within the superfamily Oestroidea not containing Mesembrinellidae, Mystacinobiidae, Oestridae, Polleniidae, Sarcophagidae, Tachinidae, and Ulurumyiidae. The constituent subfamilies are re-classified as Ameniinae (including the Helicoboscinae, syn. nov.), Bengaliinae, Calliphorinae (including Aphyssurinae, syn. nov., Melanomyinae, syn. nov., and Toxotarsinae, syn. nov.), Chrysomyinae, Luciliinae, Phumosiinae, Rhiniinae stat. rev., and Rhinophorinae stat. rev. Metallic coloration in the adult is shown to be widespread but does not emerge as the most likely ground plan feature. Conclusions Our study provides the first phylogeny of oestroid calyptrates including all blowfly subfamilies. This allows settling a long-lasting controversy in Diptera by redefining blowflies as a well-supported monophylum, and blowfly classification is adjusted accordingly. The archetypical blowfly trait of carrion-feeding maggots most likely evolved twice, and the metallic color may not belong to the blowfly ground plan.

The wasp family Chrysididae (cuckoo wasps, gold wasps) comprises exclusively parasitoid and kleptoparasitic species, many of which feature a stunning iridescent coloration and phenotypic adaptations to their parasitic life style. Previous attempts to infer phylogenetic relationships among the family’s major lineages (subfamilies, tribes, genera) based on Sanger sequence data were insufficient to statistically resolve the monophyly and the phylogenetic position of the subfamily Amiseginae and the phylogenetic relationships among the tribes Allocoeliini, Chrysidini, Elampini, and Parnopini (Chrysidinae). Here, we present a phylogeny inferred from nucleotide sequence data of 492 nuclear single-copy genes (230,915 aligned amino acid sites) from 94 species of Chrysidoidea (representing Bethylidae, Chrysididae, Dryinidae, Plumariidae) and 45 outgroup species by combining RNA-seq and DNA target enrichment data. We find support for Amiseginae being more closely related to Cleptinae than to Chrysidinae. Furthermore, we find strong support for Allocoeliini being the sister lineage of all remaining Chrysidinae, whereas Elampini represent the sister lineage of Chrysidini and Parnopini. Our study corroborates results from a recent phylogenomic investigation, which revealed Chrysidoidea as likely paraphyletic.

Background: Most insects are relatively short-lived, with a maximum lifespan of a few weeks, like the aging model organism, the fruit-fly Drosophila melanogaster. By contrast, the queens of social insects (termites, ants, some bees and wasps) can live for more than a decade. This makes social insects promising new models in aging research providing insights into how a long reproductive life can be achieved. Yet, aging studies on social insect reproductives are hampered by a lack of quantitative data on age-dependent survival and time series analyses that cover the whole lifespan of such long-lived individuals. We studied aging in queens of the drywood termite Cryptotermes secundus by determining survival probabilities over more than 15 years and performed transcriptome analyses for queens of known age that covered their whole lifespan. Results: The maximum lifespan of C. secundus queens was 13 years with a median maximum longevity of 11.0 years. Time course and co-expression network analyses of gene expression patterns over time indicated a non-gradual aging pattern. It was characterized by networks of genes that became differentially expressed only late in life, namely after an age of 10 years, which associates well with the median maximum lifespan for queens. These old-age gene networks reflect processes of physiological upheaval. We detected strong signs of stress, decline, defence and repair at the transcriptional level of epigenetic control as well as at the post-transcriptional level with changes in transposable element activity and the proteostasis network. The latter depicts an upregulation of protein degradation, together with protein synthesis and protein folding, processes which are often down-regulated in old animals. The simultaneous upregulation of protein synthesis and autophagy is indicative for a stress-response mediated by the transcription factor cnc, a homolog of human nrf genes. Conclusion: Our results show non-linear senescence with a rather sudden physiological upheaval at old-age. Most importantly, they point to a re-wiring in the proteostatis network as central for explaining the long life of social insect queens.

Insect orders have been defined and stable for decades, with few notable exceptions (e.g., Blattodea and Psocoptera). One of the few remaining questions of order-level monophyly is that of Mecoptera in respect to the phylogenetic placement of Siphonaptera (fleas). We used a large set of transcriptomic nucleotide sequence data representing 56 species and more than 3,000 single-copy genes to resolve the evolutionary history of Antliophora, including fleas (Siphonaptera), scorpionflies and relatives (Mecoptera), and true flies (Diptera). We find that fleas and mecopterans together are the sister group of flies. However, our data and/or analyses are unable to distinguish whether fleas are sister to a monophyletic Mecoptera, or whether they arose from within extant mecopteran families, rendering Mecoptera paraphyletic. We did not detect parameter bias in our dataset after applying a broad range of detection methods. Counter to a previous hypothesis that placed fleas within Mecoptera as the sister group to wingless boreids (snow fleas), we found a potential sister group relationship between fleas and the enigmatic family Nannochoristidae. Although we lack conclusive evidence, it seems possible that fleas represent the most-species rich group of modern mecopterans and that their parasitic lifestyle and morphological adaptations have simply made them unrecognizable in respect to their order-level classification.

BackgroundCladobranch sea slugs represent roughly half of the biodiversity of soft-bodied, marine gastropod molluscs (Nudibranchia) on the planet. Despite their global distribution from shallow waters to the deep sea, from tropical into polar seas, and their important role in marine ecosystems and for humans (as bioindicators and providers of medical drug leads), the evolutionary history of cladobranch sea slugs is not yet fully understood. Here, we amplify the current knowledge on the phylogenetic relationships by extending the cladobranch and outgroup taxon sampling using transcriptome data.ResultsWe generated new transcriptome data for 19 species of cladobranch sea slugs and two additional outgroup taxa. We complemented our taxon sampling with previously published transcriptome data, resulting in a final supermatrix covering 56 species from all but one accepted cladobranch superfamilies. Transcriptome assembly using six different assemblers, selection of those assemblies providing the largest amount of potentially phylogenetically informative sites, and quality-driven compilation of data sets resulted in three different supermatrices: one with a full coverage of genes per species (446 single-copy protein-coding genes) and two with a less stringent coverage (667 genes with 98.9% partition coverage and 1,767 genes with 86% partition coverage, respectively). We used these supermatrices to infer statistically robust maximum-likelihood trees. All analyses, irrespective of the data set, indicate maximum statistical support for all major splits and phylogenetic relationships on family level. The only discordance between the inferred trees is the position ofEmbletonia pulchra. Extensive testing using Four-cluster Likelihood Mapping, Approximately Unbiased tests, and Quartet Scores revealed that its position is not due to any informative phylogenetic signal, but caused by confounding signal.ConclusionsOur data matrices and the inferred trees inferred can serve as a solid foundation for future work on the taxonomy and evolutionary history of Cladobranchia. The correct placement ofE. pulchra, however, proves challenging, even with large data sets. Moreover, quartet mapping shows that confounding signal present in the data is sufficient to explain the inferred position ofE. pulchra, again leaving its phylogenetic position as an enigma.

Multiple sequence alignments (MSAs) play a pivotal role in studies of molecular sequence data, but nobody has developed a minimum reporting standard (MRS) to quantify the completeness of MSAs in terms of completely specified nucleotides or amino acids. We present an MRS that relies on four simple completeness metrics. The metrics are implemented in AliStat, a program developed to support the MRS. A survey of published MSAs illustrates the benefits and unprecedented transparency offered by the MRS.

Social insects are promising new models in aging research. Within single colonies, longevity differences of several magnitudes exist that can be found elsewhere only between different species. Reproducing queens (and, in termites, also kings) can live for several decades, whereas sterile workers often have a lifespan of a few weeks only. We studied aging in the wild in a highly social insect, the termite Macrotermes bellicosus, which has one of the most pronounced longevity differences between reproductives and workers. We show that gene-expression patterns differed little between young and old reproductives, implying negligible aging. By contrast, old major workers had many genes up-regulated that are related to transposable elements (TEs), which can cause aging. Strikingly, genes from the PIWI-interacting RNA (piRNA) pathway, which are generally known to silence TEs in the germline of multicellular animals, were down-regulated only in old major workers but not in reproductives. Continued up-regulation of the piRNA defense commonly found in the germline of animals can explain the long life of termite reproductives, implying somatic cooption of germline defense during social evolution. This presents a striking germline/soma analogy as envisioned by the superorganism concept: the reproductives and workers of a colony reflect the germline and soma of multicellular animals, respectively. Our results provide support for the disposable soma theory of aging.

The formation and spread of the Australian arid zone during the Neogene was a profoundly transformative event in the biogeographic history of Australia, resulting in extinction or range contraction in lineages adapted to mesic habitats, as well as diversification and range expansion in arid-adapted taxa (most of which evolved from mesic ancestors). However, the geographic origins of the arid zone biota are still relatively poorly understood, especially among highly diverse invertebrate lineages, many of which are themselves poorly documented at the species level. Spiny trapdoor spiders (Idiopidae: Arbanitinae) are one such lineage, having mesic ‘on-the-continent’ Gondwanan origins, while also having experienced major arid zone radiations in select clades. In this study, we present new orthologous nuclear markers for the phylogenetic inference of mygalomorph spiders, and use them to infer the phylogeny of Australasian Idiopidae with a 12-gene parallel tagged amplicon next-generation sequencing approach. We use these data to test the mode and timing of diversification of arid-adapted idiopid lineages across mainland Australia, and employ a continent-wide sampling of the fauna’s phylogenetic and geographic diversity to facilitate ancestral area inference. We further explore the evolution of phenotypic and behavioural characters associated with both arid and mesic environments, and test an ‘out of south-western Australia’ hypothesis for the origin of arid zone clades. Three lineages of Idiopidae are shown to have diversified in the arid zone during the Miocene, one (genus Euoplos) exclusively in Western Australia. Arid zone Blakistonia likely had their origins in South Australia, whereas in the most widespread genus Aganippe, a more complex scenario is evident, with likely range expansion from southern Western Australia to southern South Australia, from where the bulk of the arid zone fauna then originated. In Aganippe, remarkable adaptations to phragmotic burrow-plugging in transitional arid zone taxa have evolved twice independently in Western Australia, while in Misgolas and Cataxia, burrow door-building behaviours have likely been independently lost at least three times in the eastern Australian mesic zone. We also show that the presence of idiopids in New Zealand (Cantuaria) is likely to be the result of recent dispersal from Australia, rather than ancient continental vicariance. By providing the first comprehensive, continental synopsis of arid zone biogeography in an Australian arachnid lineage, we show that the diversification of arbanitine Idiopidae was intimately associated with climate shifts during the Neogene, resulting in multiple Mio-Pliocene radiations.

RNA interference (RNAi) refers to the set of molecular processes found in eukaryotic organisms in which small RNA molecules mediate the silencing or down-regulation of target genes. In insects, RNAi serves a number of functions, including regulation of endogenous genes, anti-viral defense, and defense against transposable elements. Despite being well studied in model organisms, such as Drosophila, the distribution of core RNAi pathway genes and their evolution in insects is not well understood. Here we present the most comprehensive overview of the distribution and diversity of core RNAi pathway genes across 100 insect species, encompassing all currently recognized insect orders. We inferred the phylogenetic origin of insect-specific RNAi pathway genes and also identified several hitherto unrecorded gene expansions using whole-body transcriptome data from the international 1KITE (1000 Insect Transcriptome Evolution) project as well as other resources such as i5K (5000 Insect Genome Project). Specifically, we traced the origin of the double stranded RNA binding protein R2D2 to the last common ancestor of winged insects (Pterygota), the loss of Sid-1/Tag-130 orthologs in Antliophora (fleas, flies and relatives, and scorpionflies in a broad sense), and confirm previous evidence for the splitting of the Argonaute proteins Aubergine and Piwi in Brachyceran flies (Diptera, Brachycera). Our study offers new reference points for future experimental research on RNAi-related pathway genes in insects.

Significance Lepidoptera play key roles in many biological systems. Butterflies are hypothesized to have evolved contemporaneously with flowering plants, and moths are thought to have gained anti-bat defenses in response to echolocating predatory bats, but these hypotheses have largely gone untested. Using a transcriptomic, dated evolutionary tree of Lepidoptera, we demonstrate that the most recent common ancestor of Lepidoptera is considerably older than previously hypothesized. The oldest moths in crown Lepidoptera were present in the Carboniferous, some 300 million years ago, and began to diversify largely in synchrony with angiosperms. We show that multiple lineages of moths independently evolved hearing organs well before the origin of bats, rejecting the hypothesis that lepidopteran hearing organs arose in response to these predators., Butterflies and moths (Lepidoptera) are one of the major superradiations of insects, comprising nearly 160,000 described extant species. As herbivores, pollinators, and prey, Lepidoptera play a fundamental role in almost every terrestrial ecosystem. Lepidoptera are also indicators of environmental change and serve as models for research on mimicry and genetics. They have been central to the development of coevolutionary hypotheses, such as butterflies with flowering plants and moths’ evolutionary arms race with echolocating bats. However, these hypotheses have not been rigorously tested, because a robust lepidopteran phylogeny and timing of evolutionary novelties are lacking. To address these issues, we inferred a comprehensive phylogeny of Lepidoptera, using the largest dataset assembled for the order (2,098 orthologous protein-coding genes from transcriptomes of 186 species, representing nearly all superfamilies), and dated it with carefully evaluated synapomorphy-based fossils. The oldest members of the Lepidoptera crown group appeared in the Late Carboniferous (∼300 Ma) and fed on nonvascular land plants. Lepidoptera evolved the tube-like proboscis in the Middle Triassic (∼241 Ma), which allowed them to acquire nectar from flowering plants. This morphological innovation, along with other traits, likely promoted the extraordinary diversification of superfamily-level lepidopteran crown groups. The ancestor of butterflies was likely nocturnal, and our results indicate that butterflies became day-flying in the Late Cretaceous (∼98 Ma). Moth hearing organs arose multiple times before the evolutionary arms race between moths and bats, perhaps initially detecting a wide range of sound frequencies before being co-opted to specifically detect bat sonar. Our study provides an essential framework for future comparative studies on butterfly and moth evolution.

It has been hypothesized that selection pressure from pathogens plays an important role in shaping social evolution. Social behaviour, in particular brood care, is associated with pathogen pressure in wood-dwelling “lower” termites. Yet, generally pathogen pressure is low in wood-dwelling termite species that never leave the nest except for the mating flight. In comparison, pathogen pressure is higher in species that leave the nest to forage, and thus constantly encounter a diversity of microbes from their environment. We hypothesized that such differences in pathogen pressure are also reflected by differences in the intensity of natural selection on immune genes. We tested this hypothesis in a phylogenetic framework, analyzing rates of non-synonymous and synonymous substitutions on single-copy immune genes. Therefore, we leveraged recent genomic and transcriptomic data from eight termite species, representing wood-dwelling and foraging species as well as 14 additional species spanning the winged insects (Pterygota). Our results provide no evidence for a role of pathogen pressure in selection intensity on single-copy immune genes. Instead, we found evidence for a genome-wide pattern of relaxed selection in termites.Life Science Identifiers (as available Zoobank)Ephemera danica:urn:lsid:zoobank.org:act:06633F75-4809-4BB3-BDCB-6270795368D5Coptotermes sp.urn:lsid:zoobank.org:pub:D6724B7F-F27A-47DC-A4FC-12859ECA0C71Blattella germanica:rn:lsid:zoobank.org:pub:1EA126BA-E9D2-4AA6-8202-26BA5B09B8ADLocusta migratoriaurn:lsid:zoobank.org:pub:D792A09E-844A-412A-BFCA-5293F8388F8CPeriplaneta americana (Blatta americana):urn:lsid:zoobank.org:act:95113A55-4C6D-4DC7-A0E5-620BACADFFE5Apis mellifera:urn:lsid:zoobank.org:act:9082C709-6347-4768-A0DC-27DC44400CB2Bombyx mori (Phalæna (Bombyx) mori)urn:lsid:zoobank.org:act:215466E3-E77F-46E9-8097-372837D7A375Drosophila melanogaster:urn:lsid:zoobank.org:act:5B39F0AA-270D-4AA8-B9A3-C36A3A265910

The beetle superfamily Dytiscoidea, placed within the suborder Adephaga, comprises six families. The phylogenetic relationships of these families, whose species are aquatic, remain highly contentious. In particular the monophyly of the geographically disjunct Aspidytidae (China and South Africa) remains unclear. Here we use a phylogenomic approach to demonstrate that Aspidytidae are indeed monophyletic, as we inferred this phylogenetic relationship from analyzing nucleotide sequence data filtered for compositional heterogeneity and from analyzing amino-acid sequence data. Our analyses suggest that Aspidytidae are the sister group of Amphizoidae, although the support for this relationship is not unequivocal. A sister group relationship of Hygrobiidae to a clade comprising Amphizoidae, Aspidytidae, and Dytiscidae is supported by analyses in which model assumptions are violated the least. In general, we find that both concatenation and the applied coalescent method are sensitive to the effect of among-species compositional heterogeneity. Four-cluster likelihood-mapping suggests that despite the substantial size of the dataset and the use of advanced analytical methods, statistical support is weak for the inferred phylogenetic placement of Hygrobiidae. These results indicate that other kinds of data (e.g. genomic meta-characters) are possibly required to resolve the above-specified persisting phylogenetic uncertainties. Our study illustrates various data-driven confounding effects in phylogenetic reconstructions and highlights the need for careful monitoring of model violations prior to phylogenomic analysis.

The Calyptratae, one of the most species-rich fly clades, only originated and diversified after the Cretaceous-Palaeogene extinction event and yet exhibit high species diversity and a diverse array of life history strategies including predation, phytophagy, saprophagy, haematophagy and parasitism. We present the first phylogenomic analysis of calyptrate relationships. The analysis is based on 40 species representing all calyptrate families and on nucleotide and amino acid data for 1456 single-copy protein-coding genes obtained from shotgun sequencing of transcriptomes. Topologies are overall well resolved, robust and largely congruent across trees obtained with different approaches (maximum parsimony, maximum likelihood, coalescent-based species tree, four-cluster likelihood mapping). Many nodes have 100% bootstrap and jackknife support, but the true support varies by more than one order of magnitude [Bremer support from 3 to 3427; random addition concatenation analysis (RADICAL) gene concatenation size from 10 to 1456]. Analyses of a Dayhoff-6 recoded amino acid dataset also support the robustness of many clades. The backbone topology Hippoboscoidea+(Fanniidae+(Muscidae+((Anthomyiidae-Scathophagidae)+Oestroidea))) is strongly supported and most families are monophyletic (exceptions: Anthomyiidae and Calliphoridae). The monotypic Ulurumyiidae is either alone or together with Mesembrinellidae as the sister group to the rest of Oestroidea. The Sarcophagidae are sister to Mystacinobiidae+Oestridae. Polleniinae emerge as sister group to Tachinidae and the monophyly of the clade Calliphorinae+Luciliinae is well supported, but the phylogenomic data cannot confidently place the remaining blowfly subfamilies (Helicoboscinae, Ameniinae, Chrysomyinae). Compared to hypotheses from the Sanger sequencing era, many clades within the muscoid grade are congruent but now have much higher support. Within much of Oestroidea, Sanger era and phylogenomic data struggle equally with regard to finding well-supported hypotheses.

Phylogenetic relationships among subgroups of cockroaches and termites are still matters of debate. Their divergence times and major phenotypic transitions during evolution are also not yet settled. We addressed these points by combining the first nuclear phylogenomic study of termites and cockroaches with a thorough approach to divergence time analysis, identification of endosymbionts, and reconstruction of ancestral morphological traits and behaviour. Analyses of the phylogenetic relationships within Blattodea robustly confirm previously uncertain hypotheses such as the sister-group relationship between Blaberoidea and remaining Blattodea, and Lamproblatta being the closest relative to the social and wood-feeding Cryptocercus and termites. Consequently, we propose new names for various clades in Blattodea: Cryptocercus + termites = Tutricablattae; Lamproblattidae + Tutricablattae = Kittrickea; and Blattoidea + Corydioidea = Solumblattodea. Our inferred divergence times contradict previous studies by showing that most subgroups of Blattodea evolved in the Cretaceous, reducing the gap between molecular estimates of divergence times and the fossil record. On a phenotypic level, the blattodean ground-plan is for egg packages to be laid directly in a hole while other forms of oviposition, including ovovivipary and vivipary, arose later. Finally, other changes in egg care strategy may have allowed for the adaptation of nest building and other novelties.

Social insects provide promising new avenues for aging research. Within a colony, individuals that share the same genetic background can differ in lifespan by up to two orders of magnitude. Reproducing queens (and in termites also kings) can live for more than 20 years, extraordinary lifespans for insects. We studied aging in a termite species, Cryptotermes secundus, which lives in less socially complex societies with a few hundred colony members. Reproductives develop from workers which are totipotent immatures. Comparing transcriptomes of young and old individuals, we found evidence for aging in reproductives that was especially associated with DNA and protein damage and the activity of transposable elements. By contrast, workers seemed to be better protected against aging. Thus our results differed from those obtained for social insects that live in more complex societies. Yet, they are in agreement with lifespan estimates for the study species. Our data are also in line with expectations from evolutionary theory. For individuals that are able to reproduce, it predicts that aging should only start after reaching maturity. As C. secundus workers are immatures with full reproductive options we expect them to invest into anti-aging processes. Our study illustrates that the degree of aging can differ between social insects and that it may be associated with caste-specific opportunities for reproduction.

Termites are important ecosystem engineers. Yet they are often difficult to identify due to the lack of reliable species-specific morphological traits for many species, which hampers ecological research. Recently, termitologists working with West African termites (West African Termite Taxonomy Initiative) convened for a workshop with the aim of beginning to address this problem. Repeated determination of the same termite samples by the most renowned taxonomists for West African termites identified the huge scale of the problem, as less than 10% of all species could be unambiguously determined to the species level. Intensive discussions and comparisons increased the identification success to around 25% at the end of the workshop. Yet many groups remained problematic and molecular markers and barcoding techniques combined with species delimitation approaches will be needed to help resolve these existing taxonomic problems. Based on the outcome of this workshop, we propose concerted initiatives to address termite taxonomy on a global scale. We are convinced that dedicated workshops on regional taxonomy that follow a similar structured approach, with repeated determination of the same sample, will help overcome the difficulties that termite taxonomy faces. This initiative can also serve as a blueprint for other taxonomical groups that are difficult to identify., SCOPUS: re.j, info:eu-repo/semantics/published

Hemipteroid insects (Paraneoptera), with over 10% of all known insect diversity, are a major component of terrestrial and aquatic ecosystems. Previous phylogenetic analyses have not consistently resolved the relationships among major hemipteroid lineages. We provide maximum likelihood-based phylogenomic analyses of a taxonomically comprehensive dataset comprising sequences of 2,395 single-copy, protein-coding genes for 193 samples of hemipteroid insects and outgroups. These analyses yield a well-supported phylogeny for hemipteroid insects. Monophyly of each of the three hemipteroid orders (Psocodea, Thysanoptera, and Hemiptera) is strongly supported, as are most relationships among suborders and families. Thysanoptera (thrips) is strongly supported as sister to Hemiptera. However, as in a recent large-scale analysis sampling all insect orders, trees from our data matrices support Psocodea (bark lice and parasitic lice) as the sister group to the holometabolous insects (those with complete metamorphosis). In contrast, four-cluster likelihood mapping of these data does not support this result. A molecular dating analysis using 23 fossil calibration points suggests hemipteroid insects began diversifying before the Carboniferous, over 365 million years ago. We also explore implications for understanding the timing of diversification, the evolution of morphological traits, and the evolution of mitochondrial genome organization. These results provide a phylogenetic framework for future studies of the group.

It has been experimentally shown that DNA methylation is involved in the regulation of gene expression and the silencing of transposable element activity in eukaryotes. The variable levels of DNA methylation among different insect species indicate an evolutionarily flexible role of DNA methylation in insects, which due to a lack of comparative data is not yet well-substantiated. Here, we use computational methods to trace signatures of DNA methylation across insects by analyzing transcriptomic and genomic sequence data from all currently recognized insect orders. We conclude that: 1) a functional methylation system relying exclusively on DNA methyltransferase 1 is widespread across insects. 2) DNA methylation has potentially been lost or extremely reduced in species belonging to springtails (Collembola), flies and relatives (Diptera), and twisted-winged parasites (Strepsiptera). 3) Holometabolous insects display signs of reduced DNA methylation levels in protein-coding sequences compared with hemimetabolous insects. 4) Evolutionarily conserved insect genes associated with housekeeping functions tend to display signs of heavier DNA methylation in comparison to the genomic/transcriptomic background. With this comparative study, we provide the much needed basis for experimental and detailed comparative analyses required to gain a deeper understanding on the evolution and function of DNA methylation in insects.

The phylogenetic relationships of the winged insect lineages - mayflies (Ephemeroptera), damselflies and dragonflies (Odonata), and all other winged insects (Neoptera) - are still controversial with three hypotheses supported by different datasets: Palaeoptera, Metapterygota and Chiastomyaria. Here, we reanalyze available phylogenomic data with a focus on detecting confounding and alternative signal. In this context, we provide a framework to quantitatively evaluate and assess incongruent molecular phylogenetic signal inherent in phylogenomic datasets. Despite overall support for the Palaeoptera hypothesis, we also found considerable signal for Chiastomyaria, which is not easily detectable by standardized tree inference approaches. Analyses of the accumulation of signal across gene partitions showed that signal accumulates gradually. However, even in case signal only slightly supported one over the other hypothesis, topologies inferred from large datasets switch from statistically strongly supported Palaeoptera to strongly supported Chiastomyaria. From a morphological point of view, Palaeoptera currently appears to be the best-supported hypothesis; however, recent analyses were restricted to head characters. Phylogenetic approaches covering all organ systems including analyses of potential functional or developmental convergence are still pending so that the Palaeoptera problem has to be considered an open question in insect systematics.

In the scorpionfly Panorpa, a recent study suggested monochromatic vision due to evidence of only a single opsin found in transcriptome data. To reconsider this hypothesis, the present study investigates opsin expression using transcriptome data of 21 species including representatives of all major lineages of scorpionflies (Mecoptera) and of three families of their closest relatives, the fleas (Siphonaptera). In most mecopteran species investigated, transcripts encode two opsins with predicted peak absorbances in the green, two in the blue, and one in the ultraviolet spectral region. Only in groups with reduced or absent ocelli, like Caurinus and Apteropanorpa, less than four visual opsin messenger RNAs have been identified. In addition, we found a Rh7-like opsin in transcriptome data derived from larvae of the mecopteran Nannochorista, and in two flea species. Peropsin expression was observed in two mecopterans. In light of these new data, we question the hypothesis on monochromatic vision in the genus Panorpa. In a broader phylogenetic perspective, it is suggested that the common ancestor of the monophyletic taxon Antliophora (Diptera, Mecoptera and Siphonaptera) possessed the full set of visual opsins, a Rh7-like opsin, and in addition a pteropsin as well as a peropsin. In the course of evolution individual opsins were likely lost in several lineages of this clade.

Background The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. Results Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. Conclusions The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda. Electronic supplementary material The online version of this article (10.1186/s12864-017-4197-1) contains supplementary material, which is available to authorized users.

Brochosomes (BS) are secretory granules resembling buckyballs, produced intracellularly in specialized glandular segments of the Malpighian tubules and forming superhydrophobic coatings on the integuments of leafhoppers (Hemiptera, Cicadellidae). Their composition is poorly known. Using a combination of SDS-PAGE, LC-MS/MS, next-generation sequencing (RNAseq) and bioinformatics we demonstrate that the major structural component of BS of the leafhopper Graphocephala fennahi Young is a novel family of 21–40-kDa secretory proteins, referred to herein as brochosomins (BSM), apparently cross-linked by disulfide bonds. At least 28 paralogous BSM were identified in a transcriptome assembly of this species, most of which were detected in BS. Multiple additional BS-associated proteins (BSAP), possibly loosely attached to the outer and inner surfaces of BS, were also identified; some of these were glycine-, tyrosine- and proline-rich. BSM and BSAP together accounted for half of the 100 most expressed transcripts in the Malpighian tubules of G. fennahi. Except for several minor BSAP possibly related to cyclases, BSM and BSAP had no homologs among known proteins, thus representing taxonomically restricted gene families (orphans). Searching in 50 whole-body transcriptome assemblies of Hemiptera found homologs of BSM in representatives of all five families of the superfamily Membracoidea (Cicadellidae, Myerslopiidae, Aetalionidae, Membracidae, and Melizoderidae), but not in other lineages. Among the identified proteins only BSM were shared in common between all 17 surveyed leafhoppers known to produce BS. Combined CHN elemental and aminoacid analyses estimated the total protein content of BS from the integument of G. fennahi to be 60–70%.