The Sirex woodwasp (Sirex noctilio; Hymenoptera: Siricidae) is among the most destructive invasive pests affecting Pinus plantations worldwide. Chemosensory systems offer promising targets for pest control strategies. The identification and characterization of chemosensation genes in non-model, economically significant insects such as S. noctilio, is an important first step towards the development of such control methods. Here we sequenced and assembled a draft genome of S. noctilio and performed RNA-sequencing of 15 olfactory and non-olfactory tissues to study the expression patterns of chemosensation-related genes. Specific genes, such as SnocOR16 and SnocSNMP1, displayed tissue- and sex-specific expression patterns, making them particularly intriguing for their potential roles in chemosensation and oviposition. As woodwasps and their related lineages form a sister group to the majority of other Hymenoptera, including Apocrita, insights into their gene repertoires are crucial for tracing the evolutionary history of chemosensory multigene families of this ecologically and economically significant insect order. This study enhances our understanding of the molecular mechanisms underlying S. noctilio chemosensation, paving the way for further research in chemical ecology and the functional characterization of S. noctilio chemosensation genes.

Despite sequencing advances, genomic resources remain uneven and rarely comparable by design. Turning today’s genomes into analysis-ready resources requires species-dense sampling anchored in natural history collections. Large-scale genome initiatives are transforming biodiversity genomics, but genomes remain taxonomically sparse and difficult to compare across species. Kapli and colleagues argue that natural history collections should play a central role in building species-dense genomic resources, creating the foundation needed for comparative and predictive biodiversity genomics.

Bibionomorpha is a large and diverse dipteran infraorder, of which its composition and family-level relationships have long been debated. In this study, we constructed a phylogenomic tree of Bibionomorpha using an extensive dataset of transcriptomic, genomic and anchored hybrid enrichment data. To further investigate the evolutionary timeline of the group, we also generated a fossil-calibrated time tree using data-driven calibration priors. Bibionomorpha, comprising Anisopodoidea, Bibionoidea, Scatopsoidea and Sciaroidea, is recovered as the sister group to Brachycera. Axymyiidae and Perissommatidae, which have been considered part of Bibionomorpha by some authorities, are instead recovered as sister to the Bibionomorpha + Brachycera. We discuss recalcitrant nodes within the infraorder, particularly regarding the placements of Bolitophilidae and Cecidomyiidae. Although tentative, phylogenetic network analysis suggests a possible reticulation event leading to Bolitophilidae, while mitochondrial data support a sex-specific hybridization event between ancestral Cecidomyiidae and Sciaridae + Diadocidiidae group. Time tree analyses suggest a Lower Triassic or deeper origin of Bibionomorpha, consistent with fossil evidence.

BACKGROUND:Dietary shifts are major evolutionary steps that shape ecological niches and biodiversity. The beetle family Coccinellidae, commonly known as ladybirds, first transitioned from a fungivorous to an insectivorous and subsequently a plant diet. However, the molecular basis of this dietary diversification remained unexplored. RESULTS:We investigated the molecular evolution of dietary shifts in ladybirds, focusing on the transitions from fungivory to carnivory (Coccinellidae) and from carnivory to herbivory (Epilachnini), by comparing 25 genomes and 62 transcriptomes of beetles. Our analysis shows that chemosensory gene families have undergone significant expansions at both nodes of diet change and were differentially expressed in feeding experiments, suggesting that they may be related to foraging. We found expansions of digestive and detoxifying gene families and losses of chitin-related digestive genes in the herbivorous ladybirds, and absence of most plant cell wall-degrading enzymes in the ladybirds dating from the transition to carnivory, likely indicating the effect of different digestion requirements on the gene repertoire. Immunity effector genes tend to emerge or have specific amino acid sequence compositions in carnivorous ladybirds and are downregulated under suboptimal dietary treatments, suggesting a potential function of these genes related to microbial symbionts in the sternorrhynchan prey. CONCLUSIONS:Our study provides a comprehensive comparative genomic analysis to address evolution of chemosensory, digestive, detoxifying, and immune genes associated with dietary shifts in ladybirds. Ladybirds can be considered a ubiquitous example of dietary shifts in insects, and thus a promising model system for evolutionary and applied biology.

Biodiversity knowledge, from genes to ecosystems, is crucial for addressing the biodiversity crisis. However, even in well-explored countries like Germany, much biodiversity remains unknown. Therefore, several research institutions are joining forces to conduct a comprehensive biodiversity inventory, combining broad taxonomic expertise with advanced technologies. By consolidating data across many organismic groups, the Unknown Germany initiative will significantly enhance conservation strategies and may serve as a model for similar efforts worldwide.

Illegal wildlife trade is a growing problem internationally. Poaching of animals not only leads to the extinction of populations and species but also has serious consequences for ecosystems and economies. This study introduces a molecular marker system that authorities can use to detect and substantiate wildlife trafficking. SNPSTR markers combine short tandem repeats with single nucleotide polymorphisms within an amplicon to increase discriminatory power. Within the FOGS (Forensic Genetics for Species Protection) project, we have established SNPSTR marker sets for 74 vertebrate species. On average, each set consists of 19 SNPSTR markers with 82 SNPs per set. More than 1300 SNPSTR markers and over 300 STR markers were identified. Also, through its biobanking pipeline, the FOGS project enabled the cryopreservation of somatic cells from 91 vertebrate species as well as viable tissues for later cell initiation from a further 109 species, providing future strategies for ex situ conservation. In addition, many more fixed tissues and DNA samples of endangered species were biobanked. Therefore, FOGS was an interdisciplinary study, combining molecular wildlife forensics and conservation tools. The SNPSTR sets and cell culture information are accessible through the FOGS database (https://fogs-portal.de/data) that is open to scientists, researchers, breeders and authorities worldwide to protect wildlife from illegal trade.

Caddisflies (Trichoptera) are among the most diverse groups of freshwater animals with more than 16 000 described species. They play a fundamental role in freshwater ecology and environmental engineering in streams, rivers and lakes. Because of this, they are frequently used as indicator organisms in biomonitoring programmes. Despite their importance, key questions concerning the evolutionary history of caddisflies, such as the timing and origin of larval case making, remain unanswered owing to the lack of a well-resolved phylogeny. Here, we estimated a phylogenetic tree using a combination of transcriptomes and targeted enrichment data for 207 species, representing 48 of 52 extant families and 174 genera. We calibrated and dated the tree with 33 carefully selected fossils. The first caddisflies originated approximately 295 million years ago in the Permian, and major suborders began to diversify in the Triassic. Furthermore, we show that portable case making evolved in three separate lineages, and shifts in diversification occurred in concert with key evolutionary innovations beyond case making.

Trade of wild-caught animals is illegal for many taxa and in many countries. Common regulatory procedures involve documentation and marking techniques. However, these procedures are subject to fraud and thus should be complemented by routine genetic testing in order to authenticate the captive-bred origin of animals intended for trade. A suitable class of genetic markers are SNPSTRs that combine a short tandem repeat (STR) and single nucleotide polymorphisms (SNPs) within one amplicon. This combined marker type can be used for genetic identification and for parentage analyses and in addition, provides insight into haplotype history. As a proof of principle, this study establishes a set of 20 SNPSTR markers for Athene noctua, one of the most trafficked owls in CITES Appendix II. These markers can be coamplified in a single multiplex reaction. Based on population data, the percentage of observed and expected heterozygosities of the markers ranged from 0.400 to 1.000 and 0.545 to 0.850, respectively. A combined probability of identity of 5.3 * 10−23 was achieved with the whole set, and combined parentage exclusion probabilities reached over 99.99%, even if the genotype of one parent was missing. A direct comparison of an owl family and an unrelated owl demonstrated the applicability of the SNPSTR set in parentage testing. The established SNPSTR set thus proved to be highly useful for identifying individuals and analysing parentage to determine wild or captive origin. We propose to implement SNPSTR-based routine certification in wildlife trade as a way to reveal animal laundering and misdeclaration of wild-caught animals.

Abstract Classifications of forest vegetation types and characterization of related species assemblages are important analytical tools for mapping and diversity monitoring of forest communities. The discrimination of forest communities is often based on β‐diversity, which can be quantified via numerous indices to derive compositional dissimilarity between samples. This study aims to evaluate the applicability of unsupervised classification for National Forest Inventory data from Georgia by comparing two cluster hierarchies. We calculated the mean basal area per hectare for each woody species across 1059 plot observations and quantified interspecies distances for all 87 species. Following an unspuervised cluster analysis, we compared the results derived from the species‐neutral dissimilarity (Bray‐Curtis) with those based on the Discriminating Avalanche dissimilarity, which incorporates interspecies phylogenetic variation. Incorporating genetic variation in the dissimilarity quantification resulted in a more nuanced discrimination of woody species assemblages and increased cluster coherence. Favorable statistics include the total number of clusters (23 vs. 20), mean distance within clusters (0.773 vs. 0.343), and within sum of squares (344.13 vs. 112.92). Clusters derived from dissimilarities that account for genetic variation showed a more robust alignment with biogeographical units, such as elevation and known habitats. We demonstrate that the applicability of unsupervised classification of species assemblages to large‐scale forest inventory data strongly depends on the underlying quantification of dissimilarity. Our results indicate that by incorporating phylogenetic variation, a more precise classification aligned with biogeographic units is attained. This supports the concept that the genetic signal of species assemblages reflects biogeographical patterns and facilitates more precise analyses for mapping, monitoring, and management of forest diversity.

Abstract All ecosystems face ecological challenges in this century. Therefore, it is becoming increasingly important to understand the ecology and degree of local adaptation of functionally important Arctic‐alpine biomes by looking at the most diverse taxon of metazoans: the Arthropoda. This is the first study to utilize metabarcoding in the Alpine tundra, providing insights into the effects of micro‐environmental parameters on alpha‐ and beta‐diversity of arthropods in such unique environments. To characterize arthropod diversity, pitfall traps were set at three middle‐alpine sampling sites in the Scandinavian mountain range in Norway during the snow‐free season in 2015. A metabarcoding approach was then used to determine the small‐scale biodiversity patterns of arthropods in the Alpine tundra. All DNA was extracted directly from the preservative EtOH from 27 pitfall traps. In order to identify the controlling environmental conditions, all sampling locations were equipped with automatic data loggers for permanent measurement of the microenvironmental conditions. The variables measured were: air temperature [°C] at 15 cm height, soil temperature [°C] at 15 cm depth, and soil moisture [vol.%] at 15 cm depth. A total of 233 Arthropoda OTUs were identified. The number of unique OTUs found per sampling location (ridge, south‐facing slope, and depression) was generally higher than the OTUs shared between the sampling locations, demonstrating that niche features greatly impact arthropod community structure. Our findings emphasize the fine‐scale heterogeneity of arctic–alpine ecosystems and provide evidence for trait‐based and niche‐driven adaptation. The spatial and temporal differences in arthropod diversity were best explained by soil moisture and soil temperature at the respective locations. Furthermore, our results show that arthropod diversity is underestimated in alpine‐tundra ecosystems using classical approaches and highlight the importance of integrating long‐term functional environmental data and modern taxonomic techniques into biodiversity research to expand our ecological understanding of fine‐ and meso‐scale biogeographical patterns.

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Scarab beetles (Scarabaeidae) are a diverse and ecologically important group of angiosperm-associated insects. As conventionally understood, scarab beetles comprise two major lineages: dung beetles and the phytophagous Pleurosticti. However, previous phylogenetic analyses have not been able to convincingly answer the question whether or not the two lineages form a monophyletic group. Here we report our results from phylogenetic analyses of more than 4,000 genes mined from transcriptomes of more than 50 species of Scarabaeidae and other Scarabaeoidea. Our results provide convincing support for the monophyly of Scarabaeidae, confirming the debated sister group relationship of dung beetles and phytophagous pleurostict scarabs. Supermatrix-based maximum likelihood and multispecies coalescent phylogenetic analyses strongly imply the subfamily Melolonthinae as currently understood being paraphyletic. We consequently suggest various changes in the systematics of Melolonthinae: Sericinae Kirby, 1837 stat. rest. and sensu n. to include the tribes Sericini, Ablaberini and Diphucephalini, and Sericoidinae Erichson, 1847 stat. rest. and sensu n. to include the tribes Automoliini, Heteronychini, Liparetrini, Maechidiini, Scitalini, Sericoidini, and Phyllotocini. Both subfamilies appear to consistently form a monophyletic sister group to all remaining subfamilies so far included within pleurostict scarabs except Orphninae. Our results represent a major step towards understanding the diversification history of one of the largest angiosperm-associated radiations of beetles.

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.

Species are the fundamental units of life and evolution. Their recognition is essential for science and society. Molecular methods have been increasingly used for the identification of animal species, despite several challenges. Here, we explore with genomic data from nine animal lineages a set of nuclear markers, namely metazoan‐level universal single‐copy orthologs (metazoan USCOs), for their use in species delimitation. Our data sets include arthropods and vertebrates. We use various data assembly strategies and use coalescent‐based species inference as well as population admixture analyses and phenetic methods. We demonstrate that metazoan USCOs distinguish well closely related morphospecies and consistently outperform classical mitochondrial DNA barcoding in discriminating closely related species in different animal taxa, as judged by comparison with morphospecies delimitations. USCOs overcome the general shortcomings of mitochondrial DNA barcodes, and due to standardization across Metazoa, also those of other approaches. They accurately assign samples not only to lower but also to higher taxonomic levels. Metazoan USCOs provide a powerful and unifying framework for DNA‐based species delimitation and taxonomy in animals and their employment could result in a more efficient use of research data and resources.

AbstractGeometric morphometrics (GM) enable the quantification of morphological variation on various scales. Recent technical advances allow analyzing complex three‐dimensional shapes also in cases where landmark‐based approaches are not appropriate. Pelvic girdle bones (basipterygia) of Sulawesi ricefishes are 3D structures that challenge traditional morphometrics. We hypothesize that the pelvic girdle of ricefishes experienced sex‐biased selection pressures in species where females provide brood care by carrying fertilized eggs supported by elongated pelvic fins (“pelvic brooding”). We test this by comparing pelvic bone shapes of both sexes in species exhibiting pelvic brooding and the more common reproductive strategy “transfer brooding,” by using landmark‐free 2D and 3D GM, as well as qualitative shape descriptions. Both landmark‐free approaches revealed significant interspecific pelvic bone variation in the lateral process, medial facing side of the pelvic bone, and overall external and internal wing shape. Within pelvic brooders, the three analyzed species are clearly distinct, while pelvic bones of the genus Adrianichthys are more similar to transfer brooding Oryzias. Female pelvic brooding Oryzias exhibit prominent, medially pointing tips extending from the internal wing and basipterygial plate that are reduced or absent in conspecific males, Adrianichthys and transfer brooding Oryzias, supporting our hypothesis that selection pressures affecting pelvic girdle shape are sex‐biased in Sulawesi ricefishes. Furthermore, both sexes of pelvic brooding Oryzias have overall larger pelvic bones than other investigated ricefishes. Based on these differences, we characterized two reproductive strategy‐ and sex‐dependent pelvic girdle types for Sulawesi ricefishes. Morphological differences between the investigated pelvic brooding genera Adrianichthys and Oryzias provide additional evidence for two independent origins of pelvic brooding. Overall, our findings add to a better understanding on traits related to pelvic brooding in ricefishes and provide a basis for upcoming studies on pelvic girdle function and morphology.

The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA–BABA test with D -statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations.

Cuticular hydrocarbons (CHCs) cover the cuticle of insects and serve as desiccation barrier and as semiochemicals. While the main enzymatic steps of CHC biosynthesis are well understood, few of the underlying genes have been identified. Here we show how exploitation of intrasexual CHC dimorphism in a mason wasp, Odynerus spinipes, in combination with whole-genome sequencing and comparative transcriptomics facilitated identification of such genes. RNAi-mediated knockdown of twelve candidate gene orthologs in the honey bee, Apis mellifera, confirmed nine genes impacting CHC profile composition. Most of them have predicted functions consistent with current knowledge of CHC metabolism. However, we found first-time evidence for a fatty acid amide hydrolase also influencing CHC profile composition. In situ hybridization experiments furthermore suggest trophocytes participating in CHC biosynthesis. Our results set the base for experimental CHC profile manipulation in Hymenoptera and imply that the evolutionary origin of CHC biosynthesis predates the arthropods' colonization of land.

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.

The evolution of complex phenotypes like reproductive strategies is challenging to understand, as they often depend on multiple adaptations that only jointly result in a specific functionality. Sulawesi ricefishes (Adrianichthyidae) evolved a reproductive strategy termed as pelvic brooding. In contrast to the more common transfer brooding, female pelvic brooders carry an egg bundle connected to their body for weeks until the fry hatches. To examine the genetic architecture of pelvic brooding, we crossed the pelvic brooding Oryzias eversi and the transfer brooding Oryzias nigrimas (species divergence time: ∼3.6 my). We hypothesize, that a low number of loci and modularity have facilitated the rapid evolution of pelvic brooding. Traits associated to pelvic brooding, like rib length, pelvic fin length, and morphology of the genital papilla, were correlated in the parental species but correlations were reduced or lost in their F1 and F2 hybrids. Using the Castle-Wright estimator, we found that generally few loci underlie the studied traits. Further, both parental species showed modularity in their body plans. In conclusion, morphological traits related to pelvic brooding were based on a few loci and the mid-body region likely could evolve independently from the remaining body parts. Both factors presumably facilitated the evolution of pelvic brooding.

Extant members of the ancient insect order of stoneflies exhibit a disjunct, antitropical distribution, with one major lineage exclusively occurring in the Southern Hemisphere and the other, with few exceptions, on the Northern continents. Here, we address the biogeographic distribution and phylogenetic relationships of stoneflies using a phylogenetic workflow that combines both transcriptomic and Sanger sequence datasets with heterogeneous taxon coverage. We used a dataset comprising 2997 genes derived from the transcriptomes of 30 species and Sanger sequences of seven genes for 498 species. The backbone phylogeny was mainly inferred from the transcriptomic data, whereas the Sanger nucleotide sequence data provided high species density for divergence time estimation and diversification analyses. Our results show that the biogeographic pattern we observe today is primarily more likely shaped by long‐distance over‐land dispersal than by vicariance. We inferred that the ancestors of extant stoneflies originated in the Northern Hemisphere approximately 265 Ma and were presumably restricted to this area due to climatic and geographic boundaries. Our analyses suggest that with the break‐up of Pangaea around 200 Ma and the associated climatic and geographical changes, two groups of stoneflies, the Anarctoperlaria and the Notonemouridae, dispersed to Gondwana and subsequently went extinct on the northern continents. Both groups likely dispersed across Gondwana before its break‐up into the modern continents. At least one member of another group of ‘northern’ stoneflies, the Acroneuriinae, seems to have migrated from North America to South America around 67 Ma. We found four major net diversification rate shifts, indicating rapid radiation patterns that hampered a robust phylogenetic placement of these stonefly groups. Our study provides the first conclusive evolutionary explanation for the unique distribution pattern of stoneflies.