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

Evolution across the adaptive landscape in a hyperdiverse beetle radiation

The extraordinary diversification of beetles on Earth is a textbook example of adaptive evolution. Yet, the tempo and drivers of this super-radiation remain largely unclear. Here, we address this problem by investigating macroevolutionary dynamics in darkling beetles (Coleoptera: Tenebrionidae), one of the most ecomorphologically diverse beetle families (with over 30,000 species). Using multiple genomic datasets and analytical approaches, we resolve the long-standing inconsistency over deep relationships in the family. In conjunction with a landmark-based dataset of body shape morphology, we show that the evolutionary history of darkling beetles is marked by ancient rapid radiations, frequent ecological transitions, and rapid bursts of morphological diversification. On a global scale, our analyses uncovered a significant pulse of phenotypic diversification proximal to the Cretaceous-Palaeogene (K/Pg) mass extinction and convergence of body shape associated with recurrent ecological specializations. On a regional scale, two major Australasian radiations, the Adeliini and the Heleine clade, exhibited contrasting patterns of ecomorphological diversification, representing phylogenetic niche conservatism versus adaptive radiation. Our findings align with the Simpsonian model of adaptive evolution across the macroevolutionary landscape and highlight a significant role of ecological opportunity in driving the immense ecomorphological diversity in a hyperdiverse beetle group.

Morphology and distribution of sensilla on head appendages in the water beetle Hygrobia hermanni (Coleoptera: Adephaga: Hygrobiidae)

Sensilla on head appendages were studied in detail for the first time in a member of the relict family Hygrobiidae (squeak beetles), closely related to Dytiscidae (diving beetles). Adult and third instar larval stage specimens of Hygrobia hermanni (Fabricius, 1775) were examined using scanning electron microscopy, focusing on antennae, palps and larval mandibles. In total, 37 sensilla subtypes are described, including 22 observed in the adult (basiconica: 3; Böhm's bristles: 2; circumvallate sensilla: 2; coeloconica: 10; ovoid placodea: 3; digitiform placodea: 2) and 16 in the larva (basiconica: 4; campaniformia: 1; chaetica: 4; coeloconica: 5; trichodea: 1; unnamed: 1). Only one subtype (of sensilla coeloconica) was shared between the adult and the larva. Autapomorphies of Hygrobiidae and Dytiscidae, and putatively shared derived characters (synapomorphies) of Hygrobiidae + Dytiscidae are discussed. Among the latter, the most remarkable is the acquisition of a special sensory field, located on the apical segment of the adult maxillary palp, subapically and postero-dorsally. This sensory field is made up of ovoid multiporous sensilla placodea otherwise present on the anterior (internal) surface of antennal segments, suggesting that in a common ancestor of Hygrobiidae and Dytiscidae, maxillary palps might have taken over enhanced capacities of longe-range molecule detection.

New Species of Paussus, Subgenus Scaphipaussus (Coleoptera: Carabidae: Paussinae), from Southeast Asia Reveal Ambiguities in Species Group Limits and High Species Diversity in the Oriental Region

Paussus, commonly known as ant nest beetles, is the most diverse genus of Paussinae (Coleoptera: Carabidae) with a very complex taxonomic history. Biodiversity research in Southeast and South Asia yields new species that can contribute to a better understanding of the morphological disparity and species-group or subgenus delimitation. Here, we describe nine new species from Southeast Asia and China: Paussus (Scaphipaussus) fencli sp. nov. (China), P. (S.) mawdsleyi sp. nov. (Borneo), P. (S.) bakeri sp. nov. (Philippines), P. (S.) jendeki sp. nov. (Laos), P. (S.) saueri sp. nov. (India), P. (S.) annamensis sp. nov. (Vietnam), P. (S.) phoupanensis sp. nov. (Laos, Vietnam), P. (S.) bilyi sp. nov. (Thailand), and P. (S.) haucki sp. nov. (Thailand). We also bring new data on P. (S.) corporaali Reichensperger, 1927 (Java) and P. (S.) madurensis Wasmann, 1913 (India). Besides formal descriptions, we provide photographs of the habitus in the dorsal and dorsolateral view, antennal club, head crest, and male genitalia if the male is available. Based on the comparison of new and earlier described species, we show that the antennae are highly diverse within the Scaphipaussus. Considering other characters, some species are placed in Scaphipaussus, but they differ from putative relatives in the antennal morphology. The presence of the frontal protuberances and crests is a more reliable character. Additional species show that Scaphipaussus is most diverse in southeastern Asia, especially in Indo-Burma. Concerning its supposed late Miocene origin, the group underwent rapid radiation. The species diversity of Scaphipaussus almost doubled in the last decade, and it is highly probable that further species will be described in the future.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01944-0.

Phylogenomic and mitogenomic data can accelerate inventorying of tropical beetles during the current biodiversity crisis

Conservation efforts must be evidence-based, so rapid and economically feasible methods should be used to quantify diversity and distribution patterns. We have attempted to overcome current impediments to the gathering of biodiversity data by using integrative phylogenomic and three mtDNA fragment analyses. As a model, we sequenced the Metriorrhynchini beetle fauna, sampled from ~700 localities in three continents. The species-rich dataset included ~6,500 terminals, ~1,850 putative species delimited at 5% uncorrected pairwise threshold, possibly ~1,000 of them unknown to science. Neither type of data could alone answer our questions on biodiversity and phylogeny. The phylogenomic backbone enabled the integrative delimitation of robustly defined natural genus-group units that will inform future research. Using constrained mtDNA analysis, we identified the spatial structure of species diversity, very high species-level endemism, and a biodiversity hotspot in New Guinea. We suggest that focused field research and subsequent laboratory and bioinformatic workflow steps would substantially accelerate the inventorying of any hyperdiverse tropical group with several thousand species. The outcome would be a scaffold for the incorporation of further data from environmental sequencing and ecological studies. The database of sequences could set a benchmark for the spatiotemporal evaluation of biodiversity, would support evidence-based conservation planning, and would provide a robust framework for systematic, biogeographic, and evolutionary studies.

Anchored Phylogenomics, Evolution and Systematics of Elateridae: Are All Bioluminescent Elateroidea Derived Click Beetles?

Simple Summary

In the era of phylogenomics, new molecular sequencing and computational techniques can aid in resolving phylogenetic relationships that were previously intractable by morphological or limited molecular data. In this study, we used anchored hybrid enrichment—designed to recover DNA sequences from hundreds of single-copy orthologous genes—to resolve the phylogeny of the Elateridae (click-beetles) and establish their placement within superfamily Elateroidea. The resulting data were compatible with published transcriptomes, allowing for integrating our dataset with previously published data. Using a wide range of analyses on these molecular data, we tested hypotheses long-debated in the morphological literature and also the robustness of our phylogenetic inferences. Our results placed the bioluminescent lampyroids (fireflies and relatives) within the click-beetles, challenging the current classification of Elateridae, Lampyridae, Phengodidae, and Rhagophthalmidae. However, despite the large amount of molecular data analyzed, a few nodes with conflicting phylogenetic signals could not be unambiguously resolved. Overall, we recovered well-resolved tree topologies that will serve as a framework for further systematic and evolutionary studies of click-beetles. This work further demonstrates that the click-beetle lineage contains not only pest wireworms, but also many species that benefit agriculture.

Abstract

Click-beetles (Coleoptera: Elateridae) are an abundant, diverse, and economically important beetle family that includes bioluminescent species. To date, molecular phylogenies have sampled relatively few taxa and genes, incompletely resolving subfamily level relationships. We present a novel probe set for anchored hybrid enrichment of 2260 single-copy orthologous genes in Elateroidea. Using these probes, we undertook the largest phylogenomic study of Elateroidea to date (99 Elateroidea, including 86 Elateridae, plus 5 non-elateroid outgroups). We sequenced specimens from 88 taxa to test the monophyly of families, subfamilies and tribes. Maximum likelihood and coalescent phylogenetic analyses produced well-resolved topologies. Notably, the included non-elaterid bioluminescent families (Lampyridae + Phengodidae + Rhagophthalmidae) form a clade within the otherwise monophyletic Elateridae, and Sinopyrophoridae may not warrant recognition as a family. All analyses recovered the elaterid subfamilies Elaterinae, Agrypninae, Cardiophorinae, Negastriinae, Pityobiinae, and Tetralobinae as monophyletic. Our results were conflicting on whether the hypnoidines are sister to Dendrometrinae or Cardiophorinae + Negastriinae. Moreover, we show that fossils with the eucnemid-type frons and elongate cylindrical shape may belong to Eucnemidae, Elateridae: Thylacosterninae, ancestral hard-bodied cantharoids or related extinct groups. Proposed taxonomic changes include recognition of Plastocerini as a tribe in Dendrometrinae and Hypnoidinae stat. nov. as a subfamily within Elateridae.

Root Digger: a root placement program for phylogenetic trees

BACKGROUND

In phylogenetic analysis, it is common to infer unrooted trees. However, knowing the root location is desirable for downstream analyses and interpretation. There exist several methods to recover a root, such as molecular clock analysis (including midpoint rooting) or rooting the tree using an outgroup. Non-reversible Markov models can also be used to compute the likelihood of a potential root position.

RESULTS

We present a software called RootDigger which uses a non-reversible Markov model to compute the most likely root location on a given tree and to infer a confidence value for each possible root placement. We find that RootDigger is successful at finding roots when compared to similar tools such as IQ-TREE and MAD, and will occasionally outperform them. Additionally, we find that the exhaustive mode of RootDigger is useful in quantifying and explaining uncertainty in rooting positions.

CONCLUSIONS

RootDigger can be used on an existing phylogeny to find a root, or to asses the uncertainty of the root placement. RootDigger is available under the MIT licence at https://www.github.com/computations/root_digger .

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.

Data Types and the Phylogeny of Neoaves

The phylogeny of Neoaves, the largest clade of extant birds, has remained unclear despite intense study. The difficulty associated with resolving the early branches in Neoaves is likely driven by the rapid radiation of this group. However, conflicts among studies may be exacerbated by the data type analyzed. For example, analyses of coding exons typically yield trees that place Strisores (nightjars and allies) sister to the remaining Neoaves, while analyses of non-coding data typically yield trees where Mirandornites (flamingos and grebes) is the sister of the remaining Neoaves. Our understanding of data type effects is hampered by the fact that previous analyses have used different taxa, loci, and types of non-coding data. Herein, we provide strong corroboration of the data type effects hypothesis for Neoaves by comparing trees based on coding and non-coding data derived from the same taxa and gene regions. A simple analytical method known to minimize biases due to base composition (coding nucleotides as purines and pyrimidines) resulted in coding exon data with increased congruence to the non-coding topology using concatenated analyses. These results improve our understanding of the resolution of neoavian phylogeny and point to a challenge—data type effects—that is likely to be an important factor in phylogenetic analyses of birds (and many other taxonomic groups). Using our results, we provide a summary phylogeny that identifies well-corroborated relationships and highlights specific nodes where future efforts should focus.

Phylogenomics of Parasitic and Nonparasitic Lice (Insecta: Psocodea): Combining Sequence Data and Exploring Compositional Bias Solutions in Next Generation Data Sets

The insect order Psocodea is a diverse lineage comprising both parasitic (Phthiraptera) and nonparasitic members (Psocoptera). The extreme age and ecological diversity of the group may be associated with major genomic changes, such as base compositional biases expected to affect phylogenetic inference. Divergent morphology between parasitic and nonparasitic members has also obscured the origins of parasitism within the order. We conducted a phylogenomic analysis on the order Psocodea utilizing both transcriptome and genome sequencing to obtain a data set of 2370 orthologous genes. All phylogenomic analyses, including both concatenated and coalescent methods suggest a single origin of parasitism within the order Psocodea, resolving conflicting results from previous studies. This phylogeny allows us to propose a stable ordinal level classification scheme that retains significant taxonomic names present in historical scientific literature and reflects the evolution of the group as a whole. A dating analysis, with internal nodes calibrated by fossil evidence, suggests an origin of parasitism that predates the K-Pg boundary. Nucleotide compositional biases are detected in third and first codon positions and result in the anomalous placement of the Amphientometae as sister to Psocomorpha when all nucleotide sites are analyzed. Likelihood-mapping and quartet sampling methods demonstrate that base compositional biases can also have an effect on quartet-based methods.[Illumina; Phthiraptera; Psocoptera; quartet sampling; recoding methods.].

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.

Data curation and modeling of compositional heterogeneity in insect phylogenomics: A case study of the phylogeny of Dytiscoidea (Coleoptera: Adephaga)

Diving beetles and their allies are an almost ubiquitous group of freshwater predators. Knowledge of the phylogeny of the adephagan superfamily Dytiscoidea has significantly improved since the advent of molecular phylogenetics. However, despite recent comprehensive phylogenomic studies, some phylogenetic relationships among the constituent families remain elusive. In particular, the position of the family Hygrobiidae remains uncertain. We address these issues by re-analyzing recently published phylogenomic datasets for Dytiscoidea, using approaches to reduce compositional heterogeneity and adopting a site-heterogeneous mixture model. We obtained a consistent, well-resolved, and strongly supported tree. Consistent with previous studies, our analyses support Aspidytidae as the monophyletic sister group of Amphizoidae, and more importantly, Hygrobiidae as the sister of the diverse Dytiscidae, in agreement with morphology-based phylogenies. Our analyses provide a backbone phylogeny of Dytiscoidea, which lays the foundation for better understanding the evolution of morphological characters, life habits, and feeding behaviors of dytiscoid beetles.

Genomic signatures accompanying the dietary shift to phytophagy in polyphagan beetles

BACKGROUND

The diversity and evolutionary success of beetles (Coleoptera) are proposed to be related to the diversity of plants on which they feed. Indeed, the largest beetle suborder, Polyphaga, mostly includes plant eaters among its approximately 315,000 species. In particular, plants defend themselves with a diversity of specialized toxic chemicals. These may impose selective pressures that drive genomic diversification and speciation in phytophagous beetles. However, evidence of changes in beetle gene repertoires driven by such interactions remains largely anecdotal and without explicit hypothesis testing.

RESULTS

We explore the genomic consequences of beetle-plant trophic interactions by performing comparative gene family analyses across 18 species representative of the two most species-rich beetle suborders. We contrast the gene contents of species from the mostly plant-eating suborder Polyphaga with those of the mainly predatory Adephaga. We find gene repertoire evolution to be more dynamic, with significantly more adaptive lineage-specific expansions, in the more speciose Polyphaga. Testing the specific hypothesis of adaptation to plant feeding, we identify families of enzymes putatively involved in beetle-plant interactions that underwent adaptive expansions in Polyphaga. There is notable support for the selection hypothesis on large gene families for glutathione S-transferase and carboxylesterase detoxification enzymes.

CONCLUSIONS

Our explicit modeling of the evolution of gene repertoires across 18 species identifies putative adaptive lineage-specific gene family expansions that accompany the dietary shift towards plants in beetles. These genomic signatures support the popular hypothesis of a key role for interactions with plant chemical defenses, and for plant feeding in general, in driving beetle diversification.