Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges.

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

Allelic resolution of insect and spider silk genes reveals hidden genetic diversity.

Arthropod silk is vital to the evolutionary success of hundreds of thousands of species. The primary proteins in silks are often encoded by long, repetitive gene sequences. Until recently, sequencing and assembling these complex gene sequences has proven intractable given their repetitive structure. Here, using high-quality long-read sequencing, we show that there is extensive variation-both in terms of length and repeat motif order-between alleles of silk genes within individual arthropods. Further, this variation exists across two deep, independent origins of silk which diverged more than 500 Mya: the insect clade containing caddisflies and butterflies and spiders. This remarkable convergence in previously overlooked patterns of allelic variation across multiple origins of silk suggests common mechanisms for the generation and maintenance of structural protein-coding genes. Future genomic efforts to connect genotypes to phenotypes should account for such allelic variation.

Phylogenomics recovers multiple origins of portable case making in caddisflies (Insecta: Trichoptera), nature's underwater architects.

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.

The de novo genome of the Black-necked Snakefly (Venustoraphidia nigricollis Albarda, 1891): A resource to study the evolution of living fossils.

Snakeflies (Raphidioptera) are the smallest order of holometabolous insects that have kept their distinct and name-giving appearance since the Mesozoic, probably since the Jurassic, and possibly even since their emergence in the Carboniferous, more than 300 million years ago. Despite their interesting nature and numerous publications on their morphology, taxonomy, systematics, and biogeography, snakeflies have never received much attention from the general public, and only a few studies were devoted to their molecular biology. Due to this lack of molecular data, it is therefore unknown, if the conserved morphological nature of these living fossils translates to conserved genomic structures. Here, we present the first genome of the species and of the entire order of Raphidioptera. The final genome assembly has a total length of 669 Mbp and reached a high continuity with an N50 of 5.07 Mbp. Further quality controls also indicate a high completeness and no meaningful contamination. The newly generated data was used in a large-scaled phylogenetic analysis of snakeflies using shared orthologous sequences. Quartet score and gene concordance analyses revealed high amounts of conflicting signals within this group that might speak for substantial incomplete lineage sorting and introgression after their presumed re-radiation after the asteroid impact 66 million years ago. Overall, this reference genome will be a door-opening dataset for many future research applications, and we demonstrated its utility in a phylogenetic analysis that provides new insights into the evolution of this group of living fossils.

Adaptive monitoring in action-what drives arthropod diversity and composition in central European beech forests?

Recent studies suggest that arthropod diversity in German forests is declining. Currently, different national programs are being developed to monitor arthropod trends and to unravel the effects of forest management on biodiversity in forests. To establish effective long-term monitoring programs, a set of drivers of arthropod diversity and composition as well as suitable species groups have to be identified. To aid in answering these questions, we investigated arthropod data collected in four Hessian forest reserves (FR) in the 1990s. To fully utilize this data set, we combined it with results from a retrospective structural sampling design applied at the original trap locations in central European beech (Fagus sylvatica) forests. As expected, the importance of the different forest structural, vegetation, and site attributes differed largely between the investigated arthropod groups: beetles, spiders, Aculeata, and true bugs. Measures related to light availability and temperature such as canopy cover or potential radiation were important to all groups affecting either richness, composition, or both. Spiders and true bugs were affected by the broadest range of explanatory variables, which makes them a good choice for monitoring general trends. For targeted monitoring focused on forestry-related effects on biodiversity, rove and ground beetles seem more suitable. Both groups were driven by a narrower, more management-related set of variables. Most importantly, our study approach shows that it is possible to utilize older biodiversity survey data. Although, in our case, there are strong restrictions due to the long time between species and structural attribute sampling.

Conservation of Three-Dimensional Structure of Lepidoptera and Trichoptera L-Fibroins for 290 Million Years.

The divergence of sister orders Trichoptera (caddisflies) and Lepidoptera (moths and butterflies) from a silk-spinning ancestor occurred around 290 million years ago. Trichoptera larvae are mainly aquatic, and Lepidoptera larvae are almost entirely terrestrial-distinct habitats that required molecular adaptation of their silk for deployment in water and air, respectively. The major protein components of their silks are heavy chain and light chain fibroins. In an effort to identify molecular changes in L-fibroins that may have contributed to the divergent use of silk in water and air, we used the ColabFold implementation of AlphaFold2 to predict three-dimensional structures of L-fibroins from both orders. A comparison of the structures revealed that despite the ancient divergence, profoundly different habitats, and low sequence conservation, a novel 10-helix core structure was strongly conserved in L-fibroins from both orders. Previously known intra- and intermolecular disulfide linkages were accurately predicted. Structural variations outside of the core may represent molecular changes that contributed to the evolution of insect silks adapted to water or air. The distributions of electrostatic potential, for example, were not conserved and present distinct order-specific surfaces for potential interactions with or modulation by external factors. Additionally, the interactions of L-fibroins with the H-fibroin C-termini are different for these orders; lepidopteran L-fibroins have N-terminal insertions that are not present in trichopteran L-fibroins, which form an unstructured ribbon in isolation but become part of an intermolecular ß-sheet when folded with their corresponding H-fibroin C-termini. The results are an example of protein structure prediction from deep sequence data of understudied proteins made possible by AlphaFold2.

Integrative taxonomy by molecular species delimitation: multi-locus data corroborate a new species of Balkan Drusinae micro-endemics.

BACKGROUND: Taxonomy offers precise species identification and delimitation and thus provides basic information for biological research, e.g. through assessment of species richness. The importance of molecular taxonomy, i.e., the identification and delimitation of taxa based on molecular markers, has increased in the past decade. Recently developed exploratory tools now allow estimating species-level diversity in multi-locus molecular datasets. RESULTS: Here we use molecular species delimitation tools that either quantify differences in intra- and interspecific variability of loci, or divergence times within and between species, or perform coalescent species tree inference to estimate species-level entities in molecular genetic datasets. We benchmark results from these methods against 14 morphologically readily differentiable species of a well-defined subgroup of the diverse Drusinae subfamily (Trichoptera, Limnephilidae). Using a 3798 bp (6 loci) molecular data set we aim to corroborate a geographically isolated new species by integrating comparative morphological studies and molecular taxonomy. CONCLUSIONS: Our results indicate that only multi-locus species delimitation provides taxonomically relevant information. The data further corroborate the new species Drusus zivici sp. nov. We provide differential diagnostic characters and describe the male, female and larva of this new species and discuss diversity patterns of Drusinae in the Balkans. We further discuss potential and significance of molecular species delimitation. Finally we argue that enhancing collaborative integrative taxonomy will accelerate assessment of global diversity and completion of reference libraries for applied fields, e.g., conservation and biomonitoring.

A hairy case: The evolution of filtering carnivorous Drusinae (Limnephilidae, Trichoptera).

The caddisfly subfamily Drusinae BANKS comprises roughly 100 species inhabiting mountain ranges in Europe, Asia Minor and the Caucasus. A 3-gene phylogeny of the subfamily previously identified three major clades that were corroborated by larval morphology and feeding ecologies: scraping grazers, omnivorous shredders and filtering carnivores. Larvae of filtering carnivores exhibit unique head capsule complexities, unknown from other caddisfly larvae. Here we assess the species-level relationships within filtering carnivores, hypothesizing that head capsule complexity is derived from simple shapes observed in the other feeding groups. We summarize the current systematics and taxonomy of the group, clarify the systematic position of Cryptothrix nebulicola, and present a larval key to filtering carnivorous Drusinae. We infer relationships of all known filtering carnivorous Drusinae and 34 additional Drusinae species using Bayesian species tree analysis and concatenated Bayesian phylogenetic analysis of 3805bp of sequence data from six gene regions (mtCOI5-P, mtCOI3-P, 16S mrDNA, CADH, WG, 28S nrDNA), morphological cladistics from 308 characters, and a total evidence analysis. All analyses support monophyly of the three feeding ecology groups but fail to fully resolve internal relationships. Within filtering carnivores, variation in head setation and frontoclypeus structure may be associated with progressive niche adaptation, with less complex species recovered at a basal position. We propose that diversification of complex setation and frontoclypeus shape represents a recent evolutionary development, hypothetically enforcing speciation and niche specificity within filtering carnivorous Drusinae.

Larval morphology and phylogenetic position of Drusus balcanicus, Drusus botosaneanui, Drusus serbicus and Drusus tenellus (Trichoptera: Limnephilidae: Drusinae).

In a recent 3-gene phylogeny of the Trichoptera subfamily Drusinae Banks, 1916 molecular data clearly correlated with the morphology and feeding ecology of larvae. The largest of three main groups, the Drusinae grazer clade, exhibits an unusual larval feeding ecology for Limnephilidae, and is the most diverse group. In this paper we describe four previously unknown Drusinae larvae from this clade: Drusus balcanicus Kumanski, 1973 (micro-endemic to Eastern Balkans); Drusus botosaneanui Kumanski, 1968 (Dinaric Western Balkans, Hellenic and Eastern Balkan, Asia Minor), Drusus serbicus Marinkovic-Gospodnetic, 1971a (micro-endemic to Dinaric Western Balkans); and Drusus tenellus (Klapálek, 1898) (Carpathians, Dinaric Eastern Balkans). Characteristically, the larvae of these species develop toothless mandibles typical for the Drusinae grazer clade. Larvae and adults were unambiguously associated by a phylogenetic approach based on two mitochondrial (mtCOI, mtLSU= 16S rDNA) and two nuclear genes (nuWG, nuCAD). In addition, information on the morphology of the larvae is given and the diagnostic features necessary for identification are illustrated.

Freshwater biodiversity and aquatic insect diversification.

Inland waters cover less than 1% of Earth's surface but harbor more than 6% of all insect species: Nearly 100,000 species from 12 orders spend one or more life stages in freshwater. Little is known about how this remarkable diversity arose, although allopatric speciation and ecological adaptation are thought to be primary mechanisms. Freshwater habitats are highly susceptible to environmental change and exhibit marked ecological gradients. Standing waters appear to harbor more dispersive species than running waters, but there is little understanding of how this fundamental ecological difference has affected diversification. In contrast to the lack of evolutionary studies, the ecology and habitat preferences of aquatic insects have been intensively studied, in part because of their widespread use as bioindicators. The combination of phylogenetics with the extensive ecological data provides a promising avenue for future research, making aquatic insects highly suitable models for the study of ecological diversification.

Genomic resources of two aquatic Lepidoptera, Elophila obliteralis and Hyposmocoma kahamanoa, reveal similarities with Trichoptera in amino acid composition of major silk genes.

While most species of butterflies and moths (Lepidoptera) have entirely terrestrial life histories, ∼0.5% of the described species are known to have an aquatic larval stage. Larvae of aquatic Lepidoptera are similar to caddisflies (Trichoptera) in that they use silk to anchor themselves to underwater substrates or to build protective cases. However, the physical properties and genetic elements of silks in aquatic Lepidoptera remain unstudied, as most research on lepidopteran silk has focused on the commercially important silkworm, Bombyx mori. Here, we provide high-quality PacBio HiFi genome assemblies of two distantly-related aquatic Lepidoptera species (Elophila obliteralis (Pyraloidea: Crambidae) and Hyposmocoma kahamanoa (Gelechioidea: Cosmopterigidae)). As a step toward understanding the evolution of underwater silk in aquatic Lepidoptera, we used our two genome assemblies and compared them to published genetic data of aquatic and terrestrial Lepidoptera. Sequences of the primary silk protein, h-fibroin in aquatic moths have conserved termini and share a basic motif structure with terrestrial Lepidoptera. However, these sequences were similar to aquatic Trichoptera in that the percentage of positively and negatively charged amino acids was much higher than in terrestrial Lepidoptera, indicating a possible adaptation of silks to aquatic environments.

The impact of global climate change on genetic diversity within populations and species.

Genetic diversity provides the basic substrate for evolution, yet few studies assess the impacts of global climate change (GCC) on intraspecific genetic variation. In this review, we highlight the importance of incorporating neutral and non-neutral genetic diversity when assessing the impacts of GCC, for example, in studies that aim to predict the future distribution and fate of a species or ecological community. Specifically, we address the following questions: Why study the effects of GCC on intraspecific genetic diversity? How does GCC affect genetic diversity? How is the effect of GCC on genetic diversity currently studied? Where is potential for future research? For each of these questions, we provide a general background and highlight case studies across the animal, plant and microbial kingdoms. We further discuss how cryptic diversity can affect GCC assessments, how genetic diversity can be integrated into studies that aim to predict species' responses on GCC and how conservation efforts related to GCC can incorporate and profit from inclusion of genetic diversity assessments. We argue that studying the fate of intraspecifc genetic diversity is an indispensable and logical venture if we are to fully understand the consequences of GCC on biodiversity on all levels.

Population genetics of ecological communities with DNA barcodes: an example from New Guinea Lepidoptera.

Comparative population genetics of ecological guilds can reveal generalities in patterns of differentiation bearing on hypotheses regarding the origin and maintenance of community diversity. Contradictory estimates of host specificity and beta diversity in tropical Lepidoptera (moths and butterflies) from New Guinea and the Americas have sparked debate on the role of host-associated divergence and geographic isolation in explaining latitudinal diversity gradients. We sampled haplotypes of mitochondrial cytochrome c oxidase I from 28 Lepidoptera species and 1,359 individuals across four host plant genera and eight sites in New Guinea to estimate population divergence in relation to host specificity and geography. Analyses of molecular variance and haplotype networks indicate varying patterns of genetic structure among ecologically similar sympatric species. One-quarter lacked evidence of isolation by distance or host-associated differentiation, whereas 21% exhibited both. Fourteen percent of the species exhibited host-associated differentiation without geographic isolation, 18% showed the opposite, and 21% were equivocal, insofar as analyses of molecular variance and haplotype networks yielded incongruent patterns. Variation in dietary breadth among community members suggests that speciation by specialization is an important, but not universal, mechanism for diversification of tropical Lepidoptera. Geographically widespread haplotypes challenge predictions of vicariance biogeography. Dispersal is important, and Lepidoptera communities appear to be highly dynamic according to the various phylogeographic histories of component species. Population genetic comparisons among herbivores of major tropical and temperate regions are needed to test predictions of ecological theory and evaluate global patterns of biodiversity.

The larvae of Drusus franzressli Malicky 1974 and Drusus spelaeus (Ulmer 1920) (Trichoptera: Limnephilidae: Drusinae) with notes on ecology and zoogeography.

Water quality monitoring is greatly dependent on identification tools for aquatic and semi-aquatic insects. Species-level identification improves resolution and precision of water quality assessment and requires comprehensive keys. With the aim of increasing the suitability of Drusinae for such applications, this paper gives a description of the hitherto unknown larvae of Drusus franzressli Malicky 1974 and Drusus spelaeus (Ulmer 1920). Information on the morphology of the larvae is given and the most important diagnostic features are illustrated. In the context of already available keys, the larvae of D. franzressli and D. spelaeus key together with Metanoea flavipennis (Pictet 1834), M. rhaetica Schmid 1956, D. improvisus McLachlan 1884, D. nigrescens Meyer-Dür 1875 and Ecclisopteryx malickyi Moretti 1991. These species are easily separated by differences in larval morphology (dorsal outline and sculpturing of pronotum, presence/absence of lateral gills at 2nd and 3rd abdominal segments, start of lateral fringe) and their distribution ranges. Drusus franzressli is endemic to the Hellenic western Balkans whereas D. spelaeus is endemic to the western Alps (Grenoble area). In addition, ecological characteristics are briefly discussed.

Holotype sequencing of Silvataresholzenthali Rázuri-Gonzales, Ngera & Pauls, 2022 (Trichoptera, Pisuliidae).

While DNA barcodes are increasingly provided in descriptions of new species, the whole mitochondrial and nuclear genomes are still rarely included. This is unfortunate because whole genome sequencing of holotypes allows perpetual genetic characterization of the most representative specimen for a given species. Thus, de novo genomes are invaluable additional diagnostic characters in species descriptions, provided the structural integrity of the holotype specimens remains intact. Here, we used a minimally invasive method to extract DNA of the type specimen of the recently described caddisfly species Silvataresholzenthali Rázuri-Gonzales, Ngera & Pauls, 2022 (Trichoptera: Pisuliidae) from the Democratic Republic of the Congo. A low-cost next generation sequencing strategy was used to generate the complete mitochondrial and draft nuclear genome of the holotype. The data in its current form is an important extension to the morphological species description and valuable for phylogenomic studies.

Characterization of the primary structure of the major silk gene, h-fibroin, across caddisfly (Trichoptera) suborders.

Larvae of caddisflies (Trichoptera) produce silk to build various underwater structures allowing them to exploit a wide range of aquatic environments. The silk adheres to various substrates underwater and has high tensile strength, extensibility, and toughness and is of interest as a model for biomimetic adhesives. As a step toward understanding how the properties of underwater silk evolved in Trichoptera, we used genomic data to identify full-length sequences and characterize the primary structure of the major silk protein, h-fibroin, across the order. The h-fibroins have conserved termini and basic motif structure with high variation in repeating modules and variation in the percentage of amino acids, mainly proline. This finding might be linked to differences in mechanical properties related to the different silk usage and sets a starting point for future studies to screen and correlate amino acid motifs and other sequence features with quantifiable silk properties.

RNA interference to combat the Asian tiger mosquito in Europe: A pathway from design of an innovative vector control tool to its application.

The Asian tiger mosquito Aedes albopictus is currently spreading across Europe, facilitated by climate change and global transportation. It is a vector of arboviruses causing human diseases such as chikungunya, dengue hemorrhagic fever and Zika fever. For the majority of these diseases, no vaccines or therapeutics are available. Options for the control of Ae. albopictus are limited by European regulations introduced to protect biodiversity by restricting or phasing out the use of pesticides, genetically modified organisms (GMOs) or products of genome editing. Alternative solutions are thus urgently needed to avoid a future scenario in which Europe faces a choice between prioritizing human health or biodiversity when it comes to Aedes-vectored pathogens. To ensure regulatory compliance and public acceptance, these solutions should preferably not be based on chemicals or GMOs and must be cost-efficient and specific. The present review aims to synthesize available evidence on RNAi-based mosquito vector control and its potential for application in the European Union. The recent literature has identified some potential target sites in Ae. albopictus and formulations for delivery. However, we found little information concerning non-target effects on the environment or human health, on social aspects, regulatory frameworks, or on management perspectives. We propose optimal designs for RNAi-based vector control tools against Ae. albopictus (target product profiles), discuss their efficacy and reflect on potential risks to environmental health and the importance of societal aspects. The roadmap from design to application will provide readers with a comprehensive perspective on the application of emerging RNAi-based vector control tools for the suppression of Ae. albopictus populations with special focus on Europe.

Methodological framework for projecting the potential loss of intraspecific genetic diversity due to global climate change.

BACKGROUND: While research on the impact of global climate change (GCC) on ecosystems and species is flourishing, a fundamental component of biodiversity - molecular variation - has not yet received its due attention in such studies. Here we present a methodological framework for projecting the loss of intraspecific genetic diversity due to GCC. METHODS: The framework consists of multiple steps that combines 1) hierarchical genetic clustering methods to define comparable units of inference, 2) species accumulation curves (SAC) to infer sampling completeness, and 3) species distribution modelling (SDM) to project the genetic diversity loss under GCC. We suggest procedures for existing data sets as well as specifically designed studies. We illustrate the approach with two worked examples from a land snail (Trochulus villosus) and a caddisfly (Smicridea (S.) mucronata). RESULTS: Sampling completeness was diagnosed on the third coarsest haplotype clade level for T. villosus and the second coarsest for S. mucronata. For both species, a substantial species range loss was projected under the chosen climate scenario. However, despite substantial differences in data set quality concerning spatial sampling and sampling depth, no loss of haplotype clades due to GCC was predicted for either species. CONCLUSIONS: The suggested approach presents a feasible method to tap the rich resources of existing phylogeographic data sets and guide the design and analysis of studies explicitly designed to estimate the impact of GCC on a currently still neglected level of biodiversity.

A new species of Silvatares (Trichoptera, Pisuliidae) from the Democratic Republic of the Congo.

A new species of caddisfly in the family Pisuliidae from the Democratic Republic of the Congo is described and illustrated herein, Silvataresholzenthali sp. nov. Based on the presence of a pair of spines on the endotheca, this species belongs to the thrymmifer group. Additionally, Silvatareslaetae is recorded for the first time from the D.R. Congo.

De Novo Genome Assembly and Annotation of an Andean Caddisfly, Atopsyche davidsoni Sykora, 1991, a Model for Genome Research of High-Elevation Adaptations.

We sequence, assemble, and annotate the genome of Atopsyche davidsoni Sykora, 1991, the first whole-genome assembly for the caddisfly family Hydrobiosidae. This free-living and predatory caddisfly inhabits streams in the high-elevation Andes and is separated by more than 200 Myr of evolutionary history from the most closely related caddisfly species with genome assemblies available. We demonstrate the promise of PacBio HiFi reads by assembling the most contiguous caddisfly genome assembly to date with a contig N50 of 14 Mb, which is more than 6× more contiguous than the current most contiguous assembly for a caddisfly (Hydropsyche tenuis). We recover 98.8% of insect BUSCO genes indicating a high level of gene completeness. We also provide a genome annotation of 12,232 annotated proteins. This new genome assembly provides an important new resource for studying genomic adaptation of aquatic insects to harsh, high-altitude environments.