Potential Contribution of Ancient Introgression to the Evolution of a Derived Reproductive Strategy in Ricefishes.

Transitions from no parental care to extensive care are costly and involve major changes in life history, behavior, and morphology. Nevertheless, in Sulawesi ricefishes, pelvic brooding evolved from transfer brooding in two distantly related lineages within the genera Adrianichthys and Oryzias, respectively. Females of pelvic brooding species carry their eggs attached to their belly until the fry hatches. Despite their phylogenetic distance, both pelvic brooding lineages share a set of external morphological traits. A recent study found no direct gene flow between pelvic brooding lineages, suggesting independent evolution of the derived reproductive strategy. Convergent evolution can, however, also rely on repeated sorting of preexisting variation of an admixed ancestral population, especially when subjected to similar external selection pressures. We thus used a multispecies coalescent model and D-statistics to identify gene-tree-species-tree incongruencies, to evaluate the evolution of pelvic brooding with respect to interspecific gene flow not only between pelvic brooding lineages but also between pelvic brooding lineages and other Sulawesi ricefish lineages. We found a general network-like evolution in Sulawesi ricefishes, and as previously reported, we detected no gene flow between the pelvic brooding lineages. Instead, we found hybridization between the ancestor of pelvic brooding Oryzias and the common ancestor of the Oryzias species from the Lake Poso area. We further detected signs of introgression within the confidence interval of a quantitative trait locus associated with pelvic brooding in O. eversi. Our results hint toward a contribution of ancient standing genetic variation to the evolution of pelvic brooding in Oryzias.

Pitfalls of the site-concordance factor (sCF) as measure of phylogenetic branch support.

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.

Phylogenetic Revision and Patterns of Host Specificity in the Fungal Subphylum Entomophthoromycotina.

The Entomophthoromycotina, a subphylum close to the root of terrestrial fungi with a bias toward insects as their primary hosts, has been notoriously difficult to categorize taxonomically for decades. Here, we reassess the phylogeny of this group based on conserved genes encoding ribosomal RNA and RNA polymerase II subunits, confirming their general monophyly, but challenging previously assumed taxonomic relationships within and between particular clades. Furthermore, for the prominent, partially human-pathogenic taxon Conidiobolus, a new type species C. coronatus is proposed in order to compensate for the unclear, presumably lost previous type species C. utriculosus Brefeld 1884. We also performed an exhaustive survey of the broad host spectrum of the Entomophthoromycotina, which is not restricted to insects alone, and investigated potential patterns of co-evolution across their megadiverse host range. Our results suggest multiple independent origins of parasitism within this subphylum and no apparent co-evolutionary events with any particular host lineage. However, Pterygota (i.e., winged insects) clearly constitute the most dominantly parasitized superordinate host group. This appears to be in accordance with an increased dispersal capacity mediated by the radiation of the Pterygota during insect evolution, which has likely greatly facilitated the spread, infection opportunities, and evolutionary divergence of the Entomophthoromycotina as well.

Transcriptomics provides a robust framework for the relationships of the major clades of cladobranch sea slugs (Mollusca, Gastropoda, Heterobranchia), but fails to resolve the position of the enigmatic genus Embletonia.

BACKGROUND: The soft-bodied cladobranch sea slugs represent roughly half of the biodiversity of marine nudibranch molluscs 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 targets for drug discovery), the evolutionary history of cladobranch sea slugs is not yet fully understood. RESULTS: To enlarge the current knowledge on the phylogenetic relationships, we generated new transcriptome data for 19 species of cladobranch sea slugs and two additional outgroup taxa (Berthella plumula and Polycera quadrilineata). We complemented our taxon sampling with previously published transcriptome data, resulting in a final data set covering 56 species from all but one accepted cladobranch superfamilies. We assembled all transcriptomes using six different assemblers, selecting those assemblies that provided the largest amount of potentially phylogenetically informative sites. Quality-driven compilation of data sets resulted in four different supermatrices: two with full coverage of genes per species (446 and 335 single-copy protein-coding genes, respectively) and two with a less stringent coverage (667 genes with 98.9% partition coverage and 1767 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 maximal statistical support for all major splits and phylogenetic relationships at the family level. Besides the questionable position of Noumeaella rubrofasciata, rendering the Facelinidae as polyphyletic, the only notable discordance between the inferred trees is the position of Embletonia 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. CONCLUSIONS: Our data matrices and the inferred trees can serve as a solid foundation for future work on the taxonomy and evolutionary history of Cladobranchia. The placement of E. pulchra, however, proves challenging, even with large data sets and various optimization strategies. Moreover, quartet mapping results show that confounding signal present in the data is sufficient to explain the inferred position of E. pulchra, again leaving its phylogenetic position as an enigma.

Monophyletic blowflies revealed by phylogenomics.

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.

Development and evaluation of a custom bait design based on 469 single-copy protein-coding genes for exon capture of isopods (Philosciidae: Haloniscus).

Transcriptome-based exon capture approaches, along with next-generation sequencing, are allowing for the rapid and cost-effective production of extensive and informative phylogenomic datasets from non-model organisms for phylogenetics and population genetics research. These approaches generally employ a reference genome to infer the intron-exon structure of targeted loci and preferentially select longer exons. However, in the absence of an existing and well-annotated genome, we applied this exon capture method directly, without initially identifying intron-exon boundaries for bait design, to a group of highly diverse Haloniscus (Philosciidae), paraplatyarthrid and armadillid isopods, and examined the performance of our methods and bait design for phylogenetic inference. Here, we identified an isopod-specific set of single-copy protein-coding loci, and a custom bait design to capture targeted regions from 469 genes, and analysed the resulting sequence data with a mapping approach and newly-created post-processing scripts. We effectively recovered a large and informative dataset comprising both short (<100 bp) and longer (>300 bp) exons, with high uniformity in sequencing depth. We were also able to successfully capture exon data from up to 16-year-old museum specimens along with more distantly related outgroup taxa, and efficiently pool multiple samples prior to capture. Our well-resolved phylogenies highlight the overall utility of this methodological approach and custom bait design, which offer enormous potential for application to future isopod, as well as broader crustacean, molecular studies.

Disentangling the aging gene expression network of termite queens.

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 many social insects (termites, ants and some bees) can live from a few years to decades. This makes social insects promising 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 a period of 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 ten 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, defense 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 of a stress-response mediated by the transcription factor cnc, a homolog of human nrf genes. CONCLUSIONS: 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 proteostasis network and stress as part of the aging process of social insect queens, shortly before queens die.

Transcriptomic Signatures of Ageing Vary in Solitary and Social Forms of an Orchid Bee.

Eusocial insect queens are remarkable in their ability to maximize both fecundity and longevity, thus escaping the typical trade-off between these two traits. Several mechanisms have been proposed to underlie the remolding of the trade-off, such as reshaping of the juvenile hormone (JH) pathway, or caste-specific susceptibility to oxidative stress. However, it remains a challenge to disentangle the molecular mechanisms underlying the remolding of the trade-off in eusocial insects from caste-specific physiological attributes that have subsequently arisen. The socially polymorphic orchid bee Euglossa viridissima represents an excellent model to address the role of sociality per se in longevity as it allows direct comparisons of solitary and social individuals within a common genetic background. We investigated gene expression and JH levels in young and old bees from both solitary and social nests. We found 902 genes to be differentially expressed with age in solitary females, including genes involved in oxidative stress, versus only 100 genes in social dominant females, and 13 genes in subordinate females. A weighted gene coexpression network analysis further highlights pathways related to ageing in this species, including the target of rapamycin pathway. Eleven genes involved in translation, apoptosis, and DNA repair show concurrent age-related expression changes in solitary but not in social females, representing potential differences based on social status. JH titers did not vary with age or social status. Our 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.

Comparative transcriptomic analysis of the mechanisms underpinning ageing and fecundity in social insects.

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 (TOR/IIS-JH-Lifespan and Fecundity) 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?'

Oxidative stress and senescence in social insects: a significant but inconsistent link?

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?'

Beyond Drosophila: resolving the rapid radiation of schizophoran flies with phylotranscriptomics.

BACKGROUND: The most species-rich radiation of animal life in the 66 million years following the Cretaceous extinction event is that of schizophoran flies: a third of fly diversity including Drosophila fruit fly model organisms, house flies, forensic blow flies, agricultural pest flies, and many other well and poorly known true flies. Rapid diversification has hindered previous attempts to elucidate the phylogenetic relationships among major schizophoran clades. A robust phylogenetic hypothesis for the major lineages containing these 55,000 described species would be critical to understand the processes that contributed to the diversity of these flies. We use protein encoding sequence data from transcriptomes, including 3145 genes from 70 species, representing all superfamilies, to improve the resolution of this previously intractable phylogenetic challenge. RESULTS: Our results support a paraphyletic acalyptrate grade including a monophyletic Calyptratae and the monophyly of half of the acalyptrate superfamilies. The primary branching framework of Schizophora is well supported for the first time, revealing the primarily parasitic Pipunculidae and Sciomyzoidea stat. rev. as successive sister groups to the remaining Schizophora. Ephydroidea, Drosophila's superfamily, is the sister group of Calyptratae. Sphaeroceroidea has modest support as the sister to all non-sciomyzoid Schizophora. We define two novel lineages corroborated by morphological traits, the 'Modified Oviscapt Clade' containing Tephritoidea, Nerioidea, and other families, and the 'Cleft Pedicel Clade' containing Calyptratae, Ephydroidea, and other families. Support values remain low among a challenging subset of lineages, including Diopsidae. The placement of these families remained uncertain in both concatenated maximum likelihood and multispecies coalescent approaches. Rogue taxon removal was effective in increasing support values compared with strategies that maximise gene coverage or minimise missing data. CONCLUSIONS: Dividing most acalyptrate fly groups into four major lineages is supported consistently across analyses. Understanding the fundamental branching patterns of schizophoran flies provides a foundation for future comparative research on the genetics, ecology, and biocontrol.

Four myriapod relatives - but who are sisters? No end to debates on relationships among the four major myriapod subgroups.

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.

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

An amendment to this paper has been published and can be accessed via the original article.

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.

A minimum reporting standard for multiple sequence alignments.

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.

Phylogenomics reveals the evolutionary timing and pattern of butterflies and moths.

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.

Identifying genetic markers for a range of phylogenetic utility-From species to family level.

Resolving the phylogenetic relationships of closely related species using a small set of loci is challenging as sufficient information may not be captured from a limited sample of the genome. Relying on few loci can also be problematic when conflict between gene-trees arises from incomplete lineage sorting and/or ongoing hybridization, problems especially likely in recently diverged lineages. Here, we developed a method using limited genomic resources that allows identification of many low copy candidate loci from across the nuclear and chloroplast genomes, design probes for target capture and sequence the captured loci. To validate our method we present data from Eucalyptus and Melaleuca, two large and phylogenetically problematic genera within the Myrtaceae family. With one annotated genome, one transcriptome and two whole-genome shotgun sequences of one Eucalyptus and four Melaleuca species, respectively, we identified 212 loci representing 263 kbp for targeted sequence capture and sequencing. Of these, 209 were successfully tested from 47 samples across five related genera of Myrtaceae. The average percentage of reads mapped back to the reference was 57.6% with coverage of more than 20 reads per position across 83.5% of the data. The methods developed here should be applicable across a large range of taxa across all kingdoms. The core methods are very flexible, providing a platform for various genomic resource availabilities and are useful from shallow to deep phylogenies.

An integrative phylogenomic approach illuminates the evolutionary history of cockroaches and termites (Blattodea).

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.

No signal of deleterious mutation accumulation in conserved gene sequences of extant asexual hexapods.

Loss of sex and recombination is generally assumed to impede the effectiveness of purifying selection and to result in the accumulation of slightly deleterious mutations. Empirical evidence for this has come from several studies investigating mutational load in a small number of individual genes. However, recent whole transcriptome based studies have yielded inconsistent results, hence questioning the validity of the assumption of mutational meltdown in asexual populations. Here, we study the effectiveness of purifying selection in eight asexual hexapod lineages and their sexual relatives, as present in the 1 K Insect Transcriptome Evolution (1KITE) project, covering eight hexapod groups. We analyse the accumulation of slightly deleterious nonsynonymous and synonymous point mutations in 99 single copy orthologue protein-coding loci shared among the investigated taxa. While accumulation rates of nonsynonymous mutations differed between genes and hexapod groups, we found no effect of reproductive mode on the effectiveness of purifying selection acting at nonsynonymous and synonymous sites. Although the setup of this study does not fully rule out nondetection of subtle effects, our data does not support the established consensus of asexual lineages undergoing 'mutational meltdown'.

Phylogenomics of the superfamily Dytiscoidea (Coleoptera: Adephaga) with an evaluation of phylogenetic conflict and systematic error.

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.