Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

A winged relative of ice-crawlers in amber bridges the cryptic extant Xenonomia and a rich fossil record

Until the advent of phylogenomics, the atypical morphology of extant representatives of the insect orders Grylloblattodea (ice-crawlers) and Mantophasmatodea (gladiators) had confounding effects on efforts to resolve their placement within Polyneoptera. This recent research has unequivocally shown that these species-poor groups are closely related and form the clade Xenonomia. Nonetheless, divergence dates of these groups remain poorly constrained, and their evolutionary history debated, as the few well-identified fossils, characterized by a suite of morphological features similar to that of extant forms, are comparatively young. Notably, the extant forms of both groups are wingless, whereas most of the pre-Cretaceous insect fossil record is composed of winged insects, which represents a major shortcoming of the taxonomy. Here, we present new specimens embedded in mid-Cretaceous amber from Myanmar and belonging to the recently described species Aristovia daniili. The abundant material and pristine preservation allowed a detailed documentation of the morphology of the species, including critical head features. Combined with a morphological data set encompassing all Polyneoptera, these new data unequivocally demonstrate that A. daniili is a winged stem Grylloblattodea. This discovery demonstrates that winglessness was acquired independently in Grylloblattodea and Mantophasmatodea. Concurrently, wing apomorphic traits shared by the new fossil and earlier fossils demonstrate that a large subset of the former "Protorthoptera" assemblage, representing a third of all known insect species in some Permian localities, are genuine representatives of Xenonomia. Data from the fossil record depict a distinctive evolutionary trajectory, with the group being both highly diverse and abundant during the Permian but experiencing a severe decline from the Triassic onwards.

The Chalcidoidea bush of life: evolutionary history of a massive radiation of minute wasps

Chalcidoidea are mostly parasitoid wasps that include as many as 500 000 estimated species. Capturing phylogenetic signal from such a massive radiation can be daunting. Chalcidoidea is an excellent example of a hyperdiverse group that has remained recalcitrant to phylogenetic resolution. We combined 1007 exons obtained with Anchored Hybrid Enrichment with 1048 ultra-conserved elements (UCEs) for 433 taxa including all extant families, >95% of all subfamilies, and 356 genera chosen to represent the vast diversity of the superfamily. Going back and forth between the molecular results and our collective knowledge of morphology and biology, we detected bias in the analyses that was driven by the saturation of nucleotide data. Our final results are based on a concatenated analysis of the least saturated exons and UCE datasets (2054 loci, 284 106 sites). Our analyses support an expected sister relationship with Mymarommatoidea. Seven previously recognized families were not monophyletic, so support for a new classification is discussed. Natural history in some cases would appear to be more informative than morphology, as illustrated by the elucidation of a clade of plant gall associates and a clade of taxa with planidial first-instar larvae. The phylogeny suggests a transition from smaller soft-bodied wasps to larger and more heavily sclerotized wasps, with egg parasitism as potentially ancestral for the entire superfamily. Deep divergences in Chalcidoidea coincide with an increase in insect families in the fossil record, and an early shift to phytophagy corresponds with the beginning of the "Angiosperm Terrestrial Revolution". Our dating analyses suggest a middle Jurassic origin of 174 Ma (167.3-180.5 Ma) and a crown age of 162.2 Ma (153.9-169.8 Ma) for Chalcidoidea. During the Cretaceous, Chalcidoidea may have undergone a rapid radiation in southern Gondwana with subsequent dispersals to the Northern Hemisphere. This scenario is discussed with regard to knowledge about the host taxa of chalcid wasps, their fossil record and Earth's palaeogeographic history.

Bridging two insect flight modes in evolution, physiology and robophysics

Since taking flight, insects have undergone repeated evolutionary transitions between two seemingly distinct flight modes1-3. Some insects neurally activate their muscles synchronously with each wingstroke. However, many insects have achieved wingbeat frequencies beyond the speed limit of typical neuromuscular systems by evolving flight muscles that are asynchronous with neural activation and activate in response to mechanical stretch2-8. These modes reflect the two fundamental ways of generating rhythmic movement: time-periodic forcing versus emergent oscillations from self-excitation8-10. How repeated evolutionary transitions have occurred and what governs the switching between these distinct modes remain unknown. Here we find that, despite widespread asynchronous actuation in insects across the phylogeny3,6, asynchrony probably evolved only once at the order level, with many reversions to the ancestral, synchronous mode. A synchronous moth species, evolved from an asynchronous ancestor, still preserves the stretch-activated muscle physiology. Numerical and robophysical analyses of a unified biophysical framework reveal that rather than a dichotomy, these two modes are two regimes of the same dynamics. Insects can transition between flight modes across a bridge in physiological parameter space. Finally, we integrate these two actuation modes into an insect-scale robot11-13 that enables transitions between modes and unlocks a new self-excited wingstroke strategy for engineered flight. Together, this framework accounts for repeated transitions in insect flight evolution and shows how flight modes can flip with changes in physiological parameters.

Illusion of flight? Absence, evidence and the age of winged insects

The earliest fossils of winged insects (Pterygota) are mid-Carboniferous (latest Mississippian, 328–324 Mya), but estimates of their age based on fossil-calibrated molecular phylogenetic studies place their origin at 440–370 Mya during the Silurian or Devonian. This discrepancy would require that winged insects evaded fossilization for at least the first ~50 Myr of their history. Here, we examine the plausibility of such a gap in the fossil record, and possible explanations for it, based on comparisons with the fossil records of other arthropod groups, the distribution of first occurrence dates of pterygote families, phylogenetically informed simulations of the fossilization of Palaeozoic insects, and re-analysis of data presented by Misof and colleagues using updated fossil calibrations under a variety of prior probability settings. We do not find support for the mechanisms previously suggested to account for such an extended gap in the pterygote fossil record, including sampling bias, preservation bias, and body size. We suggest that inference of an early origin of Pterygota long prior to their first appearance in the fossil record is probably an analytical artefact of taxon sampling and choice of fossil calibration points, possibly compounded by heterogeneity in rates of sequence evolution or speciation, including radiations or ‘bursts’ during their early history.

Early diversifications of angiosperms and their insect pollinators: were they unlinked?

The present-day ubiquity of angiosperm-insect pollination has led to the hypothesis that these two groups coevolved early in their evolutionary history. However, recent fossil discoveries and fossil-calibrated molecular dating analyses challenge the notion that early diversifications of angiosperms and insects were inextricably linked. In this article, we examine (i) the discrepancies between dates of emergence for angiosperms and major clades of insects; (ii) the long history of gymnosperm-insect pollination modes, which likely shaped early angiosperm-insect pollination mutualisms; and (iii) how the K-Pg (Cretaceous-Paleogene) mass extinction event was vital in propelling modern angiosperm-insect mutualisms. We posit that the early diversifications of angiosperms and their insect pollinators were largely decoupled until the end of the Cretaceous.

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

Paradoxes of Hymenoptera flight muscles, extreme machines

In the Carboniferous, insects evolved flight. Intense selection drove for high performance and approximately 100 million years later, Hymenoptera (bees, wasps and ants) emerged. Some species had proportionately small wings, with apparently impossible aerodynamic challenges including a need for high frequency flight muscles (FMs), powered exclusively off aerobic pathways and resulting in extreme aerobic capacities. Modern insect FMs are the most refined and form large dense blocks that occupy 90% of the thorax. These can beat wings at 200 to 230 Hz, more than double that achieved by standard neuromuscular systems. To do so, rapid repolarisation was circumvented through evolution of asynchronous stimulation, stretch activation, elastic recoil and a paradoxically slow Ca2+ reuptake. While the latter conserves ATP, considerable ATP is demanded at the myofibrils. FMs have diminished sarcoplasmic volumes, and ATP is produced solely by mitochondria, which pack myocytes to maximal limits and have very dense cristae. Gaseous oxygen is supplied directly to mitochondria. While FMs appear to be optimised for function, several unusual paradoxes remain. FMs lack any significant equivalent to the creatine kinase shuttle, and myofibrils are twice as wide as those of within cardiomyocytes. The mitochondrial electron transport systems also release large amounts of reactive oxygen species (ROS) and respiratory complexes do not appear to be present at any exceptional level. Given that the loss of the creatine kinase shuttle and elevated ROS impairs heart function, we question how do FM shuttle adenylates at high rates and tolerate oxidative stress conditions that occur in diseased hearts?

Phylogenomic Analysis of Velvet Worms (Onychophora) Uncovers an Evolutionary Radiation in the Neotropics

Onychophora ("velvet worms") are charismatic soil invertebrates known for their status as a "living fossil," their phylogenetic affiliation to arthropods, and their distinctive biogeographic patterns. However, several aspects of their internal phylogenetic relationships remain unresolved, limiting our understanding of the group's evolutionary history, particularly with regard to changes in reproductive mode and dispersal ability. To address these gaps, we used RNA sequencing and phylogenomic analysis of transcriptomes to reconstruct the evolutionary relationships and infer divergence times within the phylum. We recovered a fully resolved and well-supported phylogeny for the circum-Antarctic family Peripatopsidae, which retains signals of Gondwanan vicariance and showcases the evolutionary lability of reproductive mode in the family. Within the Neotropical clade of Peripatidae, though, we found that amino acid-translated sequence data masked nearly all phylogenetic signal, resulting in highly unstable and poorly supported relationships. Analyses using nucleotide sequence data were able to resolve many more relationships, though we still saw discordant phylogenetic signal between genes, probably indicative of a rapid, mid-Cretaceous radiation in the group. Finally, we hypothesize that the unique reproductive mode of placentotrophic viviparity found in all Neotropical peripatids may have facilitated the multiple inferred instances of over-water dispersal and establishment on oceanic islands.

Spatial distribution of campaniform sensilla mechanosensors on wings: form, function, and phylogeny

Insect wings serve two crucial functions in flight: propulsion and sensing. During flapping flight, complex spatiotemporal patterns of strain on the wing reflect mechanics, kinematics, and external perturbations; sensing wing deformation provides feedback necessary for flight control. Campaniform sensilla distributed across the wing transduce local strain fluctuations into neural signals, so their placement on the wing determines sensory information available to the insect. Thus, understanding the significance of these sensor locations will also reveal how sensing and wing movement are coupled. Here, we identify trends in wing campaniform sensilla placement across flying insects from the literature. We then discuss how these patterns can influence sensory encoding by wing mechanosensors. Finally, we propose combining a comparative approach on model insect clades with computational modeling, leveraging the spectacular natural diversity in wings to uncover biological principles of mechanosensory feedback in flight control.

Rewinding the molecular clock in the genus Carabus (Coleoptera: Carabidae) in light of fossil evidence and the Gondwana split: A reanalysis

Molecular clocks have become powerful tools given increasing sequencing and fossil resources. However, calibration analyses outcomes depend on the choice of priors. Here, we revisited the seminal dating study published by Andújar and coworkers of the genus Carabus proposing that prior choices need re-evaluation. We hypothesized that reflecting fossil evidence and the Gondwanan split properly significantly rewinds the molecular clock. We re-used the dataset including five mitochondrial and four nuclear DNA fragments with a total length of 7888 nt. Fossil evidence for Oligocene occurrence of Calosoma was considered. Root age was set based on the fossil evidence of Harpalinae ground beetles in the Upper Cretaceous. Paleogene divergence of the outgroup taxa Ceroglossini and Pamborini is introduced as a new prior based on current paleontological and geological literature. The ultrametric time-calibrated tree of the extended nd5 dataset resulted in a median TMRCA Carabus of 53.92 Ma (HPD 95% 45.01–63.18 Ma), roughly 30 Ma older than in the Andújar study. The splits among C. rugosus and C. morbillosus (A), C. riffensis from the European Mesocarabus (B), and Eurycarabus and Nesaeocarabus (C) were dated to 17.58 (12.87–22.85), 24.14 (18.02–30.58), and 21.6 (16.44–27.43) Ma. They were decidedly older than those previously reported (7.48, 10.93, and 9.51 Ma). These changes were driven almost entirely by constraining the Carabidae time-tree root with a Harpalinae amber fossil at ~99 Ma. Utilizing the nd5 dating results of three well-supported Carabus clades as secondary calibration points for the complete MIT-NUC dataset led to a TMRCA of Carabus of 44.72 (37.54–52.22) Ma, compared with 25.16 Ma (18.41–33.04 Ma) in the previous study. Considering fossil evidence for Oligocene Calosoma and Late Cretaceous Harpalini together with the Gondwanan split as a new prior, our new approach supports the origin of genus Carabus in the Eocene. Our results are preliminary because of the heavy reliance on the nd5 gene, and thus will have to be tested with a sufficient set of nuclear markers. Additionally, uncertainties due to dating root age of the tree based on a single fossil and outgroup taxon affect the results. Improvement of the fossil database, particularly in the supertribe Carabitae, is needed to reduce these uncertainties in dating Carabus phylogeny.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.

Survey of biomechanical aspects of arthropod terrestrialisation - Substrate bound legged locomotion

Arthropods are the most diverse clade on earth with regard to both species number and variability of body plans. Their general body plan is characterised by variable numbers of legs, and many-legged locomotion is an essential aspect of many aquatic and terrestrial arthropod species. Moreover, arthropods belong to the first groups of animals to colonise subaerial habitats, and they did so repeatedly and independently in a couple of clades. Those arthropod clades that colonised land habitats were equipped with highly variable body plans and locomotor apparatuses. Proceeding from their respective specific anatomies, they were challenged with strongly changing environmental conditions as well as altered physical and physiological constraints. This review explores the transitions from aquatic to terrestrial habitats across the different arthropod body plans and explains the major mechanisms and principles that constrain design and function of a range of locomotor apparatuses. Important aspects of movement physiology addressed here include the effects of different numbers of legs, different body sizes, miniaturisation and simplification of body plans and different ratios of inertial and damping forces. The article's focus is on continuous legged locomotion, but related ecological and behavioural aspects are also taken into account.

Mind the Outgroup and Bare Branches in Total-Evidence Dating: a Case Study of Pimpliform Darwin Wasps (Hymenoptera, Ichneumonidae)

Taxon sampling is a central aspect of phylogenetic study design, but it has received limited attention in the context of total-evidence dating, a widely used dating approach that directly integrates molecular and morphological information from extant and fossil taxa. We here assess the impact of commonly employed outgroup sampling schemes and missing morphological data in extant taxa on age estimates in a total-evidence dating analysis under the uniform tree prior. Our study group is Pimpliformes, a highly diverse, rapidly radiating group of parasitoid wasps of the family Ichneumonidae. We analyze a data set comprising 201 extant and 79 fossil taxa, including the oldest fossils of the family from the Early Cretaceous and the first unequivocal representatives of extant subfamilies from the mid-Paleogene. Based on newly compiled molecular data from ten nuclear genes and a morphological matrix that includes 222 characters, we show that age estimates become both older and less precise with the inclusion of more distant and more poorly sampled outgroups. These outgroups not only lack morphological and temporal information but also sit on long terminal branches and considerably increase the evolutionary rate heterogeneity. In addition, we discover an artifact that might be detrimental for total-evidence dating: “bare-branch attraction,” namely high attachment probabilities of certain fossils to terminal branches for which morphological data are missing. Using computer simulations, we confirm the generality of this phenomenon and show that a large phylogenetic distance to any of the extant taxa, rather than just older age, increases the risk of a fossil being misplaced due to bare-branch attraction. After restricting outgroup sampling and adding morphological data for the previously attracting, bare branches, we recover a Jurassic origin for Pimpliformes and Ichneumonidae. This first age estimate for the group not only suggests an older origin than previously thought but also that diversification of the crown group happened well before the Cretaceous-Paleogene boundary. Our case study demonstrates that in order to obtain robust age estimates, total-evidence dating studies need to be based on a thorough and balanced sampling of both extant and fossil taxa, with the aim of minimizing evolutionary rate heterogeneity and missing morphological information. [Bare-branch attraction; ichneumonids; fossils; morphological matrix; phylogeny; RoguePlots.]

Phylomitogenomics provides new perspectives on the Euphasmatodea radiation (Insecta: Phasmatodea)

Phasmatodea species diversity lies almost entirely within its suborder Euphasmatodea, which exhibits a pantropical distribution and is considered to derive from a recent and rapid evolutionary radiation. To shed light on Euphasmatodea origins and diversification, we assembled the mitogenomes of 17 species from transcriptomic sequencing data and analysed them along with 22 already available Phasmatodea mitogenomes and 33 mitogenomes representing most of the Polyneoptera lineages. Maximum Likelihood and Bayesian Inference approaches retrieved consistent topologies, both showing the widespread conflict between phylogenetic approaches and traditional systematics. We performed a divergence time analysis leveraging ten fossil specimens representative of most polyneopteran lineages: the time tree obtained supports an older radiation of the clade with respect to previous hypotheses. Euphasmatodea diversification is inferred to have started ~ 187 million years ago, suggesting that the Triassic-Jurassic mass extinction and the breakup of Pangea could have contributed to the process. We also investigated Euphasmatodea mitogenomes patterns of dN, dS and dN/dS ratio throughout our time-tree, trying to characterize the selective regime which may have shaped the clade evolution.

An integrative phylogenomic approach to elucidate the evolutionary history and divergence times of Neuropterida (Insecta: Holometabola)

BACKGROUND

The latest advancements in DNA sequencing technologies have facilitated the resolution of the phylogeny of insects, yet parts of the tree of Holometabola remain unresolved. The phylogeny of Neuropterida has been extensively studied, but no strong consensus exists concerning the phylogenetic relationships within the order Neuroptera. Here, we assembled a novel transcriptomic dataset to address previously unresolved issues in the phylogeny of Neuropterida and to infer divergence times within the group. We tested the robustness of our phylogenetic estimates by comparing summary coalescent and concatenation-based phylogenetic approaches and by employing different quartet-based measures of phylogenomic incongruence, combined with data permutations.

RESULTS

Our results suggest that the order Raphidioptera is sister to Neuroptera + Megaloptera. Coniopterygidae is inferred as sister to all remaining neuropteran families suggesting that larval cryptonephry could be a ground plan feature of Neuroptera. A clade that includes Nevrorthidae, Osmylidae, and Sisyridae (i.e. Osmyloidea) is inferred as sister to all other Neuroptera except Coniopterygidae, and Dilaridae is placed as sister to all remaining neuropteran families. Ithonidae is inferred as the sister group of monophyletic Myrmeleontiformia. The phylogenetic affinities of Chrysopidae and Hemerobiidae were dependent on the data type analyzed, and quartet-based analyses showed only weak support for the placement of Hemerobiidae as sister to Ithonidae + Myrmeleontiformia. Our molecular dating analyses suggest that most families of Neuropterida started to diversify in the Jurassic and our ancestral character state reconstructions suggest a primarily terrestrial environment of the larvae of Neuropterida and Neuroptera.

CONCLUSION

Our extensive phylogenomic analyses consolidate several key aspects in the backbone phylogeny of Neuropterida, such as the basal placement of Coniopterygidae within Neuroptera and the monophyly of Osmyloidea. Furthermore, they provide new insights into the timing of diversification of Neuropterida. Despite the vast amount of analyzed molecular data, we found that certain nodes in the tree of Neuroptera are not robustly resolved. Therefore, we emphasize the importance of integrating the results of morphological analyses with those of sequence-based phylogenomics. We also suggest that comparative analyses of genomic meta-characters should be incorporated into future phylogenomic studies of Neuropterida.

Reliable Confidence Intervals for RelTime Estimates of Evolutionary Divergence Times

Confidence intervals (CIs) depict the statistical uncertainty surrounding evolutionary divergence time estimates. They capture variance contributed by the finite number of sequences and sites used in the alignment, deviations of evolutionary rates from a strict molecular clock in a phylogeny, and uncertainty associated with clock calibrations. Reliable tests of biological hypotheses demand reliable CIs. However, current non-Bayesian methods may produce unreliable CIs because they do not incorporate rate variation among lineages and interactions among clock calibrations properly. Here, we present a new analytical method to calculate CIs of divergence times estimated using the RelTime method, along with an approach to utilize multiple calibration uncertainty densities in dating analyses. Empirical data analyses showed that the new methods produce CIs that overlap with Bayesian highest posterior density intervals. In the analysis of computer-simulated data, we found that RelTime CIs show excellent average coverage probabilities, that is, the actual time is contained within the CIs with a 94% probability. These developments will encourage broader use of computationally efficient RelTime approaches in molecular dating analyses and biological hypothesis testing.

Identification of Y-chromosome scaffolds of the Queensland fruit fly reveals a duplicated gyf gene paralogue common to many Bactrocera pest species

Bactrocera tryoni (Queensland fruit fly) are polyphagous horticultural pests of eastern Australia. Heterogametic males contain a sex-determining Y-chromosome thought to be gene poor and repetitive. Here, we report 39 Y-chromosome scaffolds (~700 kb) from B. tryoni identified using genotype-by-sequencing data and whole-genome resequencing. Male diagnostic PCR assays validated eight Y-scaffolds, and one (Btry4096) contained a novel gene with five exons that encode a predicted 575 amino acid protein. The Y-gene, referred to as typo-gyf, is a truncated Y-chromosome paralogue of X-chromosome gene gyf (1773 aa). The Y-chromosome contained ~41 copies of typo-gyf, and expression occurred in male flies and embryos. Analysis of 13 tephritid transcriptomes confirmed typo-gyf expression in six additional Bactrocera species, including Bactrocera latifrons, Bactrocera dorsalis and Bactrocera zonata. Molecular dating estimated typo-gyf evolved within the past 8.02 million years (95% highest posterior density 10.56-5.52 million years), after the split with Bactrocera oleae. Phylogenetic analysis also highlighted complex evolutionary histories among several Bactrocera species, as discordant nuclear (116 genes) and mitochondrial (13 genes) topologies were observed. B. tryoni Y-sequences may provide useful sites for future transgene insertions, and typo-gyf could act as a Y-chromosome diagnostic marker for many Bactrocera species, although its function is unknown.

Recalibration of the insect evolutionary time scale using Monte San Giorgio fossils suggests survival of key lineages through the End-Permian Extinction

Insects are a highly diverse group of organisms and constitute more than half of all known animal species. They have evolved an extraordinary range of traits, from flight and complete metamorphosis to complex polyphenisms and advanced eusociality. Although the rich insect fossil record has helped to chart the appearance of many phenotypic innovations, data are scarce for a number of key periods. One such period is that following the End-Permian Extinction, recognized as the most catastrophic of all extinction events. We recently discovered several 240-million-year-old insect fossils in the Mount San Giorgio Lagerstätte (Switzerland-Italy) that are remarkable for their state of preservation (including internal organs and soft tissues), and because they extend the records of their respective taxa by up to 200 million years. By using these fossils as calibrations in a phylogenomic dating analysis, we present a revised time scale for insect evolution. Our date estimates for several major lineages, including the hyperdiverse crown groups of Lepidoptera, Hemiptera: Heteroptera and Diptera, are substantially older than their currently accepted post-Permian origins. We found that major evolutionary innovations, including flight and metamorphosis, appeared considerably earlier than previously thought. These results have numerous implications for understanding the evolution of insects and their resilience in the face of extreme events such as the End-Permian Extinction.

Exploring the underwater silken architectures of caddisworms: comparative silkomics across two caddisfly suborders

Caddisfly (Trichoptera) larvae assemble a variety of underwater structures using bioadhesive silk. The order is divided into two primary sub-orders distinguished by how the larvae deploy their silk. Foraging Integripalpia larvae construct portable tube cases. Annulipalpia larvae construct stationary retreats, some with suspended nets to capture food. To identify silk molecular adaptations that may have contributed to caddisfly diversification, we report initial characterization of silk from a net-spinner genus, Parapsyche, for comparison with the silk of a tube case-maker genus, Hesperophylax. Overall, general features of silk structure and processing are conserved across the sub-orders despite approximately 200 Ma of divergence: the H-fibroin proteins comprise repeating phosphoserine-rich motifs, naturally spun silk fibres contain approximately 1 : 1 molar ratios of divalent metal ions to phosphate, silk fibre precursors are stored as complex fluids of at least two types of complexes, and silk gland proteins contain only traces of divalent metal ions, suggesting metal ions that solidify the fibres are absorbed from the aqueous environment after silk extrusion. However, the number and arrangement of the repeating phosphoserine blocks differ between genera, suggesting molecular adaptation of H-fibroin through duplication and shuffling of conserved structural modules may correspond with the radiation of caddisflies into diverse environments. This article is part of the theme issue 'Transdisciplinary approaches to the study of adhesion and adhesives in biological systems'.

Simultaneous radiation of bird and mammal lice following the K-Pg boundary

The diversification of parasite groups often occurs at the same time as the diversification of their hosts. However, most studies demonstrating this concordance only examine single host-parasite groups. Multiple diverse lineages of ectoparasitic lice occur across both birds and mammals. Here, we describe the evolutionary history of lice based on analyses of 1107 single-copy orthologous genes from sequenced genomes of 46 species of lice. We identify three major diverse groups of lice: one exclusively on mammals, one almost exclusively on birds and one on both birds and mammals. Each of these groups radiated just after the Cretaceous-Paleogene (K-Pg) boundary, the time of the mass extinction event of the dinosaurs and rapid diversification of most of the modern lineages of birds and mammals.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

Phylogenomics and the evolution of hemipteroid insects

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.

Knockdown of RNA interference pathway genes impacts the fitness of western corn rootworm

Western corn rootworm (Diabrotica virgifera virgifera) is a serious agricultural pest known for its high adaptability to various management strategies, giving rise to a continual need for new control options. Transgenic maize expressing insecticidal RNAs represents a novel mode of action for rootworm management that is dependent on the RNA interference (RNAi) pathways of the insect for efficacy. Preliminary evidence suggests that western corn rootworm could develop broad resistance to all insecticidal RNAs through changes in RNAi pathway genes; however, the likelihood of field-evolved resistance occurring through this mechanism remains unclear. In the current study, eight key genes involved in facilitating interference in the microRNA and small interfering RNA pathways were targeted for knockdown in order to evaluate impact on fitness of western corn rootworm. These genes include drosha, dicer-1, dicer-2, pasha, loquacious, r2d2, argonaute 1, and argonaute 2. Depletion of targeted transcripts in rootworm larvae led to changes in microRNA expression, decreased ability to pupate, reduced adult beetle emergence, and diminished reproductive capacity. The observed effects do not support evolution of resistance through changes in expression of these eight genes due to reduced insect fitness.

A revised phylogeny of macropathine cave crickets (Orthoptera: Rhaphidophoridae) uncovers a paraphyletic Australian fauna

Australian cave crickets are members of the subfamily Macropathinae (Orthoptera: Rhaphidophoridae). The subfamily is thought to have originated prior to the tectonic separation of the supercontinent Gondwana based on distributions of extant lineages and molecular phylogenetic evidence, although the Australian fauna have been underrepresented in previous studies. The current study augments existing multigene data (using 12S, 16S, and 28S rRNA genes) to investigate the placement of the Australian representatives within the Macropathinae and to assess divergence dates of select clades. Results suggest that the endemic Tasmanian genus Parvotettix is the sister lineage to the remaining members of the subfamily, an outcome that presents a paraphyletic Australian fauna in contrast to previous studies. All other Australian taxa represented in this study (Micropathus and Novotettix) emerged as a sister group to the New Zealand and South American macropathine lineages. Estimation of phylogenetic divergence ages among the aforementioned clades were calibrated using two methods, in absence of suitable fossil records: (i) tectonic events depicting the fragmentation of Gondwanan landmasses that invoke vicariant scenarios of present day geographic distributions; and (ii) molecular evolutionary rates. Geological calibrations place the median age of the most recent common ancestor of extant macropathines at ∼125 to ∼165 Ma, whereas analyses derived from molecular substitution rates suggest a considerably younger origin of ∼32 Ma. This phylogenetic study represents the most rigorous taxonomic sampling of the Australian cave cricket fauna to date and stresses the influence of lineage representation on biogeographic inference.

The Past Sure is Tense: On Interpreting Phylogenetic Divergence Time Estimates

Divergence time estimation-the calibration of a phylogeny to geological time-is an integral first step in modeling the tempo of biological evolution (traits and lineages). However, despite increasingly sophisticated methods to infer divergence times from molecular genetic sequences, the estimated age of many nodes across the tree of life contrast significantly and consistently with timeframes conveyed by the fossil record. This is perhaps best exemplified by crown angiosperms, where molecular clock (Triassic) estimates predate the oldest (Early Cretaceous) undisputed angiosperm fossils by tens of millions of years or more. While the incompleteness of the fossil record is a common concern, issues of data limitation and model inadequacy are viable (if underexplored) alternative explanations. In this vein, Beaulieu et al. (2015) convincingly demonstrated how methods of divergence time inference can be misled by both (i) extreme state-dependent molecular substitution rate heterogeneity and (ii) biased sampling of representative major lineages. These results demonstrate the impact of (potentially common) model violations. Here, we suggest another potential challenge: that the configuration of the statistical inference problem (i.e., the parameters, their relationships, and associated priors) alone may preclude the reconstruction of the paleontological timeframe for the crown age of angiosperms. We demonstrate, through sampling from the joint prior (formed by combining the tree (diversification) prior with the calibration densities specified for fossil-calibrated nodes) that with no data present at all, that an Early Cretaceous crown angiosperms is rejected (i.e., has essentially zero probability). More worrisome, however, is that for the 24 nodes calibrated by fossils, almost all have indistinguishable marginal prior and posterior age distributions when employing routine lognormal fossil calibration priors. These results indicate that there is inadequate information in the data to over-rule the joint prior. Given that these calibrated nodes are strategically placed in disparate regions of the tree, they act to anchor the tree scaffold, and so the posterior inference for the tree as a whole is largely determined by the pseudodata present in the (often arbitrary) calibration densities. We recommend, as for any Bayesian analysis, that marginal prior and posterior distributions be carefully compared to determine whether signal is coming from the data or prior belief, especially for parameters of direct interest. This recommendation is not novel. However, given how rarely such checks are carried out in evolutionary biology, it bears repeating. Our results demonstrate the fundamental importance of prior/posterior comparisons in any Bayesian analysis, and we hope that they further encourage both researchers and journals to consistently adopt this crucial step as standard practice. Finally, we note that the results presented here do not refute the biological modeling concerns identified by Beaulieu et al. (2015). Both sets of issues remain apposite to the goals of accurate divergence time estimation, and only by considering them in tandem can we move forward more confidently.

Epigenetics and developmental plasticity in orthopteroid insects

Developmental plasticity is a key driver of the extraordinary ecological success of insects. Epigenetic mechanisms provide an important link between the external stimuli that initiate polyphenisms, and the stable changes in gene expression that govern alternative insect morphs. We review the epigenetics of orthopteroid insects, focussing on recent research on locusts and termites, two groups which display high levels of phenotypic plasticity, and for which genome sequences have become available in recent years. We examine research on the potential role of DNA methylation, histone modifications, and non-coding RNAs in the regulation of gene expression in these insects. DNA methylation patterns in orthopteroids share a number of characteristics with those of hymenopteran insects, although methylation levels are much higher, and extend to introns and repeat elements. Future examinations of epigenetic mechanisms in these insects will benefit from comparison of tissues from aged-matched individuals from alternative morphs, and adequate biological replication.

Evaluating the Impact of Genomic Data and Priors on Bayesian Estimates of the Angiosperm Evolutionary Timescale

The evolutionary timescale of angiosperms has long been a key question in biology. Molecular estimates of this timescale have shown considerable variation, being influenced by differences in taxon sampling, gene sampling, fossil calibrations, evolutionary models, and choices of priors. Here, we analyze a data set comprising 76 protein-coding genes from the chloroplast genomes of 195 taxa spanning 86 families, including novel genome sequences for 11 taxa, to evaluate the impact of models, priors, and gene sampling on Bayesian estimates of the angiosperm evolutionary timescale. Using a Bayesian relaxed molecular-clock method, with a core set of 35 minimum and two maximum fossil constraints, we estimated that crown angiosperms arose 221 (251-192) Ma during the Triassic. Based on a range of additional sensitivity and subsampling analyses, we found that our date estimates were generally robust to large changes in the parameters of the birth-death tree prior and of the model of rate variation across branches. We found an exception to this when we implemented fossil calibrations in the form of highly informative gamma priors rather than as uniform priors on node ages. Under all other calibration schemes, including trials of seven maximum age constraints, we consistently found that the earliest divergences of angiosperm clades substantially predate the oldest fossils that can be assigned unequivocally to their crown group. Overall, our results and experiments with genome-scale data suggest that reliable estimates of the angiosperm crown age will require increased taxon sampling, significant methodological changes, and new information from the fossil record. [Angiospermae, chloroplast, genome, molecular dating, Triassic.].

A Triassic-Jurassic window into the evolution of Lepidoptera

On the basis of an assemblage of fossilized wing scales recovered from latest Triassic and earliest Jurassic sediments from northern Germany, we provide the earliest evidence for Lepidoptera (moths and butterflies). The diverse scales confirm a (Late) Triassic radiation of lepidopteran lineages, including the divergence of the Glossata, the clade that comprises the vast multitude of extant moths and butterflies that have a sucking proboscis. The microfossils extend the minimum calibrated age of glossatan moths by ca. 70 million years, refuting ancestral association of the group with flowering plants. Development of the proboscis may be regarded as an adaptive innovation to sucking free liquids for maintaining the insect's water balance under arid conditions. Pollination drops secreted by a variety of Mesozoic gymnosperms may have been non-mutualistically exploited as a high-energy liquid source. The early evolution of the Lepidoptera was probably not severely interrupted by the end-Triassic biotic crisis.

Bugs Are Not to Be Silenced: Small RNA Pathways and Antiviral Responses in Insects

Like every other organism on Earth, insects are infected with viruses, and they rely on RNA interference (RNAi) mechanisms to circumvent viral infections. A remarkable characteristic of RNAi is that it is both broadly acting, because it is triggered by double-stranded RNA molecules derived from virtually any virus, and extremely specific, because it targets only the particular viral sequence that initiated the process. Reviews covering the different facets of the RNAi antiviral immune response in insects have been published elsewhere. In this review, we build a framework to guide future investigation. We focus on the remaining questions and avenues of research that need to be addressed to move the field forward, including issues such as the activity of viral suppressors of RNAi, comparative genomics, the development of detailed maps of the subcellular localization of viral replication complexes with the RNAi machinery, and the regulation of the antiviral RNAi response.

Reconstructing the phylogeny of Blattodea: robust support for interfamilial relationships and major clades

Cockroaches are among the most recognizable of all insects. In addition to their role as pests, they play a key ecological role as decomposers. Despite numerous studies of cockroach phylogeny in recent decades, relationships among most major lineages are yet to be resolved. Here we examine phylogenetic relationships among cockroaches based on five genes (mitochondrial 12S rRNA, 16S rRNA, COII; nuclear 28S rRNA and histone H3), and infer divergence times on the basis of 8 fossils. We included in our analyses sequences from 52 new species collected in China, representing 7 families. These were combined with data from a recent study that examined these same genes from 49 species, resulting in a significant increase in taxa analysed. Three major lineages, Corydioidea, Blaberoidea, and Blattoidea were recovered, the latter comprising Blattidae, Tryonicidae, Lamproblattidae, Anaplectidae, Cryptocercidae and Isoptera. The estimated age of the split between Mantodea and Blattodea ranged from 204.3 Ma to 289.1 Ma. Corydioidea was estimated to have diverged 209.7 Ma (180.5-244.3 Ma 95% confidence interval [CI]) from the remaining Blattodea. The clade Blattoidea diverged from their sister group, Blaberoidea, around 198.3 Ma (173.1-229.1 Ma). The addition of the extra taxa in this study has resulted in significantly higher levels of support for a number of previously recognized groupings.

Light Sheet-based Fluorescence Microscopy of Living or Fixed and Stained Tribolium castaneum Embryos

The red flour beetle Tribolium castaneum has become an important insect model organism in developmental genetics and evolutionary developmental biology. The observation of Tribolium embryos with light sheet-based fluorescence microscopy has multiple advantages over conventional widefield and confocal fluorescence microscopy. Due to the unique properties of a light sheet-based microscope, three dimensional images of living specimens can be recorded with high signal-to-noise ratios and significantly reduced photo-bleaching as well as photo-toxicity along multiple directions over periods that last several days. With more than four years of methodological development and a continuous increase of data, the time seems appropriate to establish standard operating procedures for the usage of light sheet technology in the Tribolium community as well as in the insect community at large. This protocol describes three mounting techniques suitable for different purposes, presents two novel custom-made transgenic Tribolium lines appropriate for long-term live imaging, suggests five fluorescent dyes to label intracellular structures of fixed embryos and provides information on data post-processing for the timely evaluation of the recorded data. Representative results concentrate on long-term live imaging, optical sectioning and the observation of the same embryo along multiple directions. The respective datasets are provided as a downloadable resource. Finally, the protocol discusses quality controls for live imaging assays, current limitations and the applicability of the outlined procedures to other insect species. This protocol is primarily intended for developmental biologists who seek imaging solutions that outperform standard laboratory equipment. It promotes the continuous attempt to close the gap between the technically orientated laboratories/communities, which develop and refine microscopy methodologically, and the life science laboratories/communities, which require 'plug-and-play' solutions to technical challenges. Furthermore, it supports an axiomatic approach that moves the biological questions into the center of attention.

The impacts of drift and selection on genomic evolution in insects

Genomes evolve through a combination of mutation, drift, and selection, all of which act heterogeneously across genes and lineages. This leads to differences in branch-length patterns among gene trees. Genes that yield trees with the same branch-length patterns can be grouped together into clusters. Here, we propose a novel phylogenetic approach to explain the factors that influence the number and distribution of these gene-tree clusters. We apply our method to a genomic dataset from insects, an ancient and diverse group of organisms. We find some evidence that when drift is the dominant evolutionary process, each cluster tends to contain a large number of fast-evolving genes. In contrast, strong negative selection leads to many distinct clusters, each of which contains only a few slow-evolving genes. Our work, although preliminary in nature, illustrates the use of phylogenetic methods to shed light on the factors driving rate variation in genomic evolution.

Erythropoietin-Mediated Neuroprotection in Insects Suggests a Prevertebrate Evolution of Erythropoietin-Like Signaling

The cytokine erythropoietin (Epo) mediates protective and regenerative functions in mammalian nervous systems via activation of poorly characterized receptors that differ from the "classical" homodimeric Epo receptor expressed on erythroid progenitor cells. Epo genes have been identified in vertebrate species ranging from human to fish, suggesting that Epo signaling evolved earlier than the vertebrate lineage. Studies on insects (Locusta migratoria, Chorthippus biguttulus, Tribolium castaneum) revealed Epo-mediated neuroprotection and neuroregeneration. Recombinant human Epo (rhEpo) prevents apoptosis by binding to a janus kinase-associated receptor, stimulation of STAT transcription factors, and generation of factors that prevent the activation of proapoptotic caspases. Insect neurons were also protected by a neuroprotective but nonerythropoietic Epo splice variant, suggesting similarity with mammalian neuroprotective but not with homodimeric "classical" Epo receptors. Additionally, rhEpo promotes the regeneration of neurites in primary cultured insect brain neurons and after nerve crush in an in vivo preparation. In contrast to neuroprotective and regenerative effects shared with mammalian species, no evidence for a role of Epo signaling in the regulation of neuro- or gliogenesis was found in insects. Similar structural and functional characteristics of the Epo binding receptors, partly shared transduction pathways that prevent apoptosis and the functional implication in neuroprotective and neuroregenerative processes in both mammalian and insect species, suggest that Epo-like signaling was already established in their last common ancestor. Originally functioning as a tissue-protective response to unfavorable physiological situations, cell injury, and pathogen invasion, Epo was later adapted as a humoral regulator of erythropoiesis in the vertebrate lineage.

Incorporating fossils into hypotheses of insect phylogeny

Fossils represent stem and crown lineages, and their inclusion in phylogenetic reconstruction influences branch lengths, topology, and divergence time estimation. In addition, paleontological data may inform trends in morphological evolution as well as biogeographic history. Here we review the incorporation of fossils in studies of insect evolution, from morphological analyses to combined 'total evidence' node dating analyses. We discuss challenges associated with fossil based phylogenetics, and suggest best practices for use in tree reconstruction.

Advances using molecular data in insect systematics

The size of molecular datasets has been growing exponentially since the mid 1980s, and new technologies have now dramatically increased the slope of this increase. New datasets include genomes, transcriptomes, and hybrid capture data, producing hundreds or thousands of loci. With these datasets, we are approaching a consensus on the higher level insect phylogeny. Huge datasets can produce new challenges in interpreting branch support, and new opportunities in developing better models and more sophisticated partitioning schemes. Dating analyses are improving as we recognize the importance of careful fossil calibration selection. With thousands of genes now available, coalescent methods have come of age. Barcode libraries continue to expand, and new methods are being developed for incorporating them into phylogenies with tens of thousands of individuals.

A molecular palaeobiological exploration of arthropod terrestrialization

Understanding animal terrestrialization, the process through which animals colonized the land, is crucial to clarify extant biodiversity and biological adaptation. Arthropoda (insects, spiders, centipedes and their allies) represent the largest majority of terrestrial biodiversity. Here we implemented a molecular palaeobiological approach, merging molecular and fossil evidence, to elucidate the deepest history of the terrestrial arthropods. We focused on the three independent, Palaeozoic arthropod terrestrialization events (those of Myriapoda, Hexapoda and Arachnida) and showed that a marine route to the colonization of land is the most likely scenario. Molecular clock analyses confirmed an origin for the three terrestrial lineages bracketed between the Cambrian and the Silurian. While molecular divergence times for Arachnida are consistent with the fossil record, Myriapoda are inferred to have colonized land earlier, substantially predating trace or body fossil evidence. An estimated origin of myriapods by the Early Cambrian precedes the appearance of embryophytes and perhaps even terrestrial fungi, raising the possibility that terrestrialization had independent origins in crown-group myriapod lineages, consistent with morphological arguments for convergence in tracheal systems.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Breakdown of coevolution between symbiotic bacteria Wolbachia and their filarial hosts

Wolbachia is an alpha-proteobacterial symbiont widely distributed in arthropods. Since the identification of Wolbachia in certain animal-parasitic nematodes (the Onchocercidae or filariae), the relationship between arthropod and nematode Wolbachia has attracted great interest. The obligate symbiosis in filariae, which renders infected species susceptible to antibiotic chemotherapy, was held to be distinct from the Wolbachia-arthropod relationship, typified by reproductive parasitism. While co-evolutionary signatures in Wolbachia-arthropod symbioses are generally weak, reflecting horizontal transmission events, strict co-evolution between filariae and Wolbachia has been reported previously. However, the absence of close outgroups for phylogenetic studies prevented the determination of which host group originally acquired Wolbachia. Here, we present the largest co-phylogenetic analysis of Wolbachia in filariae performed to date including: (i) a screening and an updated phylogeny of Wolbachia; (ii) a co-phylogenetic analysis; and (iii) a hypothesis on the acquisition of Wolbachia infection. First, our results show a general overestimation of Wolbachia occurrence and support the hypothesis of an ancestral absence of infection in the nematode phylum. The accuracy of supergroup J is also underlined. Second, although a global pattern of coevolution remains, the signal is derived predominantly from filarial clades associated with Wolbachia in supergroups C and J. In other filarial clades, harbouring Wolbachia supergroups D and F, horizontal acquisitions and secondary losses are common. Finally, our results suggest that supergroup C is the basal Wolbachia clade within the Ecdysozoa. This hypothesis on the origin of Wolbachia would change drastically our understanding of Wolbachia evolution.

Single origin of the Mascarene stick insects: ancient radiation on sunken islands?

BACKGROUND

The study of islands as model systems plays a key role in understanding many evolutionary processes. Knowledge of the historical events leading to present-day island communities is pivotal for exploring fundamental mechanisms of speciation and adaptation. The remote Mascarene archipelago (Mauritius, Réunion, Rodrigues), considered to be the product of an age-progressive trend of north-to-south volcanic activity in the Indian Ocean, hosts a remarkably diverse, endemic and threatened concentration of flora and fauna that has traditionally been considered to be biogeographically related to Madagascar and Africa. To explore the evolutionary diversity of the Mascarene stick insects (Phasmatodea), we constructed a global phylogeny from approximately 2.4 kb of mitochondrial and nuclear sequence data of more than 120 species representing all major phasmatodean lineages.

RESULTS

Based on the obtained time-calibrated molecular tree we demonstrate that the current phasmid community of the Mascarene archipelago, which consists of members of four presumably unrelated traditional subfamilies, is the result of a single ancient dispersal event from Australasia and started radiating between 16-29 million years ago, significantly predating the age of Mauritius (8-10 million years).

CONCLUSIONS

We propose that the Mascarene stick insects diversified on landmasses now eroded away, presumably to the north of Mauritius. In consequence, ancient islands have probably persisted in the Indian Ocean until the emergence of Mauritius and not only served as stepping stones for colonisation events during sea-level lowstands, but as long-lasting cradles of evolution. These ancient landmasses most likely allowed for adaptive speciation and served as significant sources of diversity that contributed to the biomes of the Mascarene archipelago and the megadiverse Madagascar.