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.

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.

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.

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.

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.

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.

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.

The impact of sequencing depth and relatedness of the reference genome in population genomic studies: A case study with two caddisfly species (Trichoptera, Rhyacophilidae, Himalopsyche).

Whole genome sequencing for generating SNP data is increasingly used in population genetic studies. However, obtaining genomes for massive numbers of samples is still not within the budgets of many researchers. It is thus imperative to select an appropriate reference genome and sequencing depth to ensure the accuracy of the results for a specific research question, while balancing cost and feasibility. To evaluate the effect of the choice of the reference genome and sequencing depth on downstream analyses, we used five confamilial reference genomes of variable relatedness and three levels of sequencing depth (3.5×, 7.5× and 12×) in a population genomic study on two caddisfly species: Himalopsyche digitata and H. tibetana. Using these 30 datasets (five reference genomes × three depths × two target species), we estimated population genetic indices (inbreeding coefficient, nucleotide diversity, pairwise F ST, and genome-wide distribution of F ST) based on variants and population structure (PCA and admixture) based on genotype likelihood estimates. The results showed that both distantly related reference genomes and lower sequencing depth lead to degradation of resolution. In addition, choosing a more closely related reference genome may significantly remedy the defects caused by low depth. Therefore, we conclude that population genetic studies would benefit from closely related reference genomes, especially as the costs of obtaining a high-quality reference genome continue to decrease. However, to determine a cost-efficient strategy for a specific population genomic study, a trade-off between reference genome relatedness and sequencing depth can be considered.

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.

Genome size evolution in the diverse insect order Trichoptera.

BACKGROUND: Genome size is implicated in the form, function, and ecological success of a species. Two principally different mechanisms are proposed as major drivers of eukaryotic genome evolution and diversity: polyploidy (i.e., whole-genome duplication) or smaller duplication events and bursts in the activity of repetitive elements. Here, we generated de novo genome assemblies of 17 caddisflies covering all major lineages of Trichoptera. Using these and previously sequenced genomes, we use caddisflies as a model for understanding genome size evolution in diverse insect lineages. RESULTS: We detect a ∼14-fold variation in genome size across the order Trichoptera. We find strong evidence that repetitive element expansions, particularly those of transposable elements (TEs), are important drivers of large caddisfly genome sizes. Using an innovative method to examine TEs associated with universal single-copy orthologs (i.e., BUSCO genes), we find that TE expansions have a major impact on protein-coding gene regions, with TE-gene associations showing a linear relationship with increasing genome size. Intriguingly, we find that expanded genomes preferentially evolved in caddisfly clades with a higher ecological diversity (i.e., various feeding modes, diversification in variable, less stable environments). CONCLUSION: Our findings provide a platform to test hypotheses about the potential evolutionary roles of TE activity and TE-gene associations, particularly in groups with high species, ecological, and functional diversities.

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.

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.

Long-read HiFi sequencing correctly assembles repetitive heavy fibroin silk genes in new moth and caddisfly genomes.

Insect silk is a versatile biomaterial. Lepidoptera and Trichoptera display some of the most diverse uses of silk, with varying strength, adhesive qualities, and elastic properties. Silk fibroin genes are long (>20 Kbp), with many repetitive motifs that make them challenging to sequence. Most research thus far has focused on conserved N- and C-terminal regions of fibroin genes because a full comparison of repetitive regions across taxa has not been possible. Using the PacBio Sequel II system and SMRT sequencing, we generated high fidelity (HiFi) long-read genomic and transcriptomic sequences for the Indianmeal moth (Plodia interpunctella) and genomic sequences for the caddisfly Eubasilissa regina. Both genomes were highly contiguous (N50  = 9.7 Mbp/32.4 Mbp, L50  = 13/11) and complete (BUSCO complete  = 99.3%/95.2%), with complete and contiguous recovery of silk heavy fibroin gene sequences. We show that HiFi long-read sequencing is helpful for understanding genes with long, repetitive regions.

Gene Flow and Diversification in Himalopsyche martynovi Species Complex (Trichoptera: Rhyacophilidae) in the Hengduan Mountains.

The Hengduan Mountains are one of the most species-rich mountainous areas in the world. The origin and evolution of such a remarkable biodiversity are likely to be associated with geological or climatic dynamics, as well as taxon-specific biotic processes (e.g., hybridization, polyploidization, etc.). Here, we investigate the mechanisms fostering the diversification of the endemic Himalopsyche martynovi complex, a poorly known group of aquatic insects. We used multiple allelic datasets generated from 691 AHE loci to reconstruct species and RaxML phylogenetic trees. We selected the most reliable phylogenetic tree to perform network and gene flow analyses. The phylogenetic reconstructions and network analysis identified three clades, including H. epikur, H. martynovi sensu stricto and H. cf. martynovi. Himalopsyche martynovi sensu stricto and H. cf. martynovi present an intermediate morphology between H. epikur and H. viteceki, the closest known relative to the H. martynovi-complex. The gene flow analysis revealed extensive gene flow among these lineages. Our results suggest that H. viteceki and H. epikur are likely to have contributed to the evolution of H. martynovi sensu stricto and H. cf. martynovi via gene flow, and thus, our study provides insights in the diversification process of a lesser-known ecological group, and hints at the potential role of gene flow in the emergence of biological novelty in the Hengduan Mountains.

Long Reads Are Revolutionizing 20 Years of Insect Genome Sequencing.

The first insect genome assembly (Drosophila melanogaster) was published two decades ago. Today, nuclear genome assemblies are available for a staggering 601 insect species representing 20 orders. In this study, we analyzed the most-contiguous assembly for each species and provide a "state-of-the-field" perspective, emphasizing taxonomic representation, assembly quality, gene completeness, and sequencing technologies. Relative to species richness, genomic efforts have been biased toward four orders (Diptera, Hymenoptera, Collembola, and Phasmatodea), Coleoptera are underrepresented, and 11 orders still lack a publicly available genome assembly. The average insect genome assembly is 439.2 Mb in length with 87.5% of single-copy benchmarking genes intact. Most notable has been the impact of long-read sequencing; assemblies that incorporate long reads are ∼48× more contiguous than those that do not. We offer four recommendations as we collectively continue building insect genome resources: 1) seek better integration between independent research groups and consortia, 2) balance future sampling between filling taxonomic gaps and generating data for targeted questions, 3) take advantage of long-read sequencing technologies, and 4) expand and improve gene annotations.

Potamophylax coronavirus sp. n. (Trichoptera: Limnephilidae), a new species from Bjeshkët e Nemuna National Park in the Republic of Kosovo, with molecular and ecological notes.

BACKGROUND: The Western Balkans are an important hotspot of caddisfly diversity in Europe, with several microscale endemics, many of which were discovered during the recent years. The genus Potamophylax Wallengren, 1891 likely originated and diversified in Europe, with the Balkan Peninsula being one of the most important diversity hotspots. NEW INFORMATION: In this paper, we describe the new species Potamophylax coronavirus sp. n. from Bjeshkët e Nemuna National Park in the Republic of Kosovo. The new species belongs to the Potamophylax winneguthi species group and is morphologically most similar to Potamophylax juliani Kumanski, 1999, currently known only from Bulgaria and Potamophylax winneguthi Klapalek, 1902, known from Bosnia and Herzegovina and Serbia. The male of the new species differs from its most similar congeners mainly in exhibiting: (1) elongated subrectangular superior appendages in lateral view; (2) hardly acuminate, almost rounded apex of intermediate appendages; (3) differently shaped, irregular and higher inferior appendages; (4) narrow spinate area, roughly rectangular in dorsal view, slightly wider at the base and (5) different paramere shape and/or spine pattern. The new species also differs by its considerably smaller size and association with open, high altitude eucrenal zones.The uncorrected interspecific pairwise distance between P. coronavirus and other species of the P. winneguthi species group is on par with those amongst other recognised species in the group, as well as with the yet ambiguously identified taxa from the Sharr, Rila and Bajgorë Mountains.The new species is most probably a microendemic of Bjeshkët e Nemuna, thus highlighting further this area as an important hotspot of caddisfly biodiversity in Europe.

A new genus and species of mud snake from Myanmar (Reptilia, Squamata, Homalopsidae).

Based on two male and two female individuals, we describe a new genus and species of mud snake, Myanophis thanlyinensis gen. nov., sp. nov., from the vicinity of the campus of East Yangon University, Yangon, Thanlyin, Myanmar. This species differs from every other homalopsid species by the following combination of characters: (1) dorsal scales smooth, row formula 21-21-19 or 21-21-17; (2) tail short, ratio tail length/SVL 0.185-0.204 in males, 0.160-0.167 in females; (3) nasal scales separated; (4) 125-126 ventral scales in males, 120-122 in females; (5) 38-39 subcaudal scales in males, 32-34 in females; and (6) hemipenis bilobed. Its matrilineal genealogy (based on analyses of 16S and cytochrome b sequences), associates Myanophis thanlyinensis gen. nov., sp. nov. most closely with species of the genera Myrrophis and Gyiophis. The new taxon differs from the species of Myrrophis and Gyiophis by having a bilobed hemipenis (vs. unilobed). Myanophis thanlyinensis gen. nov., sp. nov. differs further from the species of Myrrophis by having 125-126 ventral scales in males and 120-122 in females (vs. 137-162 and 137-164, respectively), and 38-39 subcaudal scales in males and 32-34 in females (vs. 39-55 and 37-52, respectively). Myanophis thanlyinensis gen. nov., sp. nov. differs further from the species of Gyiophis by lacking dark blotches along flank (vs. present), and by having 21 dorsal scales rows at midbody (vs. 25). We provide an identification key to the homalopsid species known to occur in Myanmar. As a novelty to the classic holotype description and characterization, the individual has been genome sequenced by Illumina short-read technology and its genome has been assembled into a draft nuclear genome and a complete, annotated mitochondrial genome. This innovative approach comprehensively and permanently characterizes the genomic variation of the holotype.